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According to the simple theory of the rigid rotator and 
harmonic oscillator, the coefficients of J(J+1) and (v+}) in 
the energy should be inversely proportional to respectively 
the reduced mass and the square root thereof. Actually 
various small corrections, important only in compounds 
containing hydrogen or deuterium, make the mass depend- 
ence slightly different, and hence, at least in princi- 
ple, explain why spectroscopic observations may yield 
slightly incorrect mass ratios for isotopes. The cor- 
rections are due to four causes (a) anharmonicity, (b) 
interaction between vibration and electronic motion even 
when perturbations between different electronic states are 
neglected, (c) Z uncoupling or in other words perturbations 
between states differing by one unit in A, (d) interaction 
between states of equal A. The formulas for the various 
corrections are collected—some old, some new or more 


exact than previously. The relative magnitudes of the dif- 
ferent effects are estimated. Incidentally, perturbations 
between states of equal A are studied in general. They 
should increase in intensity with the vibrational quantum 
number 2, all other things being equal, in contrast to the 
perturbations due to Z uncoupling, which increase with J. 
The electronic isotope shift, or displacement in the absolute 
value of the energy is also treated. Here the correction for 
the motion of the center of gravity of the molecule, usually 
neglected, is shown to be fully coordinate in importance 
with other terms. Shifts for ultraviolet lines of hydrogen 
are calculated with the Wang wave functions, and agree 
qualitatively with experiment. A concluding section is 
included on the mean square angular momentum of Hz, 
which enters in connection with diamagnetism and the 
Stern-Gerlach effect. 





HE observed energies of singlet states of 
diatomic molecules are ordinarily repre- 
sented by a formula 


W=W*¥+hw*(v+3)+Bt*J(J+1)+---, (1) 


where v and J, are respectively, the vibrational 
and rotational quantum numbers. The coef- 
ficients w*, B*, however, do not have exactly the 
ideal values 


wo=(1/2r)(a/M)', 


characteristic of the harmonic oscillator and 
rigid rotator respectively. The notation M is used 
for the reduced mass M,M2/(Mi+M2); ro is 
the equilibrium internuclear distance for an 


‘A preliminary account of the present paper was given 
at the St. Louis meeting of the Am. Phys. Soc., Phys. Rev. 
49, 417 (1936). 


By=h?/8x2Mry2 (2) 


irrotational molecule devoid of centrifugal dis- 
tortion, and a is identified with the coefficient of 
the first term of the usual expansion 


E,(r) = Wot parc(é?+a1 +ack+---) (3) 


of the radial potential energy. Here §=(r—10)/ro. 

Let primes and double primes refer to two 
different isotopes. Both theoretically and em- 
pirically, the ratios w*’/w*’ and (B*’/B*’)! 
deviate slightly from the simple value 


wo’ /wo”’ _ (By /Bo"’)'= (M"/M’)} (4) 


obtained from (2). There are many small cor- 
rections which cause B*, w* to differ from Bo, w, 
and which hence make the apparent isotope 
ratios, as revealed by band spectra, somewhat 
different from (4). The literature is in a rather 
confused state regarding their number and 
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nature. The purpose of the present paper is to 
summarize the various correction terms which 
have been given to date, and especially to present 
certain additional terms which appear to have 
previously been overlooked although they may 
be as large as the others. We must, however, not 
fail to mention that notable theoretical work on 
isotope corrections has recently been published 
by Dieke, and the approximations which he makes 
are adequate in the examples which he studies. 
To avoid undue complication, we throughout 
confine our attention to singlet states. 

The corrections are all sufficiently small that 
they may be considered additive. In other words, 
if (a) is some effect which, taken by itself, makes 
B* acquire a value By(1+y.) different from (2), 
and if there are other disturbing influences (b), 
(c), (d), the final formula is 


Bt=Bi(ityatrwtyetya) =Bo(it+y). (5) 


Similarly we can take 
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The apparent isotope ratio (B*’/B*’’)) in the 
rotational case should differ from (4) by a 
factor 1+ }7'—}y’’. In the vibrational isotope 
effect, the analogous relation is 


w® eo" = (M"/M")W(1+8' — 8"). 


The value of 6 is not in general the same as that 
of $y. Hence the correction factor to (4) is 
unlike for rotation and vibration. The reason 
that y’ differs from y’’, or 6’ from 6”, is that y, 
6 are functions of the nuclear masses. Since y, 6 
are themselves very small quantities, the ex- 
pressions y’—~y’’ and 6’—6” will be quite neg- 
ligible unless the ratio of the reduced masses 
for the two isotopes differs markedly from 
unity. Hence all the corrections which we discuss 
are important only when the pair of isotopes is 
a hydride and a deuteride. In fact, it is only the 
discovery of deuterium that makes most of the 
present paper of other than academic interest. 
We may distinguish between four effects (a), 
(b), (c), (d), all of which make the empirical 





isotope ratios differ from the ideal values (4). 


w* =wol1+oa.+5+5.+5a]=wolL1+s]. (6) 





(a) ANHARMONICITY 


Even for a vibrator-rotator, involving only nuclear in distinction from electronic coordinates, 
there will be departures from (4) if the restoring forces are not simple harmonic. The resulting 
modifications have been studied particularly by Dunham? and Crawford.* Dunham’s formulas show 
that 

Va = (Bo?/2wo?) [15 +410; —9a2+ 1543 — 23a ;a2+21($a;?+ 3a;') |, (7) 


6a = (Bo? / 329?) (20045 — 380a,a3— 134a2?+459a 2a), (8) 


with a;--+-a4 as in (3). 
(b) DIAGONAL ELEMENTS 
All perturbation calculations on diatomic molecules start with a set of initial wave functions ¥ 
in which the nuclear and electronic variables are separated, so that‘ 


++, 27, 7)Ravl(r)usa(O, o)f(X, Y, Z). (9) 


Here ®, is a function of the coordinates x,---z; of the electrons relative to the center of gravity as 
origin ; R,,, depends only on the internuclear distance 7, and u,, is a symmetrical top function which 
involves merely the angular variables 6, ¢ specifying the spatial orientation of the molecular axis. 
The function f depends only on the coordinates X, Y, Z of the center of gravity, and has the usual 
plane wave exponential form. The factor ®, contains 7 as a parameter, and satisfies the wave equation 


[H-E,(r) ]®, =0 (10) 


Vv =@®,(%x1, 


for ‘‘clamped nuclei.” The eigenvalue of (10) for a given state ” involves r parametrically, and is 
the same as the radial potential function (3). It is usually assumed that the differential equation 
2 J. L. Dunham, Phys. Rev. 41, 721 (1932). 


’F. H. Crawford and T. Jorgensen, Jr., Phys. Rev. 47, 931 (1935). 
*R.de L. Kronig, Band Spectra and Molecular Structure, p. 42; or J. H. Van Vleck, Phys. Rev. 33, 467 (1929). 
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satisfied by R,, is that of a system of one degree of freedom with the potential (3) supplemented by 
the ordinary centrifugal term [J(J+1)—A?*]h?/8x?Mr*. This supposition, however, is not quite 
right. Instead the wave equation appropriate to R,, is 


(—h2/8x2M9)(d2R n»/dr?) +[En(1) + Pan tQnntSnrt(J24+J—A)h?/822Mr? Ri, —WR»»=0, (11) 
where Parl’) =(h?/82r?Mr*)Enyn with EZyr=Sf-- -S%,'(P2+P,?)6,dr, (12) 
Onn(r) =(h?/82?M) Sf ++ > f —&,'0°@,,/dr'dr. (13) 


The formula for S,,, is furnished later in Eq. (14). d7 is an abbreviation for the 3f-dimensional volume 
element dv,---dv;. In (12), P, is the operator 7~'D,(y,0- - -/02,—2,0---/dy,) which gives the elec- 
tronic angular momentum about the x axis, measured in multiples of 4/27, while P, is the correspond- 
ing operator for the y axis. It is assumed throughout that the z axis connects the nuclei. The dagger 
is used to denote the conjugate imaginary, and for present purposes we take n’ =n. In writing (11), 
we have used a radial wave function normalized with respect to dr (rather than to r’dr as in reference 
4), for this procedure makes the formulas more compact. 

The extra terms Pn, Qnn in (11) may be qualitatively described as representing coupling between 
vibration and electronic motion. They have already been given by Kronig and others.4 The remaining 
additional term S,, owes its origin to the fact that the center of gravity of the molecule does not 
quite coincide with the center of gravity of the nuclei, so that the latter wobbles about the former. 
This term does not appear to have been explicitly written down previously in connection with the 
radial equation, because it has commonly been assumed that the distinction between the two centers 
of gravity is of no account in view of the light mass of the electrons. However, we shall see later in 
the article that the S,,, term is often as important as the P,,, Q,, ones, and that all three should be 
considered as an organic whole in treating the coupling of electronic and nuclear motions. The 
formula for Sy’, is 


Snin(r) = — (h?/82? M142) Sf: : fo." (SV 2423 ;5,V ;-V)’,d7, (14) 


where M42 is the total mass M,+ M, of the two nuclei. The symbols V; and V,? are respectively the 
gradient and Laplacian operators in the x;, y;, 2; space and we shall let Vxyz, V’xyz have similar 
significance in the X, Y, Z coordinate system. 

To derive (14), we note that motion of the center of gravity introduces the term — (h?/87?Mi42)V?xvz 
in the Hamiltonian function. Despite the fact that & does not apparently contain X, Y, Z explicitly, 
Vxyz® does not vanish, as in differentiation with respect to X, Y, Z it must be supposed that the 
electronic coordinates are kept fixed in space rather than fixed relative to the nucleus. Hence when 
X is given a variation 6X, one has 6X,;=--:=6X,=—6X and so 


—(h?/8n2My42)V 2xy 28 = — (h2/822M 112) Rul ®,V 2xyzf—WxvzfOV Oat f(SV 2b, +22)5,V 0). (15) 


To obtain (14), one integrates over the electronic coordinates in the usual way. The cross terms 
in Vf-V®, in (15) vanish on integration in the diagonal case n’=n because the mean momentum 
of the electrons relative to the center of gravity is zero. The first term of (15) merely governs the 
uniform motion of the center of gravity, and need not be considered for our purposes. 

To proceed further, we develop Pan+Qnn+Snn as a Taylor’s series 


Pant+}QnatSnn=Bolkothiét hot? +--+] | (16) 


about the equilibrium position ro associated with (4). A perturbation calculation then shows that 
the resulting correction factors entering in (5) (6) are 


Yo = 4k Bo? /h?w?, = (2ke—3k,a 1) Bo’, h?w?. (1 7 ) 


A simple proof of (17) is as follows: After the addition of the Pan, Qun, Snn and centrifugal terms, the equilibrium 
distance is no longer quite ro, but is very approximately 






































he i ee Seca: ha a 5 tS Bm 





























H. VAN VLECK 





ro’ =ro— (ki Bo/aro) +2(Bo/aro) J (J +1) (18) 


as is seen by taking the minimum of the total effective potential in (11) and retaining only terms of lowest order. When 
we expand about the new equilibrium position, this potential becomes 


Wot Bol J(J+1) +ho]+ 4(Bo?/aro?)[— ki? +4ki J (J +1) —4J2(J+1)2]4+ fa(r—ro')? 
+ Boro *Lk2+3J(J+1) — 3a:(3ki —6J?—6J) (ro—ro’)?+°++. (19) 


In accordance with perturbation theory, we may replace (r—7o’)? by its mean value, which to a sufficient approximation® 
is the same as that of (r—7o)? for the unperturbed motion. The latter is immediately obtained from the theorem that the 
mean kinetic and potential energies for a harmonic oscillator are each half the total energy, i.e., 


1a(r—ro)? = 3p,2/M = }hwo(v+}). (20) 


Thus W = Wot BoLko— (ki Bo/hwo)? ]+ Bo(1+-ys)J(J+1) +hwo(1 +6) (v+ 3) + DI2(J+1)2?-—a(v+3)J(J+1), (21) 
where 7s, & are defined as in (17); and where 
D = —4By3/h2wo?, a = —6(Bo?/hwo)(1+a;). (22) 


Eqs. (17) are thus established. The formulas for D and a are the usual ones of the band spectroscopists for centrifugal 
expansion and the interaction of rotation and vibration.’ As an incidental to our work, it is thus seen that the a, as well 
as, of course, the D term can be derived by a quite elementary perturbation calculation rather than, for instance, em- 
ploying the Kratzer potential.” If we are interested in the effect of the various extra terms on the vibrational wave func- 
tions, this is, to a sufficient approximation, merely a shift in origin by an amount ro’—7ro given by (18).8 

Eq. (21) does not contain higher powers than the first in v+}, as we have in (19) neglected terns of the third and 
higher orders in (r—r9’) (also unimportant terms proportional to (ro’—7)?(r—7o) or (ro— 70’), etc., have been discarded 
in Eq. (19)). For most purposes, the coefficients of the higher powers of +} can be considered the same as when the 
extra terms in (11) are omitted. The coefficient of v+} is slightly modified by the anharmonicity, essentially because 
v™ in the old quantum theory is replaced in the new mechanics by a polynomial in v+} rather than simply (v+34)™. This 
is the Dunham-Crawford effect already included in (a). 

In writing (19) and (21), we have not carried the development far enough to find the small correction factors to the 
coefficient of (v+4)J(J+1) which are comparable with those in (5)—(6) and which depend on the isotope. Further cal- 
culations, in which the equilibrium distance must be obtained more accurately than in (18), show that a more accurate 
value of this coefficient than a@ as given in (20) is? , 


a* =at (aBo?/h2wo?) [Saiki + (1 +0;)-1(8k1 — 2k2— 6a;k2+4hk3—8a2k,) |. 
The corresponding correction formula to the coefficient —xhwo of (v+ }) is 
(xhw)* = xhwo+ (4Bo%x/hwo) (12a2— 15a,2)—1( — 30a3ki + 12k4— 12a2k2+78ai:42k; — 30a1k3+30012k2 — 45a;°k;). 
In the present case n’ =n, the expression Q,,, defined in (13) may be transformed into 
Qnn(r) =(h?/82?M)Onn with Onn=Sf--: f'|d6n/dr|%dr, (23) 


and is therefore positive. The relation (23) is proved from (13) by differentiating twice with respect to r the normalization 


integral 
S + SON Ondr =5,.™. (24) 


Eq. (24) is valid for all values of the parameter r since throughout the article we shall suppose the electronic wave 
functions normalized to unity regardless of the internuclear distance. The result (23).then follows immediately in the 
real case ,'=,. If &, is complex, the proof can still be carried through when use is made of the fact that without 











5 This statement does not contradict the fact that it is 
precisely the difference between ro and ro’ which causes 
the effects in question, as the shift in origin has eliminated 
first order expressions in r—7r9', and the second-order terms 
need not be evaluated as accurately. In calculating a*—a 
it is necessary, however, to distinguish between the means 
of (r—ro)? and of (r—ro’)?. 

6 Cf., for instance, R. S. Mulliken, Rev. Mod. Phys. 1, 
67 (1930). In general we use the notation approved by 
Mulliken. His anharmonicity coefficients, however, relate 
to expansion about a Kratzer rather than harmonic poten- 
tial, and so differ from ours. Thus his —2b* —2 is the same 
as our 4. 

7 The linear term in £ in (16) or in the expansion of the 
centrifugal potential no longer vanishes when the potential 
is unsymmetrical. This is the physical reason why the 
cubic term a; in (3) contributes to the quadratic part of 





(19). Neglect of this contribution explains why Morse 


originally obtained an incorrect formula for @ (Phys. Rev. 
34, 57 (1929)). 

8 Unless one is interested in extreme accuracy it is thus 
not necessary to resort to an elaborate perturbation cal- 
culation such as that of Pekeris (Phys. Rev. 45, 98 (1934)). 

9 There are also modifications in @ due to effects (a), (); 
(d) and in xw due to (a) and (d). Dunham? has treated (a). 
It can be shown that with the same hypotheses as in (28), 
(c) alters a by a factor 

1+ Bo(5+6a;)2(L?+L)/3(1+a,)hv(2; M1). 
The correction formulas due to (d) involve Bo/hy, rather 
than Bg?/h®wo? like (b). We shall not give their explicit 
form, since they are rather complicated and hard to use 


quantitatively as they require knowledge of the constants 
nan’, Lan defined before (37). 
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loss of generality, ®, may be assumed to have the structure F+iG, where F, G are real and respectively even and odd 
with respect to reflection in some given plane containing the molecular axis. Consequently, F, 0F/dr, 0?F/dr? are al! 
orthogonal to G, or derivatives thereof, and no cross terms in FG are involved in (24) or (13). 

Furthermore, (14) may be transformed by partial integration into 


Sna(r) = (h?/82?Mi42) S- ates S|%i 


Vib, |%dr (25) 


and is also positive. There is no trouble in the complex case, as V;F is orthogonal to V;G, and F to V,°G. 


To simplify (17) it is customary to assume 
that Qnn, Sn are small compared with P,,, and 
to make use of the hypothesis of pure precession, 
in which the square of the total angular mo- 
mentum is assumed to have a quantized value 
L(L+1) like in an atom. Since in reality only the 
axial component A of electronic angular mo- 
mentum is constant, this approximation is a 
crude one strictly allowable only in the limiting 
case of the ‘‘united atom”’ formed by coalescence 
of the nuclei. With the hypothesis of pure pre- 
cession, the integral =,, in (12) has a value 
L(L+1) —A? independent of 7, and then ,, ke 
n (16) owe their existence (provided we dis- 
regard Qnn, Snn) entirely to the factor 1/r? in 
(12). If we use (22), thus Eqs. (17) reduce to 
formulas 


yo= —8LL(L+1) — A? ]Bo?/h?w?, 


b= —(L(L+1) —A?Jo/he (26) 


previously given by Kronig" and Dieke," re- 
spectively. Here a is defined as in (22). It is to 
be cautioned that use of (26), even in a first 
approximation, is allowable only when the 
molecule is in a paramagnetic state L~0 in the 
united atom approximation. When L=0, it is 
not allowable to disregard Qnn, San in com- 
parison with P,,, or overlook the variation of 
the integral in (12) with r. Even when L #0, this 
procedure is questionable unless we are dealing 
with highly excited states in which the molecule 
behaves nearly like the united atom. 


(c) L UNCOUPLING 


The remaining isotope corrections (c—d) arise 
because the separation of variables (9) is not 
accurate, or in other words because rotation and 
vibration introduce interactions between dif- 
ferent electronic states. These interactions can 
be divided into two types: (c) Those connecting 


”R. de L. Kronig, Physica 1, 617 (1934). 
"G. H. Dieke, Phys. Rev. 47, 661 (1935). 





states whose A’s differ by one unit and (d) those 
joining states of equal A. 

We shall not write down the perturbing terms 
in the Hamiltonian function which introduce 
matrix elements between states differing by 
unity in A. Such terms are well known, and give 
rise to the phenomenon commonly called 
L uncoupling. They make B* different from B, 
and so influence the isotope ratios, as Professor 
Mulliken first pointed out to me.'? In the case of 
a > state, calculation shows! that 


¥e=4BodLiy| P(N; = */hv(z ; I), 
P,(llz)=f:- fon P,®sdr, 


with P, defined as after (13) and where the sum- 
mation is over all II levels. In writing (27), we 
have assumed that in electron bands, the elec- 
tronic frequencies are not appreciably modified 
by vibration. With the hypothesis of pure pre- 
cession explained above near the end of fine 
print, only one of the II states gives a perturba- 
tion, and (27) reduces to 


Ye=2BoL(L+1)/hv(= ; M1). 


The reader is referred elsewhere!‘ for the for- 
mulas in the case of II levels. Here B* is different 
for the two components of a A doublet. In fact, it 
is precisely this difference which causes the 
doubling. 


(27) 


where 


(28) 


12 It is hard to say where it was first seriously suggested 
that Z uncoupling modifies the isotope effect. Even before 
wave mechanics, Mulliken intimated that interaction 
between electronic motion and vibration might modify 
isotope behavior (Phys. Rev. 25, 126 (1925)). The requisite 
formulas of quantum mechanics were given in reference 4, 
while later papers by Mulliken repeatedly stressed the 
distinction between B* and B, including the case of 2 states 
(cf. Mulliken and Christy, Phys. Rev. 38, 92 (1931), Eq. 
(8)). Very specific applications to the isotope effect have 
been made by Dieke.?° 

13 Mulliken and Christy’s Eq. (8) has a factor 8 rather 
than 4 as in our (27) because their convention is to include 
in the II summation only one A-doublet component in a 
A=1, —1 system of representation. In such a system both 
components contribute equally, whereas only one com- 
ponent is effective if passage is made to wave functions 
of the proper final + —symmetry 

14 J. H. Van Vleck, Phys. Rev. "33, 487 (1929). 
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L uncoupling is primarily a rotational effect, 
and so one would expect it not to affect the vibra- 
tional isotope ratio. More detailed investigation 
shows this to indeed be the case, so that one can 
take 

5.=0. (29) 


(d) INTERACTION BETWEEN STATES OF EQUAL A 


We have postponed discussion of this inter- 
action until now, as its influence on w* has not 
previously been treated, and requires some 
study. Let ” and n’ be two different electronic 
states of equal A whose vibrational quantum 
numbers are v and v’. The perturbing Hamil- 
tonian matrix element connecting them is'® 


Hyn’ v30)=SR nel Pane+Qnn' t+Snn’ 
+(—th/2r7M)Qnw pr IRnwdr, (30) 


where p, is the operator hd---/2mr10r, where 
Pan's Qnn’s Snn’ are defined as in (12), (13), (14) 
and where 

Qan = SP, (db, /Or)dr. (31) 


The centrifugal term [/J(J+1)—A?*]h?/82?°Mr? 
does not appear in (30), as it makes no con- 
tribution because of the orthogonality of ®,, ,:. 
Hence the present effect (d) does not influence 
the rotational isotope behavior, and we can take 


ya=0. (32) 
On the other hand, it gives a modification in the 
vibrational case, as we shall show that 


ro” 0°,,n°|2Bo 
ba=22 0 ; (33) 
hv(n;n’) 





Here the summation is over all states of the 
same A as the given state n. When dealing with 
> states, however, only states of the same + — 
symmetry need be considered, as (30) vanishes 
if m relates to a +, and n’ to a =~ state, or vice 
versa.!6 The zero superscript in (33) means that 
(31) is to be evaluated at r=7p. 


15 See reference 4 for the mathematical basis underlying 
the construction of the matrix elements (30). Besides the 
terms in (30), there are nondiagonal matrix elements due 
to the cross products between Vf and V®, in (15), which 
no longer necessarily give a vanishing effect when n#¥n’. 
The resulting contribution, however, does not depend 
appreciably on the vibrational quantum number, but 
rather primarily on the translational energy of the mole- 
cule, and so can be omitted for our purposes. 

16 For what is meant by =* and =~ states, see R. S. 
Mulliken, Rev. Mod. Phys. 3, 92 (1931). 
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It will be noted that (33) involves Q but not 
P,Q, S despite the fact that P, Q, Sand Q appear 
additively in (30). The reason is that the 
matrix elements due to P, Q, S are smaller than 
nonvanishing elements due to (h/27M)Q by a 
factor of the order of magnitude (Bo/hwo)! as is 
seen by a simple dimensional examination. Hence 
we have not bothered to include any of the 
former in (33), nor shall we in future Eqs. (36), 
etc. The physical basis of this diversity in size 
is that R,-. varies much more rapidly with r 
than #,- does; i.e., the operator hd---/2ridr 
applied to ®,, as in (13) is equivalent in order of 
magnitude to multiplication by 4/2mro, whereas 
when applied to R,.,, it is the vibrational mo- 
mentum operator p,, which is of the order 
(Mhwo)) by (20). The larger size of the © than 
the other perturbing terms appears to have 
previously been overlooked, and explains why 
the Q terms are able appreciably to modify w* 
despite the fact that they appear only connecting 
different electronic states, and so enter finally 
only in formulas involving the electronic fre- 
quencies in the denominator. The usual Z un- 
coupling terms, to be sure, are able to modify B*, 
but then B* represents a smaller energy term 
than w*. Fortunately, the 2 terms vanish in the 
diagonal case =n’, as one can prove that ),, 
equals zero by differentiating (24) with respect 
to r. The proof applies even when ®,, is complex. 
provided one utilizes the structure of ®, explained 
after Eq. (24). Since 2,,=0, we omitted any 
mention of the 2 terms in discussing the diagonal 
effect (b). 

Proof of (33). By a common formula of perturbation 
theory, the shift in energy due to the totality of matrix 
elements (30) is 

AW = Swen (v3 vA an(v' 3 v)/hv(nv; n'v'). (34 
Usually the denominators of (34) depend but to a small 
extent on the vibrational quantum number, and may be 
replaced to a sufficient approximation by expressions 
hv(n; n’) which can be taken outside the v’ summation. 
Then (34) becomes 

AW =2n{[EeHan(v; 0’ )Hnn(v’ 30) Vhv(n;n’)}. (35 
Furthermore 
LwHan(v; v')Ann(v' ; v) 

=(—hi/2eM)?S Rut Quy PrQninPrRadr. (30 
It is clearly to be understood that the summation in (30 
does not lead to true matrix multiplication, as the initial 
and final wave functions in (33) belong to different sets 0! 
orthogonal functions because the states m and n’ have 
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different radial constants. Summation over »’ in addition 
to v’ would be necessary to give the ordinary case. How- 
ever, the writer has shown elsewhere that even without 
summing over 7’ it is allowable to use formulas such as 
36)." 

The important point is that in (36) we have an expression 
which does not require knowledge of the vibrational wave 
functions for the perturbing electronic state n’. To evaluate 

36) we expand Q,,,° (and Q,-,) as a Taylor’s series 


Quan = nn + Q n(r—7o) + eye: 


about the equilibrium position for the state . For our 
purposes, only the constant term {2°,,,,, is needed, as the 
remainder would introduce higher powers of (v+}) or of 
B/hwo. With this approximation, (35) becomes 


AW = (h?/4x?M?) p?,.,. (v3 v)Zn | Qnne|2/hv(n; nn’). (37) 


Here use has been made of the relation Q,-,=—Q na.’ 
obtained by differentiating (24) with respect to r. The 
diagonal element p?,,,(v; v7) is the same as the mean value 
of p,? and so can immediately be obtained from (20). The 
latter shows that (37) is proportional to v+} and so con- 
tributes to the w* part of (1). Formula (33) for 64 then 
follows when one uses the definatory relation (6). 


SIGNS AND RELATIVE MAGNITUDES OF THE 
VARIOUS CORRECTIONS 


First we shall consider the subject of sign. 
Usually the empirical isotope ratios are nearer 


unity than are the ideal values (4). This behavior 
is predicted by theory if y and 6 are negative 
since in all four cases the corrections to (4) 
become more important the smaller the reduced 
mass.'* Now any of the four effects (a—b—c-d) are 
more apt to give negative than positive contribu- 
tions to y or 6 and so the observed behavior is 
quite understandable. The detailed situation on 
sign for each effect is as follows. In (a), the sign of 
Y. or 6, depends by (7-8) on the a; in (3) or in 
other words on the type of anharmonicity or 
exact shape of the potential curve. In only two 
cases, viz.,!° the (2sc)?'Z and 2se2pa'Z states 
of LiH, is this information known with sufficient 
accuracy to permit the calculation of ya, 6. In 
both cases yq is negative, also 6, for 2sc2pe '%, 
but 6, is positive for (2sc)? 4Z. The values of yo 
and 6, are negative if one uses the approximation 


J. H. Van Vleck, Proc. Nat. Acad. Sci. 15, 754 (1929). 

'* In all cases (a)—(d), y and 6 involve factors Bg?/h?w¢? 
or Bo/hv and so are apparently inversely proportional to 
the reduced mass M. However, the terms kk» in (17), or 
=, © in (38) may involve M,/M:. This is shown by (43) 
or (46-7-8). Even though the mass dependence can thus 
be more complicated, the dominant factor is usually 1/M. 


,_’ F. H. Crawford and T. Jorgensen, Jr., Phys. Rev. 49, 
745 (1936). 
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(26) based on pure precession, since the spectro- 
scopic constant @ is usually positive. The same 
sign behavior will also be found true (cf. Eqs. 
(43)) when we develop approximations for other 
limiting conditions. So y, and 6 are probably 
usually negative. Both 6, and ya vanish, while 
Eqs. (27) and (33) show that y, and 64 are 
negative if the dominant frequency denominators 
are negative (i.e., relate to absorption), which is 
surely the case if the given state is the ground 
level of the molecule. 

As to the relative importance of the various 
corrections, it is hard to make a definite predic- 
tion. The effect (c) should have a preponderant 
influence if there is close by a state differing by 
one unit in A. Since 6,=0, the departures from 
(4) should then be very much greater in the 
rotational than in the vibrational isotope effect. 
Dieke* has shown that an interesting example of 
precisely this situation is provided by the 3p *= 
state of the hydrogen molecule, and that the 
observed deviation from (4) agrees remarkably 
well with that calculated from (28), which is 
here more important than any other correction 
by a factor ten or so. 

Similarly, (d) should be the most important 
effect if there is a nearby state of the same A and 
same +— symmetry as the level under con- 
sideration. Since yz=0, there would then be a 
much larger correction to the vibrational than 
to the rotational isotope effect. It would be very 
gratifying if this point could be tested experi- 
mentally. 

When there are no nearby states, so that one 
does not have the resonance cases just cited and 
instead the frequency denominators are com- 
parable with ionization potentials, it is impos- 
sibie to say which of the various corrections is 
the largest without detailed and laborious inves- 
tigation. The different effects (a—b—c—d) are then 
of the same rough order of magnitude, as we will 
later show more fully by some illustrative cal- 
culations (except for effect a) on the ground 
state of the hydrogen molecule. Too much weight 
should not be attached to the fact that (7), (8), 
and (17) are proportional to By?/h*w»?, and (27) 
and (33) to By/hv, where wy and » are respectively 
vibrational and electronic frequencies. To be 


sure, Bo?/h®wo?, being quadratic in Bo, is some- 


20(G. H. Dieke, Phys. Rev. 48, 606, 610 (1935). 
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what smaller than Byo/hv even when one is well 


away from resonance. For example, the former 
is 0.0002 and the latter 0.0006 in the case of the 
ground state of the hydrogen molecule, a rather 
extreme example of nonresonance as v=0.9 x 105 
cm~!. However, the difference in magnitude is 
not large, and we shall later see (cf. Eqs. (43)) 
that under certain limiting conditions (17) can 
be transformed into formulas not involving wo 
at all. Any disparity resulting from the diversity 
in size of Bo?/h?wo? and Bo/hy can easily be offset 
by abnormally small or large numerical coef- 
ficients in the formulas. In particular the large 
numbers occurring in (7) and (8) are to be noted 
as they may make the anharmonic effect (a) 
fully as important as any of the other corrections 
except near resonance. However, Crawford and 
Jorgensen” find that in LiH, (a) is able toaccount 
for only a fraction of the empirical isotope cor- 
rection ; namely, the empirical values of y, 6 for 
the (2sc)* 14 state are — 0.0005, —0.00019, while 
the calculated values of ya, 5, are —0.0001, 
+0.00008; for the 2sc2po'X level, y, 65 are 
—0.010, —0.027, and ya, 5. give —0.003, —0.004. 
It is interesting that y./y has about the same 
value for both levels despite the fact that one 
has a y twenty times as large as the other. The 
large value —0.027 of 6 in the upper state is 
perhaps indicative of a large effect of type (d). 
At any rate, the inadequacy of (a) alone seems 
to corroborate our statement that (a—b-—c—d) 
likely simultaneously can play an appreciable 
role. When the united atom has an angular 
momentum, (b) is more important than other- 
wise (cf. remarks following Eq. (40)). 


DIGRESSION ON PERTURBATIONS IN 
BAND SPECTRA 


At this point we shall insert a few words on 
perturbations in general in diatomic spectra, for 
the underlying mathematical analysis is the 
same as for the isotope corrections. Dieke*! has 
called attention to the fact that perturbations 
connecting states differing by one unit in A 
increase rapidly with the rotational quantum 
number J, all other things being equal, and cites 
confirmatory experimental evidence. Eq. (37) 


(4. H. Dieke, Phys. Rev. 47, 870 (1935). 


shows that interaction between states of equal 
A increases with the vibrational quantum number 
v. This fact has not apparently been stressed pre- 
viously, and it would be very interesting if appro- 
priate experimental evidence could be obtained .?'* 
The confirmation would be more difficult than in 
the preceding case AA= +1, for the dependence 
on the quantum number is masked if there is 
much alteration in the frequency denominator, 
and of course the frequency depends much more 
on the vibrational than on the _ rotational 
quantum number. In other words, in studying 
perturbations in band spectra, one is apt to be 
near resonance, and this is just the place where 
one cannot pass from (34) to (35). Effect (c), 
L uncoupling, is primarily a rotational phe- 
nomenon, whereas (d) is due mainly to vibration, 
and this distinction explains why one type of 
perturbation should vary markedly with J and 
the other with v. 


SIMPLIFICATION OF THE FORMULAS UNDER 
LIMITING CONDITIONS 


For estimating orders of magnitude, it is con- 
venient to simplify some of the formulas further. 

Centroid approximation. If we assume that all 
the frequency denominators can be regarded as 
having sensibly the same value hy, an assumption 
warranted if at all only when one is well away 
from resonance, then (27) and (33) become 


c= 2Bozo/hy, Oa = 2Bory?@o/hv, (38 


where Zo and ©» denote the integrals defined in 
(12) and (23), respectively, evaluated for r=70, 
n'=n. The first formula of (38) follows from 
(27) by simple matrix multiplication, if v is 
constant, inasmuch as P,? and P,? contribute 
equally to =. The second relation is established 
from (33) in an analogous fashion by noting that 
the Q,,, are the coefficients in the expansion of 
d®,,/dr in terms of the orthogonal functions ?, 
(cf. Eq. (31)). 

Behavior of Pun, Qnny Snn for separated and 
united atoms. Especially in connection with the 
specific isotope shift to be considered later, it is 
illuminating to study the limiting values of P, 


218 Note added in proof. Dieke finds that an example is 
furnished by the 1so2sa' state of Hz (paper presented 
at Washington meeting of Am. Phys. Soc., 1936). 
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ISOTOPE CORRECTIONS 
Q, S, defined in (12)—(13)—(14) when the nuclei 
1, 2 are infinitely far apart or coalesced together. 
The results are as follows 


Pan OF 2Qnn—>3(Miwe+ Mow,)/ M142, | 
= 2) « 9 
F = Su n—> (Mw, 4. Mow.)/M, $25 | (3 ) 
0 (P,, or 20 nn—>§ (MiZ2— M2Z;)*w142/M; M2Z%142,\ (40) 
a n n—>Wi+2- J 


Here wy), We denote the energy corrections for the 
motion of the nucleus for isolated atoms 1, 2, 
respectively, and w1,2 that for the united atom. 
Z; and Z, are the nuclear charges of atoms 1 and 
2. Eq. (39) applies only if the dissociation 
products are atoms in S states, as the proof 
assumes that the kinetic energy, etc., is prorated 
equally between the x, y, and z directions. If this 
condition is not fulfilled, the sum Pant+Qnn is 
still given correctly by (39). Similarly, if the 
united atom has an angular momentum, (40) 
should be used only to obtain the sum Pyn+Qnn, 
and furthermore in this case there is the extra 
term [L?+L—A?]/8x?Mr? not included in (40) 
(cf. Eq. (26)). Its inclusion means that P,, 
becomes infinite when r—-0 if L #0. This appar- 
ently absurd result is obtained because the per- 
turbation calculations all assume that the elec- 
tronic charge cloud rotates with the nuclei; 
actually the two will become decoupled when r 
is made too small. However, for actual values of 
r, the decoupling has usually only started, and 
the large values in the limit r=0 explain why 
P,..(r), and hence effects (b)—(c), are abnormally 
important when the united atom has an 
angular momentum. It is usually assumed that 
in (12) the integral =,,, does not vary much with 
r, but (39) shows that for large values of 7, the 
factor Z,, is proportional to r? inasmuch as 
Pan=EZnnh?/8x2Mr* is independent of r when 


—[—> oO. 


A sketch of the arguments used to obtain (39)—(40) is as 
follows : The values of S,,, for r= «& and r=0 are quite ob- 
vious, as for r=0 S,,, is simply the usual correction term for 
nuclear motion,** while for r= © it would be except that it 
has a mass factor 1/M,42 rather than 1/M, or 1/M>. The 
behavior of Pan+Qnn at r= is obtained by noting that 
when r is varied with the center of gravity held fast, the 
origin of atom 1 is displaced by an amount M2/ Mj,» and 
that of atom 2 by M,/M,42. Furthermore the P,., and Qna 
terms then combine to involve essentially the Laplacian 
operator V2xyyz, as far as each nucleus is concerned, except 
for a numerical factor —(h2/8x*)(M2?/MM%*,,2) or 





2 Cf. E. A. Hylleraas, Zeits. f. Physik 54, 364 (1929). 
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— (h?/82?)(M,?/MM?,,.). Here M is the reduced mass 
M,M2/M,42. As already stated, the apportionment be- 
tween Prn and Qna given in (39) is correct only if the 
isolated atoms are isotropic. The behavior of Pan+Qnn at 
r=0 follows on observing that if the nuclei are drawn 
infinitesimally far apart with their center of gravity pre- 
served at the origin, the change in the potential field may 
be regarded as equivalent to that of a dipole, and the 
dipole is (Z;:M2—Z2M,)/Mi42Z142 times as intense as for 
the purely translational or S effect. This accounts for a 
factor (Z;M2—Z2M;)?/M%449Z7142 in the behavior of 
Prn+Qnn as compared with S,,; the remaining factor 
M?,,2/M;Ms¢ is because of the different mass denominators 
in (12-13) and (14). It may seem strange that P,, does 
not vanish at r=0 when L=0; the reason for this paradox 
is that the center of gravity does not coincide with the 
nuclear midpoint, so that the angular momentum about 
the former need not vanish even though that about the 
latter does. 


If one neglects the cross terms j7#7 in (14), 
the expression w; equals m/M, times the kinetic 
energy of atom 1. With a similar approximation, 
Wi42 is m/My,2 times the kinetic energy of the 
united atom, where m is the mass of one electron. 
Actually, nonvanishing cross terms owe their 
existence to the fact that the complete wave 
function is not expressible as the product of 
wave functions of individual electrons, either 
because of antisymmetrization, or because the 
interelectronic distances enter explicitly as 
arguments. The calculations of Hylleraas®? show 
that the cross terms resulting from the latter 
cause, are unimportant, at least in the ground 
state of helium. The work of Hughes and Eckart*® 
indicates that in atoms the antisymmetrization 
affects the nuclear corrections appreciably only 
if at least one electron is in a p state. In mole- 
cules, it will give an appreciable contribution at 
medium internuclear distances even when elec- 
trons are only in s states, due to the overlap of 
wave functions of the different atoms. However, 
at any internuclear distance, we can neglect the 
cross terms if we are interested only in orders of 
magnitudes, and so a very simple rough fermula 
for Sn» is 


Snn=(m/Mi42)T, (41) 


where TJ is the total kinetic energy of the elec- 
trons. The virial theorem” tells us that 


T= —E,—rdE,/or, (42) 
23D. S. Hughes and C. Eckart, Phys. Rev. 36, 694 (1930). 


* Cf. J. C. Slater, J. Chem. Phys. 1, 687 (1933); W. 
Pauli, Ann. d. Physik 68, 177 (1922). 
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and so if (41) is accepted, the variation of S with 
r, as required for (16), can immediately deter- 
mine if we know the potential energy curve 
E,(r), as in (3). In (42) the origin for the energy 
E,, must be so chosen that £,, vanishes when the 
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in symmetrical molecules the P and Q terms are 
negligible compared with S at small internuclear 
distances, provided L=0, and so for estimating 
orders of magnitude, it is natural to neglect P, Q 
in comparison with S, and use (41-42) and (3). 


Then (17) becomes 
Yo = —2m/M 442. 


molecule is blown completely to pieces (in- 


dividual isolated charges). Eqs. (40) show that 65= —3m/2Mis2. (43) 





CALCULATION OF Pun, Qnn, Snn FOR WANG HYDROGEN FUNCTIONS 


The values of Pan, Qnn, Snn can be computed for any value of 7 for the normal state of the hydrogen 


molecule if one use the well-known Wang” wave functions 
@,, = A (e~ ("1a+720) 4 e- b(ribt+r2a)) | (44) 


Here 71, is the distance of electron a from nucleus 1, etc., and the normalization constant A has the 
value 
1/A?=(27?/b*)[1+e-2°"(4b*r? + br +1)? ]. 


(45) 
By simple integrations, most easily made in elliptic coordinates, one finds that 
Pn=Bo(42?r92A?/18b4) [3 —e-?°"(b?r? + 3br +3) (br +1) J+ [(Mi— Me) /2M 42 |? Bog, 
Qnn= —Bo(rod log A /dr)?+[2(Me?+ My?) /M? 142 laBo—4(Mi1Me2/M?142)BBo, 
Snn=Bog(Mi M2,’ M*142) +4(M1Me2/M?*142)(a+ 8) Bo, 
g = (8279°A?/3b4)[ 1 +e-?°"(br+ 1) (3b? ++. br+1) ], 
a= (2?79?A?/3b*)(1+3(br+b'r?)*e-?*” 
B = (97b?792A?/3b4) (—b?r? +-br+1)(1+br+ 1b?r?)e-25", 


(46) 
(47) 
(48) 
(49) 
(50) 
(51) 


These formulas conform to (39-40) in the limiting cases.*® 


ESTIMATE OF THE VARIOUS TERMS FOR THE GROUND STATE OF THE HYDROGEN MOLECULE 


The shape of the potential curve is not known with sufficient precision for the ground level of 
hydrogen so that ya, 6. can be evaluated. However, the preceding calculations enable us to determine 
roughly the corrections due to the other causes (b—c—d). We shall let primes refer to He, double 
primes to HD, and triple primes to De. When Wang’s optimum value 1.166/a» is substituted for 
the screening constant 6 at the actual internuclear distance r=1.406a(a)=h?/477e?m), (46) and 
(47) give 


0 =20 ’=0.394, =’ =0.486, Q,'= Oo’ =0.335/ro?, Oo’ =0.377/ro?. (52) 


For the centroid frequency v, following Frank and V. V.,?" we shall use the value 1.23R, where R 
is the Rydberg constant. Then (38) yields 
v-' = —0.00035, y,.’’= —0.00032, 
6a’ = —0.00030, 6,’’= —0.00025, 


y."’ = —0.00017, 
64’ = —0.00015, 


while from (43) we have 
yo = —0.00054, y,’’ = —0.00036, 


6,’ = —0.00041, 6,’’ = —0.00027, 


yo’ = —0.00027, 
6,’ = — 0.00020. 


25S. C. Wang, Phys. Rev. 31, 579 (1928). 
26 The values of w; and wi,2 yielded by the Wang func- 
tions are, respectively, h?b?/82?M; and 2b?h?/87?M, 40. 


27 A. Frank and J. H. Van Vleck, Proc. Nat. Acad. Sci. 
15, 539 (1929). This paper gives details of the calculation 
of Zo for He. 
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Eqs. (54) should be regarded as probably giving only upper limits to the absolute values of +, and 6, as (43) assumes 
that Pan+Qnn can be neglected in comparison with S,». Since the w’s are positive, Eqs. (39-40) show that Pran+Qnx is 
larger for r= © than for r=0 whereas the reverse is true of S,,. Hence 0(Pnn+Qnn)/dr is doubtless opposite in sign to 
OSnn/Or. SO Pan+Qnn gives a contribution to &; in an opposite direction to that from S,,, (cf. Eq. (16)) and so by (17) 
tends to offset the S,,, effect in y, and in at least part of 5. Of course one can reckon yz, 6, without disregarding Pnn+Qnn 
in comparison with S,,, by expanding (46)—(47)—(48) about the equilibrium position. If one neglects the variation of 
b with 7, then one obiains a value 0.0004 for +» which is of the opposite sign from that given in (54). In all likelihood yp, 
is actually negative since (39)—(40) show that Prn+QnatSn. is greater for r=0 than for r= ~,% making 0(Pant+Qnn 
+5S,,)/8r presumably negative. On the other hand (46)—(47)—(48) yield a smaller value of Prn+Qnant+Snn. for r=0 
than for r= © if the same value of b is used in both limits. This fact makes it clear why expansion of (46)—(47)—(48) 
gives a yp of the wrong sign if } is treated as constant. A more refined calculation, in which } varies with r, does not 
seem repaying at the present time. Even the leading term koBo in (16) as obtained from (46)—(47)- (48) involves some 
error because b was regarded as independent of r in computing the integral ©,,, in (23). However, the error is probably 
not as serious as in the higher 7 derivatives involved in connection with k;, ke in (16), and Wang’s success in computing 
wy without letting 6 depend on ,r is indicative that the variation of b with 7 is not too important in lower derivatives. 


Unfortunately, no adequate experimental data on the ground state of the hydrogen molecule are 
available to test (53)—(54). For excited states, Dieke!! calculates from (26) isotope corrections not 
a great deal larger than (53)—(54), being still of the order of magnitude y’—~y’’~10~‘. The effect 
(26) considered by him is present only when the united atom has a nonvanishing angular momentum 
L. The ground state of the hydrogen molecule has Z = 0, and our main purpose in presenting (53)—(54) 
is to show that even without the anharmonic effect (a), not all of the isotope correction is connected 
with the existence of an L +0. 


THE ELECTRONIC ISOTOPE SHIFT 


The difference in the constant terms in the energy for two isotopes is usually called the electronic 
isotope shift. It consists mainly in the difference in the constant part koBo of (16), or in other words 
in Pan t+QnuntSnn evaluated at r=70, for the two isotopes. Dieke has had some success in explaining 


approximately the observed isotope shifts by neglecting Q,,, S,, in comparison with P,,, and 
assuming for the latter the value 


Pun=Bo[L(L+1) —A?] (55) 


furnished by the hypothesis of pure precession. A somewhat better approximation is probably to 
add to (55) the value of S,,, furnished by (41)—(42). As noted by Dunham? and Dieke,"! the anhar- 
monic correction (a) also makes a contribution 


Yoo = (Bo/8) (3a2—7a,/4) (56) 


to the constant part of the energy. In the cases studied by Dieke,'! the effect of (56) is unimportant. 
The inclusion of .S,,,, however, alters the calculated values somewhat, as shown in the following table. 
Our various formulas, though ostensibly for singlets, can be applied to triplets in hydrogen since 
here the multiplet width is negligible. 


He—-HD HD—Dz2 
CALCULATED OBs. CALCULATED Oss. 


(55)—(56) (41)—(55)-(56) (55)—(56) (41)-(55)—(56) 
3p 711 —2s 35 7.60 6.07 cm™ 5.01 ; 7.66 6.89 cm- 4.94 
3p *X —2s 3y 13.00 11.92 13.23 12.17 11.64 12.28 

Clearly it is not allowable to neglect Q,,+.5S,, in comparison with P,,, unless at least either the 
initial or final state has L=0. Even when this condition is met, it is no longer allowable to overlook 
Qn tSnn if the line is in the ultraviolet. By (41), the S,, effect is more important the shorter the 
wave-length, as then there is more difference, and hence less cancelation, between the values of (41) 
for the initial and final states. Examples are the lines, all in the ultraviolet, terminating on the ground 
state of the hydrogen molecule. To make an approximate calculation, we can use (46)—(47)—(48) 


*S If one neglects the cross terms in (14), then use of observed energies and the virial theorem shows that the 
actual values of w, and w.,2 for H and *He: are, respectively, R/1840 and 5.8R/3680, where R is the Rydberg constant. 
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for the ground state, and assume that in the upper state, the effect of the inner electron can be com- 
puted by using the Finkelstein-Horowitz”® wave function C(e~*"*+e-*™) for H2t+. The latter yields 


Pan t+QnantSnn=Bol(C?aro?/d) — (rod log C/dr)? ] (57) 
with d*/27C?=1+e-4"(3d’r?+dr+1). The optimum value of dr is 2.46. From (46)—(47)—(48) one finds 
Ground state H, HD 2 
PantQnnatSnn 114 85 57 cm"! 


The corresponding values for H.+, HD*t, D2* furnished by (57) are, respectively, 59, 44, 29 cm™!. 
In the case of He, the apportionment of 114 between Pan, Qna, and Sn» is, respectively, 23, 20, and 
71 cm; Eq. (41) would give 69 cm! for S,,. (In HD, the apportionment between Pan, Qan, Sun 
is different, but the total effect is inversely proportional to the reduced mass M, as one sees by 
adding (46), (47), (48).) To obtain the total shift, one must add the anharmonic correction obtained 
from (56), and, in the upper state, the effect of the outer electron as given by (55). The corresponding 
isotope shifts for He—HD are shown in the following table. It is seen that the largest contribution 
comes from (46)—(47)—(48)—(57), i.e., from the S,, terms, while (55) alone does not even give the 
right sign 


H.—HD Calculated Eq. (55) (56) (46)—(47)—(48)-(57) Total H2—HD Obs. 3(H2— Dz) Obs. 
2p MI —1s 15 +7.4 —2.7 —13.8 —9.1 —8 cm —12 cm™ 
2p 12 —15 15 +9.7 —0.8 —13.8 —4.9 —12 — 2 


The discrepancies here may be due to the very considerable experimental error, estimated by 
Dieke as 5 to 10 cm~, and to the inaccuracies of the Wang and Finkelstein-Horowitz wave func- 
tions.*® The 2p '2 state has an equilibrium distance 1.31 10-8 cm whereas 2 'II has almost exactly 
that 1.06 10-* of H.*+. Our H+ treatment of the inner electron should thus apply better to 2p 'I], 
and possibly this circumstance explains the greater disagreement for the transition involving 2!p '¥. 
At any rate, the observed shifts are qualitatively understandable. In the table we have also included 
half the observed displacement between Hz and De». This should equal the H2—HD shift according 
to our formulas, and doubtless more generally. The discrepancy between the measured value for 
H:—HD and half that for H,—Dz is hence indicative of a considerable experimental error. 








THE MEAN SQUARE ANGULAR MOMENTUM OF 
THE NORMAL He MOLECULE 


In connection with the diamagnetism of He, as 
well as with the determination of nuclear spin 
by means of Stern-Gerlach experiments on the 
hydrogen molecule, it is necessary to evaluate the 
integral Z,, defined in (12). Miss Frank and the 
writer obtained 0.394 with the Wang wave func- 
tions (cf. Eq. (52)). Wick*! obtained =,,=0.02 
with other wave functions De~‘"1t%+"+r2) | and 
suggested that our value was only an upper 
limit. Eq. (39) makes it clear why Wick obtained 
such a low result. The type of wave function 

279B. N. Finkelstein and G. E. Horowitz, Zeits. f. 
Physik 48, 118 (1928). 

39In applying (56), we use the values of a, a2 corre- 
sponding to the values of a, x given in Miss Sponer’s 
tables. Uncertainty in these constants may be a subsidiary 
source of error. The experimental value for the Ho-HD 
shift used in the table is by Jeppesen, as revised by Dieke;"4 
that for H:—Dz» is furnished by very recent work of 


Jeppesen (paper at Washington meeting, 1936). 
3G. C. Wick, Zeits. f. Physik 85, 25 (1933). 


which he used makes P,, vanish at r=~, 
whereas the Wang wave function gives the 
correct limit. Hence we do not believe that our 
value is necessarily much, if any too high. It is 
well below the asymptotic value 4b?79? charac- 
teristic of large r, which yields Z,,=0.9 at the 
actual internuclear distance. 

In adapting calculations on the Stern-Gerlach 
effect to HD rather than He, allowance must be 
made not merely for the altered reduced mass, 
but also for the somewhat different value of En» 
(cf. Eq. (52)), which arises because the center 
of gravity no longer coincides with the nuclear 
midpoint. In calculations on diamagnetism, 
however, it is not necessary to consider this 
distinction, as the formula for the susceptibility 
is invariant of the origin. 

The writer is indebted to Professors F. H. 
Crawford and R. S. Mulliken for valuable dis- 
cussions. 
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Viscosity and Diffusion Coefficients of Atomic Hydrogen and Atomic Deuterium* 


I. Ampur, Massachusetts Institute of Technology 
(Received March 16, 1936) 


The viscosity of atomic hydrogen has been recalculated from Harteck’s experimental values 
for the relative viscosity of atomic-molecular hydrogen mixtures, on the assumption that atomic 
hydrogen like molecular hydrogen may be regarded as a van der Waals gas. The numerical re- 
sults confirm this assumption and permit the evaluation of diffusion coefficients for H-H2, D-Dz, 
H-D», D-H, and D-H binary mixtures as well as the viscosity of atomic deuterium. 





N reaction kinetics investigations involving 
atomic hydrogen, it is often necessary to know 
the value of Dj2, the diffusion coefficient of atomic- 
molecular hydrogen mixtures, in order to esti- 
mate the extent of a wall reaction (determined by 
the rate at which atoms diffuse to the wall). In 
dynamic systems, the rate at which atoms 
disappear by diffusion from a given volume 
element must be considered as well as the rate at 
which atoms are being removed by pumping. A 
knowledge of the value of Dj» is required for such 
calculations. 

It is the purpose of the present communication 
to reevaluate the existing data for the viscosity 
of atomic-molecular hydrogen mixtures in order 
to obtain the value of ;, the viscosity of atomic 
hydrogen, from which Dj. may be calculated. 

Harteck! has determined the viscosities of 
atomic-molecular hydrogen mixtures of known 
composition by applying Poiseuille’s equation to 
the flow of such mixtures through a small tube. 
His measurements, in which molecular hydrogen 
was used as the reference gas, yielded relative 
values, nmix/n2z, at —80°, 0° and 100°C. In 
calculating the relative viscosity of atomic 
hydrogen, 7/n2, from his experimental results 
Harteck used a formula derived by Puluj? which 
gives the relative viscosity of one component of a 
binary mixture in terms of the relative viscosities 
of the other component and of the mixture, and 
the composition. From the relative values, 71/12, 
obtained in this way, Harteck attempted to 
evaluate 4; at the three temperatures from the 
known temperature variation of m2. If the vis- 
cosity of atomic hydrogen varies with tempera- 





* Contribution from the Research Laboratory of Physical 
paty Massachusetts Institute of Technology, No. 


1 P. Harteck, Zeits. f. physik. Chemie 139, 98 (1928). 
*J. Puluj, Wien. Ber. 79, 97, 745 (1879). 


ture according to the Sutherland’ formula we 
have the relation 


1+(C/273.16)(T/273.16)! 
1+C/T 





Np = No 
=6(T/273.16)!, (1) 


where 77 is the viscosity of atomic hydrogen at 
the absolute temperature T and 7;, that at 
0°C (273.16°A). C is the Sutherland constant for 
atomic hydrogen. Harteck found that at the 
three temperatures ® had the same value within 
his experimental error and concluded that the 
value of C was extremely small. He appears, 
however, to have forgotten the (7/273.16)! 
factor since he gives for the value of 7; at 0°C the 
mean of the three ® values. One gets the impres- 
sion from this that the viscosity of atomic 
hydrogen does not vary with temperature whereas 
inclusion of the (7'/273.16)! factor shows that m 
actually does increase with temperature. 

It has been considered advisable to recalculate 
Harteck’s results since, in addition to the error 
indicated above, Puluj’s formula does not give 
results which are in agreement with experiment. 

In the present calculation it will be assumed 
that both atomic and molecular hydrogen are 
van der Waals gases (gases whose molecules are 
hard elastic spheres of invariable diameter with 
attractive forces varying as an inverse power of 
the distance). This is equivalent to the assump- 
tion that the viscosities of these gases can be 
represented by the Sutherland formula. For 
molecular hydrogen this is actually the case at 
the temperatures under consideration. Formulas 
for the viscosities of the mixtures and of the pure 
gases considered in these calculations will be 
consistent with the above assumption. In addi- 


3 W. Sutherland, Phil. Mag. 36, 507 (1893). 
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tion, the recombination of atomic hydrogen to 
molecular hydrogen need not be considered since 
the rate of recombination is negligibly slow 
compared to the rate of momentum transfer 
which determines the viscosity. 

Chapman‘ has obtained an expression for the 
viscosity of a binary gas mixture which, ex- 
pressed in its general form, is independent of the 
laws of force governing the attraction and 
repulsion of the molecules. The derivations of 
Chapman have been subject to corrections by 
Enskog® giving the expression® 


Ax*?+ B2x(1—x)+C(1—<x)? 


, @ 
Ax?/m,+D2x(1—x)+(C/n2)(1—x)? ™ 





Qmix >= 


in which 
x=mole fraction of component 1 (atomic hydrogen), 

1 —x=mole fraction of component 2 (molecular hydrogen), 
A=1+kM,/M2, (3) 
C=1+kM2/M,, (4) 

Lg? 4 
B=(1-—k +5-(—+-), (5 
) 2G\m nz ) 
_ (M+ M2)P?kG 1 
- 2M; M, 2Gnin2 


1 (M,+ M2) PD,» 
where -= ’ 
G 2 RT 








(7) 


Dy» is the diffusion coefficient (at the pressure P) 
for the mixture, R, the gas constant, M, and Mo, 
the molecular weights of the respective com- 
ponents. 

For elastic spheres with an attractive field (van 
der Waals molecules) 


pit 2eu/3T 


k=0.60 , 
1+C\2/T 


(8) 


where Cio, the Sutherland constant for the 
mixture, is given by™* 


C = 7v(CiC2)} 
. [4((d:/de)}+(de/dy)!) p> 


v is the exponent in the expression for the 





(9) 


4S. Chapman, Phil. Trans. Roy. Soc. A211, 433 (1912); 
A216, 279 (1915); A217, 115 (1917). 

5 D. Enskog, (a) Inaug. Diss. Uppsala, 1917; (b) Arkiv. 
for Mat., Astron. och Fysik 16, No. 16 (1921). 

°M. Trautz, Ann. d. Physik (5) 18, 816 (1933). 
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attractive force, F=a/r’, and y a constant which 
is a function of v having the following values: 

4 5 7 9 
1.073 1.028 0.968 0).929. 


v 3 
BY 1.140 


C, and Cy, are the Sutherland constants for the 
pure gases and d, and de, the respective molecular 
diameters. In order to use the above expressions 
in actual calculations it is necessary to know the 
value of Dj». As a first approximation Enskog* 
gives a formula in which it is considered as being 
independent of the composition of the mixture. 
Actually the value changes slightly with compo- 
sition. If we let the approximate value be D,, 
then 
D2=Dof (x) = Dof'(a), 


where a is the degree of dissociation of molecular 
hydrogen into atoms, in which case x = 2a/(1+a) 
and 1—x=(1—a)/(i+a). The actual depend- 
ence of Dj. upon composition is quite small. 
Enskog** has evaluated f(x) for the case of 
elastic spheres without attractive forces and 
finds that the calculated dependence is confirmed 
by experiments on real gas mixtures indicating 
that f(x) is not appreciably different for a gas 
of van der Waals molecules than for a gas of 
elastic spheres without attractive forces. Ac- 
cording to Enskog, when (1—x)/x=0(a=1) 





M;? ) 
2(6M ?+8M,M.+15M,’) 
= 1.0061), 


Drs= Do 1+ 


and when x/(1—x) =0(a=0) 





Dr=Do( 1 + 


M,? ) 
2(15M,?+8M,M.+6M,?) 
= 1.0364D>». 


Thus, the value of Diz. changes by only three 
percent over the entire composition range. In 
view of this small variation it will be sufficiently 
accurate for the present purpose to let 


f' (a) =a+ba, 


where a and 6, evaluated from the fact that 
f’ (a) =1.0016 for a=1 and 1.0364 for a=0, are 
equal to 1.0364 and —0.0303. Using Enskog’s 
expression for Dy for a van der Waals gas 
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mixture we obtain finally 
3R?T3 
8(27)4N((ditd2)/2)?(1+Ci2/T) 





Dy2= 


1 1,31 
x —+—) 1.0364 — 0.03030), (10) 


M, Me 


where N is Avogadro’s number. Upon substituting 
the value of Dj. from Eq. (10) into Eq. (7) we 
obtain 

3(M,+ M2)!R'T?(1.0364 —0.0303a) 


ma di+de \2 
16(2n)!(atie)'n(——) (14 
2 


To evaluate 1/G it is necessary to know the 
values of d; and dy. The quantity d. may be 
computed directly from measured values of np. 
In order to determine d,; we shall assume as a 
first approximation that we are dealing with 
elastic spheres without attractive forces and 
solve for 7; (using Harteck’s measured values of 
Nmix/72). For this calculation it is possible to use 
rigorous expressions (analogous to Eq. (2)) 
derived by Enskog.** In such expressions Cj. 
obviously does not appear so that by a series of 
successive approximations it is possible to 
calculate consistent values of 4; and d;. In this 
manner it is found that at 0°C, d2=2.7110-8 
cm and d,;=2.49X10-8 cm. The value of ne in 
this calculation is given by the relation 


5.882 X10°T? 
1+35.8/T 





(12) 


obtained by F. G. Keyes (unpublished data) by 
applying the method of least squares to the 
existing viscosity data for hydrogen. Eq. (12), 
however, is the form assumed by the viscosity 
expression for attracting elastic spheres. The 
rigorous expression given by Enskog** is 


mee) 1 
TT 


. (13) 
16 Nd2(1+C/T) 





n= 


As a second approximation in determining d; we 
may determine d, from Eq. (13) (in this case, the 
diameter does not change with temperature), and 
assume that the ratio d,/dz is the same for 
attracting elastic spheres as for nonattracting 
elastic spheres. This assumption is equivalent to 
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the assumption that the Sutherland constant for 
atomic hydrogen is essentially the same as that 
for molecular hydrogen. We may check our 
approximate value of d; (as well as C,) when we 
have finally calculated 7; from Eq. (2), by means 
of Eq. (13). This almost amounts to an inde- 
pendent check since Eq. (2) defines the compo- 
sition dependence and Eq. (13) the temperature 
dependence of the viscosity coefficient. 

In Keyes’ expression for y2, the Sutherland 
constant, 35.8 is much smaller than the values 
found in the literature which vary from 70 to 83. 
There are two reasons for preferring 35.8 as the 
correct value. First, this value was obtained by 
the method of least squares applied to over 50 
viscosity values in the temperature range — 200° 
to +300°C, whereas the high values in the 
literature are calculated from a small number of 
viscosity values over a narrow temperature 
range. Second, the relation between the Suther- 
land constant and equation of state data clearly 
indicates that the small value is more reasonable. 

We may write the equation of state for a van 
der Waals gas at not too great pressure in the 
form’ 


P(V—Bo)=RT, 


where By is given by the expression 


By= -24N f (e-PolkT 1 )p%dr, (14) 
0 


In evaluating the integral it must be remembered 
that under the assumption of a van der Waals gas 
the intramolecular potential energy ‘‘Pot’’ be- 
comes infinite on contact of two molecules 
making the integrand zero for all values of r from 
0 to d, the diameter of a single molecule. A series 
integration of Eq. (14) gives the result 


2 2raN 
Bo=-1Nd*— 


+ (15) 
3 18d*kT 


if the potential energy is represented by —a/6r‘° 
in accordance with an attractive force a/r’. 
Since 


By=B—A/RT (16) 


we have the relations 


7F. G. Keyes, Chem. Rev. 6, 175 (1929); F. G. Keyes, 
Li _ and H. T. Gerry, Proc. Am. Acad. 70, 319 
(1936). 
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TABLE I. 








m2 X106 1/G X10 k A B , D m x 108 
69.0 84.4 0.572 1.286 1.040+42.210-°/m 1.525 X104+0.612/n, 65.3 
86.0 104.8 0.579 1.290 1.030+52.410-°/m 1.243 X104+0.609/n, 69.4 

103.7 126.1 0.584 1.292 1.024+63.010-°/m 1.042 X104+0.608/n, 83.2 


T(°A) — amix/22 a 

193.16 0.966 0.650 
273.16 0.883 , 0.637 
373.16 0.886 0.591 











1 za N? 
A=- , 

9 d 
where higher series terms in Eqs. (15) and (16) 
have been neglected. To evaluate C in terms of 8 
and A we make use of the relation derived by 
Enskog** 


2 
B=-7rNd}, 
3 


25Na 
C=—————-, 
3(v—1)d’"!R 


where, for the present case of »=7(F=a/r’), 
6=4X0.5165, so that 


C=0.5165Na/3X6Rd°=2.1X10-°A/B. (17) 


For hydrogen the values of A and @ are 0.2001 
X10" ergsXcc per mole and 20.96 cc per mole, 
respectively,® so that C calculated from Eq. (17) 
has the value 20.0. 

From Eq. (12) we find that 72.=86.010-° 
poises at 273.16°A, so that dz calculated from 
Eq. (13) is 2.55 10-8 cm. The assumed value for 
d,, therefore, is (2.55 2.49/2.71) X10-8 cm or 
2.34X10-* cm. In order to evaluate Cj. from 
Eq. (9) C,; is determined from the relation 
indicated by Keyes,’ namely, that for nonpolar 


gases 
= KP», 


where P> is the molecular polarization and K, a 
constant which is independent of the gas. For 
molecular hydrogen P»)=2.00'° while for atomic 
hydrogen the calculated value is 1.68" so that 


C,=35.8 X1.68/2.00 = 30.1. 


Substituting the proper values in Eq. (9) we find 
that Ci.=31.6. We are now in position to 
calculate 7, from Eq. (2). Table I contains the 


8J..A. Beattie and O. C. Bridgeman, Proc. Am. Acad. 
63, 229 (1928). 

°F. G. Keyes, Zeits. f. physik. Chemie, Cohen Festband, 
709 (1927). 

1° This value is a weighted mean of many experimental 
values derived from dielectric constant measurements. 
References are given in 
Electric and Magnetic Susceptibilities (Oxford Press, 1932), 


p. 67. 
u J. H. Van Vleck, reference 10, p. 205. 


J. H. Van Vleck’s Theory of 


results of this calculation as well as the values of 
the constants calculated from Eqs. (3), (4), (5), 
(6), (8) and (11). The values of nmix/n2 and @ are 
taken from Harteck’s paper.' 

It is now possible to use the calculated values 
of m1 to determine C; and d,; from Eq. (13) asa 
check upon the values assumed for use in Eq. (2). 
For this purpose only the 7; values at 273.16° 
and 373.16°A will be used, since the value at 
193.16°A may be in considerable error as pointed 
out by Harteck.! At this temperature, the 
capillary used in the experiments became coated 
with a sheath of ice because of the fact that the 
atomic-molecular hydrogen mixture contained 
water vapor to poison the walls against recombi- 
nation. In addition, Harteck states that at this 
low temperature the enhanced recombination 
prevented the gas mixture from assuming the 
temperature of the surrounding bath (—80°C). 
Both effects would tend to make nmix/n2 (and 
consequently, 7;) too high, in agreement with the 
high experimental value, 0.966. It is estimated 
that this value may be in error by as much as 
ten percent. Solving Eq. (13) simultaneously for 
C, and d; by substitution of the last two 7; values 
of Table I, we find that C;=28.2 and di=2.42 
X10-§ cm, in good agreement with the pro- 
visionally assumed values 30.1 and 2.34X10°° 
cm. It is now possible to make a second calcu- 
lation of 7; values from Eq. (2), using the values 
of C; and d,; obtained as a result of the first 
calculation. These approximations may be re- 
peated until subsequent variations in C, and 4: 
produce changes of less than one percent in the 
values of ;. This accuracy is at least as great as 
that of the experimental nmix/n2 values at 273.16° 
and 373.16°A. Table II summarizes the results of 
a series of such successive approximations, where 
approximation I is the calculation represented 
by Table I. Further approximations result in 7 
values that vary within the limits of the last 
three above approximations; accordingly, we 
may take as our final values for C; and d; those 
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TABLE II. 
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TABLE III. 








m X 106 m X 106 m X 106 
(193.16°A) (273.16°A) (373.16°A) 
65.3 69.4 83.2 
66.7 70.9 85.4 
66.1 70.2 84.3 
66.4 70.8 85.1 
66.0 70.3 84.5 
66.2 70.5 84.9 


ad, X 108 Ci 


30.1 
28.2 
34.0 
30.2 
$1.7 
31.4 


Approx. 





I 2.34 
II 2.42 
III 2.37 
IV 2.40 


\ 
VI 


2.38 
2.39 








values which most closely reproduce the mean of 
the last three sets of values (values at 193.16°A, 
of course, are not considered although they have 
been calculated for the sake of completeness). 
On this basis the best values are dj =2.39 10-8 
cm and C,=30.6, corresponding to 7; values of 
70.5X10-* and 84.7 KX10-°. 

On the basis of Harteck’s experimental results, 
therefore, we may represent the viscosity of 
atomic hydrogen from 0° to 100°C by 





4.746 X10-°T! 
1+30.6/T 


In like manner, we may calculate the diffusion 
coefficient for a mixture of atomic and molecular 
hydrogen from Eq. (10). For dj=2.39X10-§ cm 
and C,= 30.6, C12 calculated from Eq. (9) is 31.9 
so that the expression for the diffusion coefficient 


becomes 
374.9T?(1.0364 —0.0303a) 


P(i+31.9/T) 


= 





where P is expressed in dynes/cm’. 

Kinetic studies on the relative recombination 
rate of atomic hydrogen and atomic deuterium” 
indicate that the diameters and attractive forces 
are equal for hydrogen and deuterium atoms. 
The coefficient of viscosity for atomic deuterium, 
therefore, should be v2 times that of atomic 
hydrogen because of the greater mass of the 
deuterium atom, or 


6.711 10-73 
 1430.6/T | 





In the case of molecular hydrogen and deuterium 
this V2 ratio for the coefficient of viscosity has 
been experimentally verified"* indicating in this 


” 1, Amdur, J. Am. Chem. Soc. 57, 856 (1935). 
4°71, Amdur, J. Am. Chem. Soc. 57, 588 (1935). 


MIXTURE DIFFUSION COEFFICIENT 





265.17 4(1.0364 —0.0303x) 
P(1+31.9/T) 

342.27 4(1.0559 —0.0541x) 
P(1+31.9/T) 

306.174(1.0172 —0.0x) 

P(1+31.9/T) 

400.47 1(1.0364 —0.0303x) 

P(1+30.6/T) 





D-D, 





H—D, 





D—H: 


H—D 











case, also, an equality of molecular diameters and 
attractive forces for the two isotopes. With this 
information it is possible to calculate the 
diffusion coefficients for all possible binary mix- 
tures of atomic and molecular hydrogen and 
deuterium by substituting the proper values for 
M, and M, in Eq. (10) and in the Enskog 
expressions giving the concentration dependence 
of the diffusion coefficient. The result of such 
substitution is given in Table III. In each case x 
represents the mole fraction of the first (lighter) 
component. Since the concentration dependence 
is small, it makes little difference that in the 
above expressions a, the degree of dissociation, 
has been replaced by x, the mole fraction. The 
quantity, a, was used previously since slightly 
better agreement with Enskog’s rigorous elastic 
sphere formula was obtained thereby. In view 
of the fact that we have assumed that the 
concentration dependence for our van der Waals 
gas is the same as that for a gas of hard elastic 
spheres, and since the use of a or x gives values 
within 0.2 percent on either side of those obtained 
from Enskog’s formula, it seems desirable to 
make Eq. (10) general by the use of x. 
Although Harteck’s relative viscosity values 
furnish an experimental basis for the present 
expressions only from 0° to 100°C, it seems quite 
probable that use of the formulas beyond this 
temperature range will give fairly accurate (one 
percent or better) values for the viscosity and 
diffusion coefficients. This prediction is based on 
the fact that the Sutherland expression for 
molecular hydrogen (with the corrected value 
of C) holds from —200° to +300°C and that 
atomic hydrogen may be regarded as a van der 
Waals gas with a diameter and attractive force 
almost equal to that of molecular hydrogen. 
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The Free Energy of Iodine and Hydrogen Iodide from Spectroscopic Data 


G. M. Murpny, Chemical Laboratory of Columbia University 
(Received March 30, 1936) 


Partition functions have been calculated for HI, I. and I 
atom between 298.1 and 1500°K from spectroscopic data 
and the results have been fitted to a free energy equation 
of the usual thermodynamic form which may be differ- 
entiated and integrated to give change of heat content, 
specific heat and entropy for chemical reactions involving 
HI, Is, I and Hy. The experimentally determined equi- 


librium constants for the reactions I,—~2I and HI—}H, 
+I, have been compared with the theoretical equi- 
librium constants. For the latter reaction, it is shown that 
the results of Bodenstein may be in error due to incorrect 
temperatures. Reliable values have been given for the 
heat of dissociation of iodine and hydrogen iodide at 
ah. 





ALCULATION of thermodynamic functions 
from spectroscopic data may be made with 
complete confidence since the statistical weights 
of the degenerate energy levels are now well 
understood. Moreover, the measurements of 
molecular spectroscopy are so exact that in most 
cases these calculated functions will far exceed 
the experimental ones in precision. The requisite 
partition functions have now been summed and 
discussed for a large number of molecules! 
including Cl. and HCl,? Bre and HBr? but the 
corresponding data for I, and HI‘ have not been 
given in detail. 

In connection with experiments now com- 
pleted in this laboratory on the rate of decompo- 
sition of deuterium iodide, it became necessary 
to calculate the equilibrium concentration of DI 
at various temperatures and we have therefore 
obtained the partition functions for I,, HI and I 
atom and determined the free energy decrease 
and other thermodynamic functions for the 
reaction 


(1) 


and the corresponding deuterium reaction. This 
paper is concerned with the hydrogen case while 
the corresponding data for DI as well as a 
discussion of the rates of the two reactions will 
appear later. 


HI-}H2+ 31s 


1 For reviews of the subject, see Giauque, J. Am. Chem. 
Soc. 52, 4808 (1930); Zeise, Zeits. f. Elektrochemie 39, 
758, 895 (1933); 40, 662, 885 (1934). The latter author 
has compiled the available data for both diatomic and 
polyatomic molecules. 
as and Overstreet, J. Am. Chem. Soc. 54, 1731 

’ Gordon and Barnes, J. Chem. Phys. 1, 692 (1933). 

4 The free energy function for HI at a few temperatures 
has been given by Urey and Rittenberg, J. Chem. Phys. 
1, 137 (1933). Zeise (reference 1) has given the same 
functions for I, but has not discussed the experimental 
data critically as we do here. 


The method of calculating thermodynamic 
functions from spectroscopic data is well known 
and reference may be made to the reviews cited 
above for complete details. The partition func- 
tions are given by Q= Zpexp (—E/kT), where 
p is the statistical weight of the energy level and 
summation is made over the required number of 
electronic, vibrational and rotational states. The 
energy, as obtained from spectroscopic data is a 
complicated function of v and J(J+1), the 
vibrational and rotational quantum numbers, 
respectively. Once the Q functions and their 
derivatives are known the free energy, entropy 
and other properties follow by the well-known 
methods of statistical mechanics and thermo- 
dynamics. The laborious task of summing over 
all the energy levels has been simplified by 
several methods. We have found the procedure 
of Kassel® useful, especially for such a molecule 
as I, where a large number of vibrational levels 
must be summed. 

The molecular constants for I, were taken from 
Jevons® and may be found in Table I along with 
those for HI, which were recalculated from the 
recent experiments of Nielsen and Nielsen,’ who 
measured the fundamental and first harmonic 


TABLE I. 








I, HI 


Bo 3.73 K10™ 6.4576 

Do —4.536 10-9 —2.1013 10° * 
a 12 x0 line XI" 
B 1.63 X10- 
We 214.26 

X We 0.592 





2309.55 
39.735 








5 Kassel, J. Chem. Phys. 1, 576 (1933); Phys. Rev. 43, 
364 (1933). 

6 Jevons, Report on Band-Spectra of Diatomic Molecules 
(University Press, Cambridge, 1932). 

7 Nielsen and Nielsen, Phys. Rev. 47, 585 (1935). 
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TABLE II. 








Q —(F° —Eo°)/T 


I, 
4,375.64 
10,356.64 
24,109.80 
36,595.51 
79,854.51 
HI 
32.6307 
54.6738 
89.0072 
114.0486 
187.4235 


I 





54.193 
58.474 
62.487 
64.425 
67.988 


42.418 
46.012 
49.315 
50.916 


1500 53.916 


38.225 
40.794 
43.128 


44.237 
46.250 


298.1 

500 

800 
1000 
1500 








absorption bands of HI in the infrared. Their 
results are in agreement with the work of 
Czerny,® Salant and Sandow® and Kirkpatrick." 
In accordance with the practice of Giauque and 
co-workers, we have used the natural constants 
given in the Jnternational Critical Tables. The 
partition functions were calculated for both HI 
and I, at 298.1, 500, 800, 1000 and 1500°K and 
are given in Table II. Many intermediate values 
were also summed in order to test the interpola- 
tion formula which is described below. For the 
iodine atom, the normal state is an inverted ?P 
level, with a doublet separation” of 7600 cm~ 
and the contribution from this level is negligible 
up to several thousand degrees. The partition 
function for I atom is therefore independent of 
temperature in this range and is Q=4. 

From the Q functions, the free energy for the 
molecule or atom may be obtained from the 
expression 


— (F°— E,°)/T = (3/2)R In M+(5/2)R In T 


+R In Q—7.267, (2) 


where M is the molecular or atomic weight and 
E,° is the energy at the absolute zero. Following 
the usual custom, which is valid for temperatures 
above 25°C, the nuclear spin contributions have 
been excluded but for the homo-nuclear molecule, 


* Czerny, Zeits. f. Physik 44, 235 (1927). 

9 Salant and Sandow, Phys. Rev. 37, 373 (1931). 

Kirkpatrick, Phys. Rev. 49, 104 (1936). 

"Int. Crit. Tab., Vol. I, p. 16. The numerical constant 
of the free energy equation was taken as 7.267; R=1.9869; 
he/k =1.4324. 

” Turner, Phys. Rev. 27, 397 (1926). 
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I,, the Q function as calculated must be divided 
by two which has been done for the results 
given in Table II. The atomic weights for 
hydrogen and iodine were taken as 1.00778 and 
126.932, respectively.'* 

Since interpolation for other temperatures is 
fairly laborious, the freeenergy functions obtained 
at the five temperatures mentioned above were 


fitted to an equation of the form 
—(F°—E,°)/T=a/T+6 log T+cT+dT?+1, (3) 


which is suggested by thermodynamic considera- 
tions.'4 The constants in this equation are given 
for HI, I, and I atom as well as for He, the latter 
being taken from the free energy functions of 
Giauque!® at the same five temperatures as for 
the other molecules. The results will be found 
in Table III. It is probably not safe to extrapo- 
late these equations beyond 298.1 or 1500°K but 
over this entire range, the free energy function 
(3) is reproduced to within +0.001 cal./mole as 
tested by calculation and comparison with about 
15 points within this temperature range. 

If the value of AEZ,° is known for any particular 
reaction, it may be combined with equations of 
the form (3) for the individual reactants to 
obtain an equation of the same form, representing 
the free energy change of the reaction 


—AF°/T= — A( F°— E,°)/T — AE, / : i 


=RinkK, (A) 


where K is the equilibrium constant at unit 
pressure. In accordance with the thermodynamic 


TABLE III. 








He 


114.572 
16.7648 
—5.410 10° 
1.67 10" 
— 17.2832 


HI 


— 108.672 
14.6592 
7.649 x 10-4 
—4.0 107° 
6.2820 








I, I 
180.778 — 28.052 
19.8028 11.1292 

4.442 x10~4 1.907 X10" 
—6.87 x10-5 —4,09 x10 
4.4605 10.7277 








13 Aston, Mass-Spectra and Isotopes (Longmans, Green 
and Co., New York, 1933). ' 

14 See, for example, Lewis and Randall, Thermodynamics 
(McGraw-Hill Book Co., New York, 1932), particularly 
Chapter XV. 

18 Giauque, J. Am. Chem. Soc. 52, 4816 (1932). 
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notation, we write Eq. (4) in the form 


AF®/T=AH,/T— AY oI|n T—3AT,T 
—1ATT?+I. (5) 


This equation, except for the uncertainty in the 
quantity, AE,° will give the correct free energy 
change of the reaction to +0.001 cal./mole over 
the entire temperature range from 298.1 to 
1500°K. It may also be differentiated and inte- 
grated to give changes in heat content, specific 
heat and entropy 


AH = AH) + AMT +43AF,T2+1AT2T#+---, (6) 
AC, =ATo + AM THAT2T?2+- + -, (7) 
AS=AT)(In T+1)+AT,7+24MT?—J. (8) 


We have not investigated the precision of any of 
these expressions by direct summation except 
the entropy equation which we discuss later but 
it seems likely that they will all give results of 
sufficient precision for most applications, al- 
though they are not as exact as the free energy 
equation. Their chief advantage lies in the 
avoidance of the tedious process of summation 
which is required for determining the entropy 
and other quantities from spectroscopic data. 


THE DISSOCIATION OF IODINE MOLECULE 


The decomposition of molecular iodine into 
atoms has been studied by Starck and Boden- 
stein,'!® Braune and Ramstetter!? and DeVries 
and Rodebush.'*§ Gibson and Heitler!® have 
shown that the experimental results of Starck 
and Bodenstein are in agreement with the 
statistical calculations. They assumed the iodine 
molecule to be a rigid rotator and a simple 
harmonic oscillator and added a slight correction 
term to account for anharmonicity but the exact 
partition functions as calculated here do not 
change their results by a very large amount. 
The value of AE,° for this reaction has been 
given by Brown” from extrapolation of spectro- 
scopic data as 1.535+0.001 ev or 35,395+23 
cal./mole.*! This result when substituted into 

1®Starck and Bodenstein, Zeits. f. Elektrochemie 16, 
961 (1910); 22, 327 (1916). 

17 Braune and Ramstetter, Zeits. f. physik. Chemie 102, 
480 (1922). 

, me _— and Rodebush, J. Am. Chem. Soc. 49, 656 
' 19 Gibson and Heitler, Zeits. f. Physik 49, 465 (1928). 


20 Brown, Phys. Rev. 38, 709 (1931). 
21 We take 1 ev = 23,085.5 cal./mole. 


MURPHY 


Eq. (4) is definitely inconsistent with the 
equilibrium data. Gibson and Heitler have 
rejected a similar result of 1.532 ev given by 
Kuhn” as too small and have preferred to take 
1.544+0.003 ev by extrapolation of the bands of 
Mecke.”* Because of this uncertainty, we have 
used Eq. (4) and the experimental equilibrium 
constants to determine AEF,°, giving a result of 
35,795+100 cal./mole. While this is definitely 
higher than any of the previous values it seems 
necessary if all of the experimental data are to 
be given equal weight. Combining this value of 
AE,° with the proper constants from Table III 
and using Eq. (5) we obtain values for the 
coefficients of that equation as given in Table 
IV. A plot of log K against 1/7 for all the 


TABLE IV. 








I,—>2I 


36.032 
1.0664 
— 1.256x10-4 
— 7.86 x10-§ 
— 16.9949 


HI—}H2+ 31. 


715.48 
1.5741 
— 1.6266 x10°° 
2.973 1077 
12.6934 











experimental data is given in Fig. 1, which is 
seen to be a straight line. A more certain value 
of AE,° will only result in a shift of the whole 
curve parallel to itself in either direction. The 
equilibrium constant is given by 


K,=4a°P/(1—a*), (9) 


where a is the degree of dissociation of iodine 
and P is the pressure in atmospheres. 

Using Eq. (6), we obtain for the heat of 
dissociation of iodine into atoms at 298.1°K 


AHo9s.1 = 36,343 +100 cal./mole 
and from Eq. (8) 
ASogg.1 = 24.097 cal./degree per mole. 


By direct summation at this temperature, we 
have obtained So9g.1(I2) = 62.29; Soos.1(1) = 43.19; 
ASo9g.1= 24.09 in satisfactory agreement with 
Eq. (8). The value for I; is in exact agreement 
with the calculations of Giauque,?4 who used an 
approximate method for determining the parti- 
tion functions. 
*% Kuhn, Zeits. f. Physik 39, 77 (1926). 


23 Mecke, Ann. d. Physik 71, 103 (1923). 
24 Giauque, J. Am. Chem. Soc. 53, 507 (1931). 
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Fic. 1. The equilibrium constant for the reaction: I,—>21. Theoretical curve from Eq. (5) and 
Table IV; experimental points from various authors. 


THE DISSOCIATION OF HYDROGEN IODIDE 


The thermal equilibrium between hydrogen, 
iodine and hydrogen iodide has been studied by 
Bodenstein®® and is generally considered to be 
the most carefully investigated gaseous equi- 
librium to be found in the literature. In addition 
to the earlier experiments of Hautefeuille?® and 
Lemoine,?’ the equilibrium has also been studied 
at 1022 and 1217°C by von Falckenstein?* and 
near room temperatures by Stegmiiller,?® using 
an electromotive force method. More recently, 
Urey and Rittenberg*® have studied the equi- 
librium with mixtures of hydrogen and deuterium 


*® Bodenstein, Zeits. f. physik. Chemie 13, 56 (1894); 
22, 1 (1897); 29, 295 (1899); Habilitationsschrift, Heidel- 
berg, 11-36, W. Engelmann, Leipzig, 1899. 

** Hautefeuille, Comptes rendus 64, 608 (1867). 

*7 Lemoine, Ann. chim. phys. (5) 12, 145 (1877). 

**von Falckenstein, Zeits. f. physik. Chemie 68, 270 
(1909); 72, 113 (1910). 

** Stegmiiller, Zeits. f. Elektrochemie 16, 85 (1910). 
(i934y enone and Urey, J. Am. Chem. Soc. 56, 1885 


and have shown that the constants obtained 
from the spectroscopic data for the reaction 


H.+2DI—D.+2HI (10) 


are in agreement with the experimentally deter- 
mined ones. The latter authors have also given 
two measurements with hydrogen iodide alone. 
For this reaction, the value of AF,° is still 
more uncertain than for the iodine decomposi- 
tion. We have therefore followed the same 
procedure as above to determine this quantity 
although we have used only the data of Boden- 
stein and the two points of Rittenberg and Urey. 
Following a suggestion of Haber,*! however, we 
first correct the measured equilibrium constants 
for the dissociation of iodine molecules into 
atoms which is not very great but should not be 
neglected at the highest temperatures. If & is 


the concentration of hydrogen or iodine at 


31 Haber, Thermodynamics of Technical Gas-Reactions 
(Longmans, Green and Co., New York, 1908). 
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equilibrium and a@ is the percentage decomposi- 
tion of iodine into atoms, the corrected equi- 
librium constant for the reaction (1) is given by 


K’=&(1—a)'/2(1—£)=K(1—a)}, (11) 


where K is Bodenstein’s constant. 

A more serious difficulty is the uncertainty of 
Bodenstein’s temperatures. His experiments 
cover the temperature range between 283 and 
508°C and several thermostats were used. A 
number of temperatures were obtained with a 
sulfur vapor bath, the temperature being varied 
by boiling the sulfur under reduced pressures. 
The temperature was determined by measuring 
the vapor pressure of the sulfur with a man- 
ometer, the pressure-temperature curve for sulfur 
having been previously determined by means of 
calibrated mercury thermometers. Comparison 
with the present vapor pressure curve” of sulfur 
shows some divergence from the temperatures 
as given by Bodenstein and since he gave the 
observed pressure in each case, it is possible to 
make corrections. The present vapor pressure 
curve of sulfur may be in error by as much as 3° 
at the lowest temperature but since Bodenstein’s 
temperature differs at this point by 10°, the 
corrected temperatures as given in Table V seem 
more nearly correct than those given by Boden- 
stein. The boiling point of mercury was used 
for another temperature and this is apparently 
correct while PS; was used in another thermostat 
and we have changed this temperature from 
508° to 514°C in accordance with the Jnter- 
national Critical Tables, although this is more 
uncertain than the sulfur points. All of his other 
temperatures were measured with a mercury 
thermometer and naturally nothing can be said 
of their errors, if any. The two points of Urey 
and Rittenberg mentioned above are apparently 
correct since the temperatures were directly 
measured with a thermocouple. The corrected 
temperatures and equilibrium constants cor- 
rected in accordance with Eq. (11) for the 
dissociation of iodine are given in Table V. 

Using these corrected temperatures only, we 
have then used Eq. (4) to calculate AK,° as 
before. The value obtained is 


AE,° = 971.83+413.65 cal./mole. 


® Int. Crit. Tab., Vol. III, p. 201. 
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TABLE V. 








TBod. log K’ (exp.) log K’ (calc.) 





—0.8121 
—0.8268 
—0.8377 
—0.8523 
— 0.8630 
—0.8736 
— 0.8826 
—0.8826 
—0.8925 
—0.9179 
—0.9361 
—0.9588 
—0.9765 


781 —0.8002 
760 — —0.8163 
— —0.8328 
716 —0.8513 
700 — 0.8618 
683 —0.8765 
— — 0.8840 
666 — 0.8876 
647 — 0.9004 
629 —0.9179 
601 —0.9350 
575 —0.9552 
556 —0.9639 








TBod, are Bodenstein’s experimentally determined temperatures in 
°K, TVeorr, are the temperatures which could be corrected as explained 


above. 


* Data of Rittenberg and Urey. 


This constant has then been used in Eq. (5) to 

give the free energy decrease of the reaction, 

the constants in the equation being found in 

Table IV. The same remarks that were made 

above as to precision apply also in this case. 
Using Eq. (6) we obtain 


AH7o9s.,5= 1113415 cal./mole. 


If we take thé heat of sublimation of iodine as 
7438 cal./mole as given by Giauque,”4 we obtain 
—6325+15 cal./mole for the heat of dissociation 
of HI into hydrogen and solid iodine. This may 
be compared with the calorimetric determina- 
tions of Thomsen and Berthelot** who found 
6000 and 6300 cal./mole, respectively. A recent 
calorimetric determination by Giinther and 
Wekua* for the reaction 
HI(g)+3Cle(g)HCl(g)+312(s) (12) 
gave AH =27,966+50 cal./mole at room temper- 
ture. Assuming the latter to be 25°C and com- 
bining with the heat of formation of HCl, given by 
Giauque and Overstreet? as 22,028+10 cal./mole 
gives 5942+50 cal./mole for the heat of decom- 
position of HI. This value cannot be combined 
with the spectroscopic data to give agreement 
with the experimental equilibrium constants” 
and it seems evident that some unsuspected source 


338 Lewis and Randall, reference 14, page 525. 

“4 Giinther and Wekua, Zeits. f. physik. Chemie A154, 
193 (1931). 

4% Zeise, reference 1, has also reached a similar conclusion. 
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Fic. 2. The equilibrium constant for the reaction: HI—}H2+ 312. Theoretical curve from Eq. 


14 
YT x1o7 


16 1.8 


(5) and 


Table IV; experimental points from various authors. 


of error was present in these calorimetric determi- 
nations. 
Using Eq. (8) we find 


ASoo3.1= — 2.623 cal./degree per mole. 


Summation gives for HI at the same temperature 
the value —49.36 and combining with S(He) 
= 31.23 from Giauque’’ and the entropy of I; 
as given above, we obtain ASoo3.;= —2.60, in 
satisfactory agreement with our interpolation 
formula (8). Giauque** using a rigid molecule 
for HI, obtained S(HI) =49.8 and since he used 
the same values as above for He and Is, he finds 
AS=—3.04. The difference is due entirely to 
his value for HI. Experimentally, Giauque and 
Wiebe, ** find Seos.1(H1) = 49.5. 

The equilibrium constants as calculated by 
Eq. (4) and the constants given in Table IV 
will be found in Table V and a plot of log K 
against 1/T in Fig. 2. We have also included 
there the results of von Falckenstein and Steg- 
miiller. The former were apparently incorrectly 


1929) 0 and Wiebe, J. Am. Chem. Soc. 51, 1441 


changed by the author for the decomposition of 
I, and the possibility of a reaction between Pt 
and I, was not considered as pointed out by 
Jellinek and Uloth,*? which undoubtedly ac- 
counts for the divergence from the calculated 
curve, although the values of K’ as plotted here 
do not help the agreement very much. The 
results of Stegmiiller are probably in error 
because of neglect of liquid junction potentials 
as stated by Lewis and Randall.** On the other 
hand, the points of Bodenstein and Rittenberg 
and Urey are in good agreement with the 
theoretical curve. Consideration of the errors 
involved shows that an error of 0.5 percent in 
the analysis would account for the divergence of 
the point which lies furthest from our theoretical 
curve. However, if one drew the best curve 
through all of the experimental points, the slope 
of the line would be quite different from the one 
we have chosen from the statistical calculations. 
Although Fig. 2 has a slight curvature and is not 
a straight line like Fig. 1, the spectroscopic data 
can certainly not be so adjusted to give the 


37 Jellinek and Uloth, Zeits. f. anorg. Chemie 151, 161 
(1926). 
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extreme curvature required by the experimental 
points but can only be shifted parallel to itself 
because of the uncertainty in AE,°. The slight 
curvature in Fig. 2, mentioned above is probably 
due to the anharmonic terms and the rotational 
stretching terms in the energy expression for HI 
which also explains why Fig. 1 is a straight line 
since these quantities are nearly negligible for a 
heavy molecule like iodine. 


KEITH BREWER 








Note added in proof: Further consideration shows that 
the attempted correction of Bodenstein’s temperatures 
may not be justified as his experimental temperatures with 
the sulfur baths fit the theoretical curve as well as the 
corrected ones. However, this does not affect any of the 
calculations made above since either set of temperatures 
will give the same AE,° within the probable error stated. 
The doubtful points are thus the ones at the lowest and 
highest temperatures which may be explained by errors in 
analysis or what is more likely some further uncertainty 
in the experiments. 








JUNE, 1936 










Abundance ratio measurements are necessarily compli- 
cated by the possibility of an isotope effect at the source. 
In the emission of positive thermions of lithium, potassium 
and rubidium the extent to which the ion ratio represents 
the atom ratio within the source can be estimated by ob- 
serving the change in the abundance ratio as a function 
of the percent of alkali evaporated. In the case of lithium 
an initial ratio of Li?7/Li®=11.60 was observed; the abun- 
dance ratio followed the computed ratio for ideal evapora- 






N accurate estimation of isotope abundance 
ratios from mass-spectrographic measure- 
ments necessarily involves the question of an 
isotope effect at the positive ion source. The 
extent to which the ion ratio represents the true 
abundance ratio is amenable to a comparatively 
simple test in the case where positive thermions 
are used as the ion source. 

In the mass spectrograph which was described 
recently,! heated platinum disks impregnated 
with alkali were used as the source of positive 
ions. These ion emitters were prepared by coating 
the platinum with a water paste made from the 
salt to be tested. These coated disks can be used 
directly as positive ion emitters, although plati- 
num impregnated with alkali is in general the 
more satisfactory source. This impregnation is 
readily accomplished by heating the disk to dull 
redness in contact with the adhering coating for 
about half an hour; the coating is then removed 
by scraping and washing with distilled water. 


1A. Keith Brewer, J. Am. Chem. Soc. 58, 365 (1936). 
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2 A. Keith Brewer, Phys. Rev. 47, 571 (1935). 
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tion until the lithium was 40 percent exhausted, beyond 
which it fell below the ideal value. The abundance ratios 
of K®/K"=14.22 and Rb®/Rb* =2.59 remained almost 
unchanged during the process of exhaustion. It is apparent 
therefore that a correction factor needs to be applied to 
the ion abundance ratio for lithium, but is unnecessary for 
potassium and for rubidium. The atomic weights com- 
puted from the abundance ratios are Li=6.939, K = 39.094 
and Rb =85.46. 


, 


During the heating process the platinum dis- 
solves sufficient alkali to become an excellent 
emitter of positive ions. 

The test for a preferential isotope effect in the 
emission of alkali ions was made by observing 
the changes in the abundance ratio as a function 
of the fraction of the alkali evaporated. This 
fraction was obtained from a determination of 
the time and current necessary for the complete 
exhaustion of the less abundant isotope; the 
abundance ratio was measured at regular inter- 
vals during the exhaustion period. 


RESULTS 


The procedure just described yields the abun- 
dance ratio as a function of the less abundant 
isotope evaporated from the platinum. Curve 4, 
Fig. 1 shows a representative run obtained with 
lithium. The initial abundance ratio, Li7/Li* 
= 11.6, is the same as that reported previously.’ 
During the exhaustion period it was necessary 
















ISOTOPE EFFECT IN EVAPORATION OF IONS 


LITHIUM 
A- Observed 
B- Computed for ideal Evaporation 


Percent Li® Evaporated 


Fic. 1. 


to raise the disk temperature to obtain sufficient 
current to measure the abundance ratio; the 
approximate temperatures at which the ratios 
were determined are shown accompanying the 
points. Temperature, in itself, appears to have 
little effect on the abundance ratio since it 
remained the same over the range of positive 
ion currents that could be investigated. 

The results obtained with potassium are shown 
in line A, Fig. 2. The abundance ratio, K*®/K* 
= 14.22, is essentially the same as that obtained 
from most mineral sources. This ratio remained 
almost unchanged throughout the evaporation 
period, dropping to 14.1 for about 90 percent 
emitted, and to 13.8 for the last traces of 
potassium left in the platinum. 

The initial abundance ratio obtained from 
rubidium is Rb*®*/Rb’?=2.59+0.01. This is 
essentially the same as 2.59+0.04 obtained 
previously? but is lower than the value of 
Rb®/Rb§’7=2.68+0.02 recently reported by 
Nier,* using a different type of ion source. 
Since the previous value was obtained with an 
entirely different spectrograph and rubidium 
salt, it seems definite that the abundance ratio 
from this type of source is very close to Rb*°/Rb*? 
= 2.59. 

The abundance ratio as a function of the 
percent of Rb*’? evaporated is shown in Fig. 3. 
No measurable change was observed during the 


(1934 
*A. O. Nier, Phys. Rev, 49, 272 (1936). 


034) Keith Brewer and P. D. Kueck, Phys. Rev. 46, 894 


POTASSIUM 


A- Observed 
B - Computed for Ideal Evaporation 


Percent K* Evaporated 


Fic. 2. 


exhaustion period. The conspicuous fact to be 
observed in the above results is that potassium 
and rubidium show almost no isotope effect 
while lithium shows a marked shift in the 
abundance ratio during the evaporation period. 


DISCUSSION OF RESULTS 


The thermionic emission of positive ions from 
platinum impregnated with alkali offers two 
possibilities for the exhibition of an isotope 
effect: First, through a difference in the rate of 
diffusion to the surface, and second, through a 
difference in the rate of escape from the surface 
as positive ions. 

Assuming that the alkali atoms are at liberty 
to diffuse freely through the platinum, the 
number N of atoms striking unit area of the 
bounding surface in unit time with a velocity 
component between u and u+du normal to the 
surface is 


N(u)du=n(m/2rkT) ue !?*T du, (1) 


where 1 is the concentration of alkali atoms per 
unit volume of platinum and m is their mass.° 

The total number of atoms striking the 
bounding surface is the integral of the above 
expression from 0 to «© while the number 
escaping as ions is the integral between % and 
o where u% is the minimum velocity for escape. 
Thus N+, the number of positive ions escaping 
per second per sq. cm is 


5L. B. Loeb, Kinetic Theory of Gases (McGraw-Hill). 
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RUBIDIUM 
nt 
2 A- Observed 
=> B - Computed for Ideal Evaporation 
« 
l 
0 10 20 30 40 50 60 70 80 90 100 
Percent Rb®” Evaporated 
Fic. 3. 
Nt+=n(kT/2am)be—m?/2k7, (2) 


This equation can be applied directly to the 
case of lithium, a similar equation being written 
for each isotope. This gives 


n;/N¢= N;z+/Net(m;/me) de—(meue?—mjuz?) /2kT (3) 


for the abundance ratio of the isotopes in the 
metal (;/m¢) as a function of the ratio of the 
positive ion currents (N;+/N,t). 

From this expression it will be seen that if 
perfect mixing of the lithium atoms takes place 
within the platinum the observed ion ratio must 
be multiplied by the square root of the ratio of 
the isotope masses to correct for the difference 
in the collision frequency with the bounding 
surface, and also by an exponential factor to 
correct for a possible difference in the work 
functions. 

While it cannot be stated with certainty that 
the work function for the two isotopes is the 
same, the difference must be negligible since the 
ionization potential is the same for each to 
within one-millionth of a volt. This being the 
case the exponential factor becomes equal to one 
and can be neglected. In consequence it seems 
apparent that the ratio of the isotope atoms 
leaving the surface as positive ions is proportional 
to the concentration at the bounding surface. 
Any isotope effect that may exist, therefore, 
comes from the difference in the collision fre- 
quency with the surface. 

Assuming the exponential factor in Eq. (3) is 
equal to unity and that perfect mixing takes 
place in the platinum disk used in these experi- 


KEITH BREWER 


ments it is possible to compute the abundance 
ratio that should exist for the escaping ions at 
different stages of the evaporation of the alkali. 
Lines B in Figs. 1, 2 and 3 are so computed for 
various percentages of the lesser isotopes evapo- 
rated. A comparison of the observed A lines with 
the computed B lines gives, in consequence, the 
extent to which perfect mixing and ideal diffusion 
take place within the platinum source. 

A possible explanation of the fact that lithium 
exhibits an isotope effect while potassium and 
rubidium do not may be found in the dimensions 
involved. The radius of the lithium ion is 0.6A, 
potassium is 1.33A, and rubidium is 1.48A; the 
radius of the platinum atom in the metal is 
1.39A. A solid solution of potassium or rubidium 
in platinum, since the dimensions are essentially 
the same, will, therefore, be of the replacement 
type in which the alkali atom replaces a platinum 
atom in the metal lattice. In this case diffusion 
takes place step by step from lattice to lattice 
with little or no opportunity for mixing; an 
isotope effect, in consequence, is not to be 
expected. The lithium atom, on the other hand 
is small compared to the separation of the 
platinum atoms; here interstitial solution is 
possible in which the lithium atoms do not enter 
into the metal structure but actually penetrate 
the platinum lattice. The lithium atoms are in 
this case free to move more or less at random 
through the platinum mass, thus permitting a 
more ideal mixing than is possible in.a replace- 
ment solution where the dissolved atoms are 
confined to fixed positions. This mixing is not 
perfect, however, since Fig. 1 shows that the 
computed and observed curves begin to diverge 
before the alkali is half exhausted. It has been 
pointed out that not all atoms colliding with the 
bounding surface will escape as positive ions, 
but only those of high energy content. This 
means that lithium atoms will tend to concen- 
trate in the surface region thus congesting the 
interstitial space and obstructing the freedom 
of collision. 

Results obtained with coated lithium alumi- 
num silicate emitters lead to a conclusion of 
similar character. The abundance ratio remained 
constant at 11.60 for some time, then rose to a 
final value of 13.88 after prolonged heating above 
the fusion point. 
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The various methods that have been employed 
for the separation of isotopes have been discussed 
by Aston.® The present results suggest some of 
the conditions under which thermionic emission 
can be utilized for this purpose. It is evident 
that when the atoms to be emitted as thermions 
do not have to conform to any definite lattice 
structure within the emitter conditions ap- 
proaching ideal evaporation may be realized. 
Emission from fused salts show an appreciable 
isotope effect. In the case of metallic emitters 
the positive ions which have been held in 
interstitial solution exhibit a preferential emis- 
sion of the light isotope ions while those held in 
replacement solution show little if any effect. 

The values for the abundance ratios which 
have just been presented may now be employed 
in estimating the atomic weights of these ele- 
ments by the method recently described in 
detail for potassium.? The atomic weight of 
39.094 obtained for potassium remains un- 


‘F. W. Aston, Mass-Spectra and Isotopes (Edward 
Arnold & Co., London, 1933). 
7 A. Keith Brewer, J. Am. Chem. Soc. 58, 370 (1936). 
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changed. With 1.00027 for the conversion factor 
from the physical to the chemical scale and —8.2 
for the packing fraction, the atomic weight of 
rubidium comes out to be 85.46. This value 
compares favorably with 85.44 accepted by the 
committee on atomic weights and 85.48 recently 
obtained by Archibald, Hooley and Phillips.*® 

In estimating the atomic weight of lithium it 
is necessary to apply a correction for the isotope 
effect of free evaporation to the abundance ratio 
obtained from fresh filaments. Assuming this 
correction is initially that for ideal evaporation 
the abundance ratio becomes 11.60¥ 7/6= 12.52. 
With this corrected value for the abundance 
ratio and 7.0146 and 6.0146 obtained by Bain- 
bridge® for the isotope masses the atomic weight 
comes out to be 6.939. The most recent values 
for the isotope masses’ of Li®=7.0180 and 
Li® = 6.0167 lead to 6.942 for the atomic weight. 
Both these values are in close agreement with 
the accepted value of 6.940. 

8 Archibald, Hooley and Phillips, J. Am. Chem. Soc. 58, 
70 (1936). 


9K. T. Bainbridge, Phys. Rev. 44, 56 (1933). 
10M. L. Oliphant, Nature 137, 396 (1936). 
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Quantum-Mechanical Treatment of Helium Hydride Molecule-Ion HeH* 


J. Y. Beacu, Gates Chemical Laboratory, California Institute of Technology 
(Received March 23, 1936) 


Helium hydride molecule-ion has been treated by the 
variation method, including ionic and polarization terms 
in the variation function and varying the effective nuclear 
charge. The equilibrium internuclear distance is found to 
be 1.57a». Fitting a parabola to the lower part of the best 
potential energy curve leads to a value of 0.43 megadyne/ 
cm for the force constant and a value of 2800 cm™ for the 


INGLY charged helium hydride molecule-ion 
has been known! for many years from mass 
spectroscopic investigations of helium hydrogen 
mixtures. Its position with respect to that of 
H*H?H!+ in the mass spectrograph has been used 
to determine accurately the mass of the deu- 
teron.? Spectroscopic work has not shown the 
existence of this ion. 
'F. W. Aston, Isotopes, Both editions; T. R. Hogness 


and E. C. Lunn, Phys. Rev. 26, 44 (1925). 
*K. T. Bainbridge, Phys. Rev. 44, 57 (1933). 


fundamental vibration frequency. The minimum energy 
obtained for the system is —5.844Wy, which, when sub- 
tracted from the best value for the energy of the helium 
atom using the same type of variation function, gives a 
value for the dissociation energy of 2.02 ev. The doubly 
charged ion, HeH**, is found to be unstable. 


The only previous quantum-mechanical dis- 
cussion of HeH? is a variation treatment made 
by Glockler and Fuller.* They consider two 
problems: (1) The interaction of a proton with 
an excited helium atom, the wave function for the 
electrons being ¥;=ae~*"! and Yo= (b+ cre)e—F"?; 
and (2) the interaction of a hydrogen atom with 
a singly ionized helium atom, the wave functions 


3G. Glockler and D. L. Fuller, J. Chem. Phys. 1, 886 
(1933). 
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for the electrons being y= ae~*"*! and yW2= be-8”2. 
When in the first case a and 6 are put equal to 
2 and 0.7, respectively, the interaction is re- 
pulsive at all distances. The functions y; and ye 
are the 1s and 2s hydrogen-like wave functions, 
the helium atom being approximately in the 
first excited state. In case two a and 8 are put 
equal to 2 and 1, the resulting potential energy 
curve has a minimum of 8.1 ev at 1.3a . In this 
calculation the unsymmetrical Sugiura-type 
integral was neglected. As pointed out by Glock- 
ler and Fuller, including this integral will raise 
the energy by an undetermined amount. 

The above calculation is not evidence for the 
existence of HeH+. The dissociation energy, 
8.1 ev, obtained by making the above approxi- 
mation, is with respect to the dissociation 
products hydrogen atom and helium ion, while 
the difference in the ionization potentials of 
hydrogen and helium is 10.9 ev. Under these 
conditions the molecule would be unstable with 
respect to dissociation into a helium atom and a 
hydrogen ion. To prove the stability of the 
molecule it is necessary to show that for some 
intermediate internuclear distance the energy of 
the molecule is less than the energy of the 
helium atom. 

The variational integral 


W= fV*HWVdr/ Sf V*Vdr 


gives an upper limit for the energy of a system. 
W is an approximate wave function. H is the 
Hamiltonian operator and in this case is 


h? 2e* 2e? 


i 





When a variation function of the form V= }-c,y, 


n=1 
is used, the condition that Wbea minimum with re- 
spect to the constants c, is that the determinantal 
equation |Hi;— Ai;W| =0 (Eq. (1)) be satis- 
fied,4 where H;;= SviHbjdr and A= S vibjdr. 
The approximate energy for the ground state 


4See for example Pauling and Wilson, Introduction to 
Quantum Mechanics. 
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will be the lowest root of this equation. The 
constants, c,, are obtained by solving the set of 
™ 


linear equations }0c;(H;;— Ai;W), i=1, 2, +++ m. 
j=1 


Most of the integrals required in evaluating J//,; 
and A;; were used by Rosen® and Weinbaum*® in 
discussing the hydrogen molecule. Dr. Weinbaum 
has not published his integrals but kindly placed 
them at my disposal. The integrals used in this 
calculation, not tabulated by Rosen,’ are given 
in the appendix. 

The variation functions used here are com- 
posed of the following functions: 


a(1 ) = (Z*/mao*) te—Zrai/a0 
b(1) =(Z3/may?) be 2701/0, 
a’ ( 1 ) = (Z5/mrao*)*r a1 cos 6,1;e7 2721/20 
and the same functions for electron (2). Each 
function is normalized to unity. a(1) is the 
ground state wave function for a hydrogenic 


atom of nuclear charge Z, and satisfies the 
equation 


h? Ze? 
(-- vit——Ja(l) = —2*Wya(1) 


822m al 


where Wy = (e?/2ao). a’(1) is a first excited wave 
function for a hydrogenic atom of nuclear charge 
2Z and satisfies the equation 


h? 2Ze? 
(- vit——Ja'(l) = -2°Wna (1), 


This form of a’(1) was used by Rosen® in intro- 
ducing polarization into the hydrogen molecule 
calculation. It is very convenient as it leaves 
the exponents in a(1) and a’(1) identical making 
the calculation of integrals much simpler. In 
the present calculation the same effective nuclear 
charge, Z, is used in each term of the variation 
functions. This is done to eliminate the necessity 
of calculating the unsymmetrical Sugiura integral 
and other integrals. This procedure was adopted 
by Pauling’ in treating the helium molecule-ion 


He,t. 


5 N. Rosen, Phys. Rev. 38, 2099 (1931). 
6S. Weinbaum, J. Chém. Phys. 1, 593 (1933). 
7 Linus Pauling, J. Chem. Phys. 1, 56 (1933). 
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It is evident that the stable state of the 
molecule is a singlet state so we shall require the 
variation function to be symmetric in the spatial 
coordinates of the two electrons (and anti- 
symmetric in the spin coordinates). In the 
following paragraphs we shall calculate the 
energy of the systems for several such variation 
functions. It will be found possible, by taking a 
linear combination of the lowest lying states for 
simple wave functions, to obtain a calculated 
value of the energy of the molecule that is 
lower than the experimental value for the helium 
atom. 


1. NorMAL HELIUM ATOM AND PROTON 


The most stable state for the system for a 
single wave function of the types we are con- 
sidering, Eq. (2), is the state helium atom and 
proton. It will therefore be necessary, if the 
lowest value of the energy is desired, to include 
this state in the calculation. Taking as the elec- 
tronic wave function of the system Y=a(1)a(2), 
Eq. (1) for the energy becomes 


w= { weHvde 


27 1 
— {2z*— 244 —— 2e-(1 +-)]| Wu, 
8 p 


where p=ZR/ap. This interaction is repulsive 
at all distances and is shown by curve C, Fig. 1 
for the case that Z=27/16. When p equals 
infinity this expression becomes the familiar 
expression for the energy of a helium atom for 
hydrogen-like wave functions with variable effec- 
tive nuclear charge, W=(2Z?—27Z/4)Wzy. This 
energy is a minimum when Z=27/16 then 
W=-—5.695Wy. The experimental value and 
the value calculated by Hylleraas* for the helium 
atom is — 5.807 Wy. 


2. HYDROGEN ATOM AND HELIuM ION 
The symmetric variation function, 
W=a(1)b(2)+a(2)d(1), 


corresponds to an electron pair bond between 
the two atoms. The energy, is shown by curve B 


*E. A. Hylleraas, Zeits. f. Physik 65, 209 (1930). 


HYDRIDE 


MOLECULE-ION 





D 


1 
D ' 2 
R/a, 











Fic. 1. Potential energy curves. A, W=a(1)a’(2) 
+a’(1)a(2); B, W=a(1)b(2)+0(1)a(2); C, V=a(1)a(2); 
D, VW = a(1)a(2) + ki(a(1)b(2) + b(1)a(2)) + ke(a(1)a’(2) 
+a’(1)a(2)). 


in Fig. 1 for the case that Z=27/16=1.69. This 
energy can be minimized with respect to Z for 
each value of p. When this is done a value of 
W only 0.04 percent lower than that for Z = 27/16 
is obtained, the effective charge then being 
equal to 1.72. When p equals infinity the energy 
becomes W=(2Z?—6Z) Wy which is a minimum 
when Z=3/2 and W=—4.5Wy. The correct 
value for the sum of the energies of a hydrogen 
atom and a helium ion is —5Wy. This dis- 
crepancy is due to the fact that we have taken 
the same effective charge on the two atoms. 
This calculation is the same as that made by 
Glockler and Fuller except that here the same 
effective nuclear charge is used on both atoms 
and consequently all the integrals can be calcu- 
lated. This curve indicates the existence of a 
metastable attractive state for HeH*. 


3. Exc!iTED HELIUM ATOM AND PROTON 


The variation function 
W=a(1)a’(2)+a’(1)a(2), 


in conjunction with a(1)a(2), will be used to 
represent a polarized helium atom. The potential 
energy corresponding to the variation function 
a’(1)a(2)+<a(1)a’(2) alone is shown by curve A 
in Fig. 1 for the case that Z=27/16. This curve 
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is higher than that for the corresponding triplet 
state, a’(1)a(2) —a(1)a’(2), but the triplet state 
does not combine with the more important 
singlet states, nondiagonal matrix elements being 
equal to zero, so cannot be used. When the 
energy for the singlet state is minimized with 
respect to the effective nuclear charge, the 
energy curve is lowered by about 0.35Wy, the 
corresponding value of the effective nuclear 
charge at the minimum being 1.26. 


4, COMBINATION OF §§1 AND 3 


Using VW =a(1)a(2)+(a’(1)a(2) +a(1)a’(2)) as 
a variation function, the energy appears as a 
root of a quadratic secular equation, Eq. (1). 
The result of solving the equation is shown by 
curve B in Fig. 2 for the case that the effective 
nuclear charge is 27/16. It is not possible to 
minimize easily with respect to the effective 
nuclear charge so the quadratic is merely solved 
for several values of the effective charge. It is 
found that the energy is worse for any value of 
the effective charge much different than 27/16. 
The minimum energy for this curve is 
W=-—5.750 at R=1.60a. The normalized 

















-5.7 re 
B 
W/W 
Cc 
-S8r ain 
D 
1 ] | 
te) ' 2 3 4 
Rao 


Fic. 2. Potential energy curves. A, Y=a(1)a(2); B, 
W = a(1)a(2) + R(a(1)a’(2) + a’(1)a(2)); C, © = a(1)a(2) 
+ k(a(1)b(2) + b(1)a(2)); D, W = a(1)a(2) + ki(a(1)b(2) 
+a(2)b(1)) +h2(a(1)a’(2)+a’(1)a(2)). The dashed line 
represents the experimental value for the energy of the 
normal helium atom. 
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function is 


W =0.99a(1)a(2) +0.13(a(1)a(2) +a(1)a’(2)). 





The best value of the energy of the system at R 
equals infinity for this type of variation function 
is —5.695Wy, giving for the dissociation energy 
of the molecule 0.055Wy, or 0.74 ev. The stabili- 
zation resulting from the use of this variation 
function can be attributed to a polarization of 
the helium atom. 





5. COMBINATION OF §§1 AND 2 


The energy corresponding to the variation 
function Y=a(1)a(2)+k(a(1)b(2)+4(1)a(2)) is 
shown in curve C in Fig. 2 for the case 
that Z=27/16. The normalized function is 
0.97a(1)a(2) +0.24(a(1)b(2) + 6(1)a(2)). It is not 
easy to minimize with respect to the effective 
charge so again the quadratic secular equation 
is solved for several values of the effective 
charge, the calculated energy being much higher 
when the effective charge is much different than 
27/16. The lowest value obtained is W= —5.805 
at R=1.60ao, which, when subtracted from the 
value —5.695Wy for the helium atom, gives a 
dissociation energy of 0.110Wy or 1.49 ev. The 
form of the variation function suggests that this 
is due to the partial formation of a covalent bond. 

















6. COMBINATION OF §§1, 2 AND 3 


The best variation function used is of 
the form W=a(1)a(2)+2,(a’(1)a(2) +a(1)a’(2)) 
+ ko(a(1)b(2)+6(1)a(2)). On the basis of calcu- 
lations 5 and 6, in which it was found that the 
best value of the effective charge was not far 
from 27/16, calculations were made only for the 
case that Z = 27/16. The energy appears as a root 
of a cubic secular equation and is shown by curve 
D, Fig. 2 and by curve D, Fig. 1. The normalized 
function is W=0.97a(1)a(2)+0.10(a’(1)a(2) 
+a(1)a’(2))+0.20(a(1)b(2) +(1)a(2)). The low- 
est energy obtained is —5.844Wy at R=1.57a». 
This, when subtracted from the helium atom 
value of —5.695Wy, gives a value for the 
dissociation energy of 0.149Wy or 2.02 ev. This 
dissociation energy is obtained by using the 
energy of the helium atom calculated from a 
variation function of the same type as the one 
used for the undissociated molecule. This pro- 
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cedure assumes that the error in the energy is 
the same at the equilibrium distance and at 
infinity. A lower limit for the dissociation energy 
is obtained by subtracting the minimum energy 
for the molecule from the experimental value, 
—5.807Wuy, for the helium atom. The result is 
0.037Wy or 0.50 ev. This positive result proves 
conclusively that the molecule is stable. The 
value 2.02 ev, however, is probably much closer 
to the actual value of the dissociation energy. 

Fitting a parabola to the bottom part of the 
curve leads to a value for the force constant of 
(0.43 megadyne/cm and a fundamental vibration 
frequency of 2800 cm". 


THE MoLecuLe-Ion HeHtt 


Doubly charged helium hydride molecule-ion 
was treated by the same type of variation 
function, W=a(1)+,b(1)+.a'(1). The result- 
ing potential energy curve shows no minimum 
and differs only slightly from —4Wy+e?/R. The 
one electron bond between the two atoms is not 
sufficiently strong to overcome the Coulomb 
repulsion. Another molecule in which there is 
Coulomb repulsion between the atoms, the 
helium molecule ion, He:**, has been treated by 
Pauling’ including ionic terms in the variation 


HYDRIDE 


MOLECULE-ION 


function. In this case the electron pair bond 
formed is strong enough to overcome the 
Coulomb repulsion of the two atoms, there being 
a minimum in the potential energy curve. The 
instability of HeH** is due to the lack of 
degeneracy of the two states He+H* and He**H. 


DISCUSSION 


From a comparison of curves B and C, Fig. 2, 
it is seen that the introduction of polarization 
into the variation function is about half as 
effective in stabilizing the molecule as the intro- 
duction of the covalent bond. It is therefore 
qualitatively correct to say that two-thirds of 
the stabilization is due to the formation of a 
covalent bond and one-third is due to polarization 
of the helium atom. The strength of the bond, 
2.02 ev, is surprisingly great in view of the small 
coefficients of the covalent bond term and the 
polarization term, 0.20 and 0.10, in the variation 
function given in section 6. It is roughly half as 
strong as the bond in the hydrogen molecule, 
4.10 ev, calculated® from the same type of 
variation function. 

I wish to thank Professor Linus Pauling for 
making many valuable suggestions during the 
preparation of this paper. 


APPENDIX 


The simpler integrals are omitted. The integrals, Moo:, oand Moo, 0, are evaluated by the use of the Neumann expansion 
in confocal elliptic coordinates. The functions H(m, n, p) and S(m, n, p) are defined by Rosen.® 


M woo, meee dv dv. = 


» a?(1)a(2)b(2) 2{ ( 
ciuininenamatinineniecoenidiniueaietio Pie e p 
ri2 do 
»a?(1)a’(2)b(2) 


Moo , 0 = 
ri2 2 dy 480 


) 


gs 
dv\dv2.= £ Movo, 2+ — £. [A {| —225H(5, 4, 0) +135H(5, 2, p)+195H(3, 4, p) —117H(3, 2, p) 


+225.5(3, 5, p) —60.S(1, 5, p) +225.5(4, 4, ») —120.S(2, 4, ») —195.8(3, 3, p) +72S(1, 3, p) —135.S(4, 2, p) 
+725S(2, 2, p)} +B{375H(5, 4, p) —225H(5, 2, p) —195H(1, 4, p) +117H(1, 2, p) —375.8(3, 5, p) +100S(1, 5, p) 
—375S(4, 4, p) —125.8(2, 4, p) +120S(0, 4, p) +225.8(4, 2, p)+75.5(2, 2, p) —72.S(0, 2, p) +195.S(3, 1, p) 

—72S(1, 1, p)}+C{ —375H(3, 4, ») +225H(1, 4, ») +225H(3, 2, ») —135H(1, 2, p) +375.5(2, 4, p) —100.S(0, 4, p) 
+3755S(3, 3, p) —225.S(3, 1, p) —100S(1, 3, p) —225S(2, 2, p) +60S(0, 2, p)+60S(1, 1, p)}] 





»a(1)a’(1)a(2)b(2 Z p® 
Maw, =f a(1)a"(1)a(2)0( dias 5 Mow, += — [A {9H(5, 2, p) —3H(5, 0, p) +9H(3, 2, p) —3H(3, 0, p) —9.S(4, 2, p) 
ao 


Tie 


+35S(4, 0, p) —9S(1, 5, p)+3.5(2, 0, p)}+B{—27H(5, 2, 9) +9H(5, 0, p) —9H(1, 2, p)+3H(1, 0, p) 
+27S(4, 2, p) —95S(4, 0, p)+9S(2, 2, p) —3.S(2, 0, p)+27S(1, 5, ) —3.S(1, 1, p)} +C{27H(3, 2, p) —9H(3, 0, p) 
—9H(1, 2, p) +3H(1, 0, p) —27S(2, 2, p)+9S(2, 0, p)—27S(1, 3, p)+9S(1, a p)} 1, 


where 


+1 +1 atl 
A= J xe "dx, B= ft, x'e’*dx, and c=. xe? dx, 
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Dissociation of Excited Iodine Molecules 


E. RABINOWITCH AND W. C. Woon, Sir William Ramsay Laboratories of Inorganic and Physical Chemistry, 
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(Received March 2, 1936) 


By means of a sensitive optical method described .else- 
where, we determined the quantum yield of the process 


I,+hv—>I+I1 


in mixtures of I, vapor with several hundreds of mms of 
He, A, He and Nao, in three spectral regions: (1) In the 
continuum, i.e., on the short wave-length side of the con- 
vergency limit at 4990A; (2) in the band region imme- 
diately behind this limit, on its long wave-length side; 
and (3) in the band region further away from the limit. 
If the quantum yield is ascumed to be yi:=1.00 in the 


first region, the values y2=1.14 and y3=1.10 are ob- 
tained for the two other regions, showing that all excited 
I, molecules dissociate by collisions with foreign mole- 
cules. The quenching of I, fluorescence by foreign gases 
is thus practically never due to the dissipation process 
(1) I,*4+X—I,+X but entirely to the dissociation process 
(2) Ip*+X—-I1+1+X. The process (2) apparently occurs 
in most cases at the very first kinetic collision of excited 
iodine with a foreign molecule; only in the case of helium 
ten or more collisions are necessary to produce dissociation 
according to (2). 





1. INTRODUCTION 


HEN a molecule of iodine vapor absorbs 

light of wave-length greater than 4990A, 

an excited Iz molecule is produced, of which the 

natural mean lifetime is of the order of 10-7 sec. 

In sufficiently diluted vapor, the return into the 

normal state takes place under emission of fluo- 

rescence. The fluorescence is strongly quenched 

by an increase in the iodine pressure or by addi- 

tion of a foreign gas. Two mechanisms of quench- 
ing are possible: 


Ip*+X—1.4+-X, (1) 
1*+X14+14+X. (2) 


In the first case, the activation energy is dissi- 
pated in the form of kinetic, vibrational and 
rotational energy of the collision partners; in the 
second case it is mainly transformed into poten- 
tial chemical energy of free iodine atoms. Quench- 
ing by dissociation is of course only possible if the 
activation energy is greater than the dissociation 
energy of the iodine molecule. This is, however, 
the case in the whole of the visible spectrum and 
even in the near infrared, up to 8200A. Some 
indirect arguments speak for the predominant 
part played by reaction (2): 


(1) The high efficiency of quenching, as compared 
for instance with that observed in the resonance 
fluorescence of mercury vapor. This latter re- 
mains practically unaffected by nitrogen, helium 
or argon. up to a pressure of several hundred 


millimeters. The intensity of iodine fluorescence 
is reduced to one-half its original value already 
by a few millimeters of nitrogen or argon and by 
something like 15-20 mm of helium. The elec- 
tronic transition is of the same kind in both gases 
(?II—' in Iz, *P->'S in Hg). The probability of 
quenching by dissipation can therefore be sup- 
posed to be of the same order in both cases. It is 
reasonable tg ascribe the strong quenching of I; 
fluorescence to the additional possibility of 
quenching by dissociation, which is impossible in 
the case of mercury. (Mercury fluorescence is 
strongly affected by hydrogen, probably due to 
the process Hg*+H:—-Hg+H+H, which is 
analogous to (2).) 

(2) Magnetic quenching. I, fluorescence can be 
quenched by a magnetic field. In this case, no 
process of the kind (1) is possible, because of the 
lack of an acceptor for the activation energy. A 
process analogous to (2) is, however, possible, 
because the magnetic field can be substituted for 
the particle X in a ‘‘catalytical’’ dissociation of 
the excited I, molecule. This is in fact the expla- 
nation which was given to the magnetic effect by 
Turner! and Van Vleck.’ 

(3) From observations of the quantum yield of 
the bromine-hydrogen reaction Jost? concluded 
that the mechanism of this reaction is probably 
the same in the band region as in the continuous 
region of the bromine spectrum. It is therefore 

1L. A. Turner, Zeits. f. Physik 65, 464 (1930). 

2J.H. Van Vleck, Phys. Rev. 40, 544 (1932). 


3 W. Jost, Zeits. f. physik. Chemie 134, 92 (1928); (B) 3, 
95 (1929). 
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probable that a great number of primarily excited 
Brz molecules dissociate into atoms by collisions 
with He and Bre2 molecuies. All these results 
together form a strong indirect evidence for the 
high probability of the reaction (2) and its pre- 
ponderant role in the quenching of iodine fluo- 
rescence. The following experiments give a direct 
and quantitative proof of this assumption. 


2. EXPERIMENTAL 


In several previous papers’: ° we described a 
method of measuring the concentration of free 
atoms in strongly illuminated halogen vapors (or 
solutions). The method consists in observing the 
reversible decrease in the molecular extinction of 
the vapor which occurs when the absorption cell 
is exposed to intense illumination. The scheme of 
the apparatus is shown in Fig. 1. In discussing the 
results obtained with I, solutions in carbon tetra- 
chloride and hexane® we mentioned that the 
quantum yield of dissociation of Iz appears to be 
the same on both sides of the convergency limit 
at 4990A, thus indicating a dissociation of 
primarily excited I, molecules by collisions with 
the molecules of the solvent. The measurements 
in solution were, however, very rough indeed, 
because of the extreme smallness of the dissocia- 
tion effect in a medium in which the recombina- 
tion can occur at each collision of two iodine 
atoms. 

















Fic. 1. LZ, parallel light beam from a 1500-watt auto- 
matic carbon arc; A, quartz vacuum cell (12 X2X2 cm); 
TP, thermopile; FL, 6-v filament lamp; PC:, PC2, selenium- 
iron photo-cells, for measuring the light intensity J, and 
(after compensating the currents i and 7’) the intensity 
change AI occurring when A is illuminated by JL; g, 
Zernicke galvanometer (Type Zc, sensitivity 1 <10-! amp. 
per mm). 


— 


‘E. Rabinowitch and H. L. Lehmann, Trans. Faraday 
Soc. 31, 689 (1935). 

°E. Rabinowitch and W. C. Wood, Trans. Faraday Soc. 
32 (1936). 
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We now repeated these experiments with iodine 
vapor. The dissociation effect is in this case large 
enough (of the order of 0.1 percent) to be exactly 
measurable by our optical arrangement, even 
when filtered light is used for illumination instead 
of the total light of the carbon arc. The most 
difficult task is not to measure the dissociation 
effect itself, but to determine with a sufficient 
accuracy the number of quanta absorbed per sec. 
in 1 cc of the vapor in each of the spectral 
regions used. 

(a) The spectral regions. Fig. 2 shows the energy 
distribution in the three spectral regions which 
were used for the experiments. No. 1 lies com- 
pletely in the continuous region of the I, spec- 
trum, No. 2 in the band region close to the con- 
vergency limit, No. 3 in the same region further 
away from the limit. The energy distribution 
curves were obtained by means of a Hilger quartz 
double monochromator with a thermopile fixed 
behind its exit slit. (In the case of insufficient 
light intensity a selenium-iron photo-cell was 
substituted for the thermopile. Its spectral sensi- 
tivity was calibrated by comparison with the 
thermopile.) A certain difficulty is caused by the 
CN bands at 3860 and 4160A, the intensity of 
which varies even when the arc itself is burning 
quite steadily. An aesculin filter was used to de- 
crease, as far as possible, the intensity especially 
of the 3860-band. 

(b) The mean extinction coefficients a were 
determined in two ways: 


(A) By calculation, using the energy distri- 
bution curves given in Fig. 2, and the extinction 
curve of iodine vapor communicated by us in a 
previous paper.® The dotted curves in Fig. 2 show 
distribution of the light intensity absorbed by 
iodine in the three regions. The mean extinction 
coefficients are 


a= fal dr/ fI,dyx. (3) 
They were calculated by integrating graphically 
the curves given in Fig. 2. 

(B) The values so obtained were checked ex- 
perimentally, with a 6-v filament lamp as light 
source. It was impossible to use the carbon arc 
itself for exact measurements of a, because of its 
insufficient steadiness. In the regions No. 2 and 
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Fic. 2. Spectral distributions. Region No. 1: Filters Wratten No. 2 (aesculin)+No. 32+1 cm saturated CuSO, solu- 
tion. Region No. 2: Jena glass GG11 (2 mm)+1 cm Cu(NOs)2 solution (2 vols. saturated solution, 1 vol. HO). Region 
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No. 3: Wratten filter No. 21+11 mm CuSO, (1 vol. saturated sol., 1 vol. H»O). 


No. 3, the substitution of the filament lamp for 
the arc was permissible, as confirmed by the good 
agreement between calculated and experimental 
values of a in Table I. In region No. 1, however, 
the energy distribution in the spectrum of the 
filament lamp is rather different from that in the 
arc spectrum, partially because of the difference 
in temperature, but more especially because of 
the absence of the CN bands in the lamp. We 
attempted therefore to measure @ for this par- 
ticular region with the arc as light source as well. 
Single measurements gave strongly divergent 
results, because of the variations of the light 
intensity during the freezing out of the iodine 
vapor in the absorption cell. By averaging 40 
single measurements, however, an a-value was 
obtained which was in good agreement with the 
calculated value (see Table I). Both values were 
10 percent lower than the one obtained in ex- 
periments with the filament lamp. This difference 


TABLE I. Mean extinction coefficients. 
[a =(760/pd) log (Io/Z) ] 








SPECTRAL REGIONS 
1 No. 2 No. 3 


30.6 15.9 


14.7** 
i3.0°" 


No. 
8.8 


— 94" 
7s” 





1. Calculated 
2. Observed with 6-v 
lamp (16°C) 
(0°C) 
3. Observed with 25-amp. 
carbon arc 8.9* 


29.8** 
oF a 








* Measured in vacuum. 
** Measured in a helium atmosphere, p =500 mm. 


can be accounted for by the CN bands, which 
increase the intensity of the arc in the region be- 
low 4300A, where the absorption by iodine 
becomes negligibly small. 

In determining the ‘“‘calculated”’ values of a, 
the “‘limiting’’ extinction curve’ was used in 
the band region. This curve is obtained by 
using very ldw optical densities or by adding a 
sufficient amount of a foreign gas, in order to 
broaden the absorption lines of iodine so as to 
avoid an excessive absorption in the centers of 
the lines. The proof that these limiting conditions 
were realized in our experiments is given by the 
two sets of “experimental” values in Table I, 
obtained in a helium atmosphere (6=500 mm) 
with iodine vapor saturated at 16° (p(I2) =0.14 
mm) and 0° (p(I:)=0.030 mm), respectively. 
Their identity shows the validity of Beer’s law, 
which, in regions of band absorption can hold 
only under the ‘‘limiting’’ conditions. 

In the following, the ‘‘calculated”’ values of a 
(8.8; 30.6 and 15.9) were used as the most 
reliable ones. 

(c) The number of quanta absorbed was de- 
termined by calculating the mean wave-length 
of the light absorbed by iodine in each of the 
three regions. Graphical integration was again 
used. 


A= Sraylydd/ fapTydn. 


The values obtained in this way 


(4) 


4=4720A 5180A 5660A 
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were only slightly different from the wave- 
lengths corresponding to the maxima of the three 
absorption curves in Fig. 2 (4750, 5150 and 
5580A). The total absorption of the illuminating 
light in the cell was only 0.8 percent in region 
No. 1, 2.8 percent in region No. 2, and 1.5 percent 
in region No. 3. The number of quanta Nj, 
absorbed in 1 cc per sec. was therefore simply 


(5) 


where L means the intensity of the illuminating 
light (which was measured with a thermopile- 
galvanometer arrangement). The value of the 
constant in (5) is irrelevant, because we are con- 
cerned only with the comparison of the dissocia- 
tion effects caused by the same number of ab- 
sorbed quanta in the three spectral regions. 


Ni,=const. Lad, 


3. CALCULATION OF RESULTS 


The change A/J in the transmitted light inten- 
sity J (Fig. 1) caused by illumination L is propor- 
tional to the stationary decrease in concentration 
of I. molecules, i.e., to the concentration of I 
atoms in the illuminated gas. In experiments with 
Br. vapor* we made it plausible, that the re- 
combination of free halogen atoms occurs in 
agreement with the kinetic rules derived theo- 
retically for the elementary process 


Br+Br+X-—Br.4+ X. (6) 


We have been able since to prove this assumption 
by exact measurements with iodine, which results 
shall be published in a forthcoming paper. In 
order to obtain a practically pure homogeneous 
equilibrium: 


I.+hvI+I, 
I+I+X-I1.4+-X, 


(7a) 
(7b) 


it is necessary to use pressures of X sufficiently 
high to prevent the walls from playing an ap- 
preciable role in the recombination.* We found 
that these conditions were realized under our 
geometrical conditions (and the light intensities 


—. 


* It is of course possible to work in the region of mixed 
homogeneous-heterogeneous recombination as well, but in 
this case the light intensities must be adjusted so as to 
give exactly the same numbers of quanta absorbed in all 
the three spectral regions, because under these conditions 
4 is an unknown function of L (between const. XZ and 
const. x L4), 
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used), at pressures above 150 mm with all gases 
except He and Hg, where pressures of the order of 
400-500 mm are to be used. In the region of 
homogeneous equilibrium (7), the stationary 
concentration of iodine atoms is proportional to 
the square root of their rate of production. This 
rate is: 

+d[I]/dt=2yNi,, (8) 
where y is the quantum yield of dissociation. 
Therefore: 


27 Ni,=const. (AJ)’, (9) 


or on account of (5) 


y=const. (AJ)?/Lan. (10) 


Our purpose is to compare the rate of production 
of atoms in region No. 2 (and No. 3) to that in 
region No. 1, where the quantum yield is surely 
yi=1. According to (10) the quantum yield is: 
yo= ((AT)o?/(AD);?)(Liasd1/Lea@ede2). (11) 

A corresponding equation holds for y3. 
In Table II, results are given of experiments 


TABLE II. Dissociation effect in different spectral regions 
p(I2) =0.14 mm. 








No. 3 
A =5660 | 
@=15.9 


No. 2 | 
\ =5180 


L Lad Al 

5.40 4.86 13.9] 1. 
4.32 8.9/1. 
40 99 
9.8 10.8) 
11.8 17.4/1. 














with helium, hydrogen, argon and nitrogen as 
foreign gases, added to 0.14 mm of iodine vapor. 
The last two columns show the values of y2 and 
73 obtained in different gases. 


The magnitude of the dissociation effect observed in 
these experiments (which we mention for orientation, 
although it is irrelevant for the calculation of y) is given by: 

AT Alls], ? ain 10. 

I {iIe] 760 
The values involved were @ (extinction coefficient for the 
photo-cell light J)=20; p (pressure of iodine) =0.14; d 
(length of the cell)=12 cm. The relative changes of the 
light intensity AJ/I were of the order of 2X10~. This 
means 

ATI. ] 2x10 760 


= =0.002. 
[Iz] 20x0.14x12X2.3 





About 0.2 percent of the I, molecules were thus dissociated 
into atoms in the stationary state during illumination. 
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The values of y:, and ys shown in Table II, are 
all equal to 1 within the limits of the experimental 
error. We may thus conclude, that beginning from 
the convergency limit at 4990A and up to the red 
limit of region No. 3, 1.e., up to 6200A, practically 
very excited I, molecule dissociates into atoms by 
collisions with molecules of helium, argon, nitrogen 
or hydrogen. 

The quenching of iodine fluorescence by all 
these gases is thus entirely due to the dissociation 
reaction (2). The velocity of (2) is therefore equal 
to that of the quenching process. The quenching 
curves, measured by Franck and Wood’ and 
more recently by Berg and Réssler* show that in 
case of all gases except helium the quenching 
occurs at the very first (or one of the first few) 
kinetic collisions. This is shown by the following 
calculation: The quenching constants for A, Ne 
and Oz are, according to Berg, of the order of 
1.5-2. These are the constants of the Stern- 
Volmer equation J,=J,/(1+ const. X p.), J) meas- 
uring the intensity of fluorescence without 
quenching, J, that in presence of », mm of the 
gas X. (The value of J) must be extrapolated to 
p(12) =0 in order to account for the self-quench- 
ing by iodine.) These values of the constants 
mean that J,=J)/2 at p,=—ca. 0.7 mm=7 X10 
atmos. The fluorescence is reduced to 1/2 its 
original value when the mean time between exci- 
tation and emission is equal to the mean time 
between excitation and quenching. The first 
being about 10-7 sec. ® in the case of Is (IIo), the 
second must be of the same order as well. The 
mean kinetic collision-interval is 10~ sec. at 1 
atmos.; and 10-7 at 0.001 atmos., and is thus 
identical with the mean ‘“‘quenching time.’’ The 
I,* molecules are thus unable to stand gas-kinetic 
collisions with A, Ne or Oz» without dissociation. 
They are somewhat stabler against collisions with 
helium molecules. A He pressure of 20 mm is 
needed for reducing the intensity of fluorescence 


7 J. Franck and R. W. Wood, Verh. d. D. phys. Ges. 13, 
78 (1911). 

8 W. Berg, Zeits. f. Physik 79, 89 (1933); F. Réssler, Zeits. 
f. Physik 96, 251 (1935). 

9 We think that the lifetime of J2(*IIo) must be nearer 
10-7 than 10-8 sec. because of the analogy with the °P, 
state of Hg, whose lifetime is well known to be 1077 sec. 
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Fic. 3. Potential curves of the iodine molecule. 


by 50 percent in green and something of the order 
of 70 mm to reduce it by the same amount in the 
red. In our experiments at 500 mm helium no 
trace of fluorescence was recognizable and the 
dissociation even in the long wave region No. 3 
was as complete in helium as in all other gases. 

The mechanism of the secondary dissociation 
of excited I,* molecules is probably the one indi- 
cated by Turner,’ and theoretically elucidated by 
Van Vleck.’ It is illustrated by Fig. 3, taken from 
Van Vleck’s paper. The transition from the 
“stable”’ excited state *IIp+ to the unstable repul- 
sive state *IIp~ is prohibited in a free iodine 
molecule; the symmetries which make it im- 
possible are, however, destroyed by a magnetic 
or an inhomogeneous electric field. The activated 
molecules I,* dissociate therefore under the influ- 
ence of a magnetic field or of the inhomogeneous 
molecular fields which act on them during the 
collisions. The transition becomes more difficult 
when the potential curves *IIp+ and *IIp- are more 
widely separated,—this is why the fluorescence is 
less sensitive in the red than in the green region. 

At the end, we wish to express our heartiest 
thanks to Professor F. G. Donnan for giving us 
the opportunity of working in the Sir William 
Ramsay Laboratories and for his interest and 
invaluable support during the work. 
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Free Energy and the Rate of Chemical Reactions 


The Relation Between the Rate and Oxidation-Reduction Potentials 


HAROLD GERSHINOWITZ, Department of Chemistry, Columbia University 
(Received March 12, 1936) 


The theories of the rates of chemical reactions previously developed by Eyring, Rice and the 
author give a relationship between the rate of reaction and the free energy of formation of the 
activated state. With the aid of this relation it is possible to explain the experimental data which 
have shown that there is a connection between the rates of oxidation-reduction reactions and 


the potentials of the reagents used. 


I. INTRODUCTION 


OME years ago, Conant and his co-workers! 

pointed out that there was a definite relation- 
ship between the rate of an irreversible oxidation- 
reduction reaction and the oxidation-reduction 
potentials of the reagents that are used. With the 
aid of some assumptions concerning the mecha- 
nism of such reactions, they were able to derive 
an equation relating the relative rates of reaction 
with two different reagents with the normal 
potentials of these reagents. It was found that 
the agreement with experimental data was satis- 
factory. The studies of Conant and his collabo- 
rators have been continued and extended by 
Chow,? who also has found that a relationship 
between oxidation reduction potentials and the 
rates of some reactions can be expressed by a 
simple empirical equation. This work has made 
it evident that there exists some connection 
between the kinetic rate constant and thermo- 
dynamic quantities. Such a conclusion is also 
indicated by the data in the fields of acid and 
basic catalysis. 

Until quite recently, there seemed to be no 
theoretical justification of the existence of such 
relationships. The collision theory of the rates of 
chemical reactions had no connection with 
thermodynamics. Within the past year, however, 
Eyring, Rice and the author* have developed 
theories of the absolute rates of chemical re- 
actions which serve to relate the data of kinetics 

‘Conant, Chem. Rev. 3, 1 (1926); Conant and Pratt, 
J. Am. Chem. Soc. 48, 3178, 3220 (1926), and earlier 
papers. 

* Chow and Kammerling, J. Biol. Chem. 104, 69 (1934); 


Chow, J. Am. Chem. Soc. 56, 894 (1934); ibid. 57, 1437, 
1440 (1935), 

*(a) Eyring, J. Chem. Phys. 3, 307 (1935); (b) O. K. 
Rice and H. Gershinowitz, ibid. 2, 857 (1934) and (c) ibid. 
3. 479 (1935). 


with thermodynamic quantities. O. K. Rice and 
the author have shown how one may determine 
the absolute rates of unimolecular decomposi- 
tions and bimolecular associations from the 
entropy change involved in the _ reactions. 
Wynne-Jones and Eyring‘ have related the ab- 
normal rates of some reactions in solution to the 
entropy of formation of an ‘“‘activated complex.”’ 
The rate can also be simply related to the free 
energy change involved in the formation of an 
activated state. Since the free energy change and 
the e.m.f. are directly proportional, it is possible 
to obtain a theoretical derivation of a relation 
between the rate of a reaction and the oxidation- 
reduction potentials of the reagents involved. 
In the following, we shall borrow freely from 
the notation and terminology of Eyring and his 
co-workers as well as those of Rice and the 
author, in accordance with the method used in 
reference 3c. 


II. THE RELATION BETWEEN THE SPECIFIC RATE 
CONSTANT AND OXIDATION-REDUCTION 
POTENTIALS 


It has been shown that the rate of a chemical 
reaction is given by the concentration of the 
reacting systems in a certain fraction of the total 
phase space that is available to the systems, 
multiplied by the velocity with which the 
systems are passing through this region.* The 
configuration of the reacting molecules which 
corresponds to this region of phase space may 
be called the activated complex. That is, if one 
wishes to describe the particles in the systems by 


4 (a) Wynne-Jones and Eyring, J. Chem. Phys. 3, 492 
(1935); see also (b) Evans and Polanyi, Trans. Faraday 
Soc. 21, 875 (1935). 
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means of their coordinates and momenta, the 
concept of position in phase space is used. If one 
prefers to describe the particles by means of 
their coordinates and potential energies, one 
speaks of the activated complex. It has been 
shown that the two methods give the same 
results. The velocity of these systems or activated 
complexes, expressed in units of cells of phase 
space of size h per second, is kT/h, where T is 
the absolute temperature and k and ih are the 
Boltzmann and Planck constants, respectively. 
If (C) is the concentration of activated complexes 
per unit length h, the specific rate of any re- 
action is given by 


ky =(C)RT/h=(C)-2, (1) 


where we write v for kT /h in order to simplify 
the writing of subsequent equations. 

Now if we represent the equilibrium reaction 
by which the activated complex is formed as 


nA+mB2C (2) 
the equilibrium constant for the reaction is 
K=(C)/(A)"(B)”. (3) 


At unit concentration of reactants K=(C) and 
the specific rate constant becomes‘ 


ky =K-v. (4) 
From thermodynamics we get the relation 
In K= —AF*/RT, (5) 


where AF* is the free energy of formation of the 
activated state. Substituting for K in (4),5 


ki=exp (—AF*/RT)-v. (6) 


In general, AF* is not the same as AF, the free 
energy change for the total reaction. The calcu- 
lations of the absolute rates of unimolecular 
decompositions and bimolecular associations in 
the gas phase*>:-* indicated that the entropy 
change for the formation of the activated state 
is, in most cases, equal to zero for the decom- 
positions and equal to the entropy change for 
the total reaction in the case of the associations, 
i.e., the activated state has approximately the 
same entropy as the normal associated molecule. 
The calculations of Wynne-Jones and Eyring‘ 


5A similar relation between the free energy and the 
reaction rate has been obtained by V. K. La Mer (J. Chem. 
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show that for many reactions in solution, also, 
the entropy change for the formation of the 
activated complex is equal to the total entropy 
change in the reaction. For the free energy 
change in such reactions to be equal to the free 
energy of formation of the activated complex, it 
would be necessary that the energy of formation 
of the activated state be equal to the energy of 
the reaction. This will be so only if the energy 
of activation of the reverse reaction is zero, 
since the energy change for the reaction is equal 
to the difference of the activation energies of the 
forward and reverse reactions.* We see, there- 
fore, that the free energy of formation of the 
activated state will rarely be equal to the free 
energy change for the total reaction.’ It is this 
fact that has obscured the relation between free 
energy and the rate of reaction. 

For a reversible reaction, we have the relation 


AF* = —nFE*, (7) 


where 7 is the number of equivalents, F the 
Faraday and E* the e.m.f. for the reaction by 
which the activated complex is formed. Sub- 
stituting (7) in (6) we have, finally, 


,ki=exp (nFE*/RT)-v. (8) 


With Eq. (8) as a basis, we can now determine 
the relation between the potential of a reagent 
and the rate of a reaction. For organic oxidations, 
Conant and Pratt proposed the following mecha- 
nism 

ROH+B@RO-+ BH 
RO--—R=0. (9) 
| 


This type of reaction can readily be treated by 
the theoretical method that we have just pro- 
posed. An added advantage of such a treatment 
is that, in order to obtain relative rates, it is not 
necessary to assume as particular a mechanism 


Phys. 1, 291 (1933), Eq. (10) et seq.), by means of thermo- 
dynamic analogies from Tolman’s general expression {or 
the temperature coefficient of the rate of reaction. The 
most important difference lies in the fact that in the 
expression we have derived there is no arbitrary constant, 
since we have obtained directly the integrated form of the 
expression for the specific rate constant so that there is no 
question of an undetermined constant of integration. 

* Gershinowitz and Rice, J. Chem. Phys. 2, 273 (1934). 

7 There does seem to be at least one reaction for which 
the free energy change that determines the rate of the 
reaction is the same as the total change in free energy. 
This is the autoxidation of hydroquinone, catalyzed by 
manganous salts in acid solutions, studied by La Mer and 
Temple (Proc. Nat. Acad. Sci. 15, 191 (1929)). 
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as that which Conant and Pratt used to explain 
the relation between e.m.f. and the rate of 
reaction. If one does however, postulate a 
definite mechanism, it becomes possible to calcu- 
late not only the relative rates, but also the 
absolute values of the rate constants. It is also 
possible to extend the theoretical equations to 
include the more complex systems studied by 
Chow. 

Let us assume that we have a substance A 
which is to be oxidized by some reagent B, and 
that the activated complex C consists of some 
oxidized form of A, say A’, and B’, the reduced 
form of B. At this point it is not necessary to 
know the further mechanism of the reaction, or 
whether the oxidized form of the molecules A is a 
single molecule or a polymolecular complex. All 
that we need know is that the rate determining 
step for the reaction is the passing of this complex 
through a certain region of phase space, no 
matter whether the motion represents a uni- 
molecular rearrangement or some bimolecular or 
polymolecular change. The change in potential 
involved in the formation of the activated 
state is 


E*=Ea:—Eat+Ep —Exp. (10a) 


At unit concentration of reactants 
Ep: —Ep=Eos, (10b) 


where Eog is the normal oxidation-reduction 
potential of B. Hence, from Eq. (7), 


ki=exp (nF(Ea4,-—Eat+Eor)/RT)-v. (11) 


If instead of B, we used some other oxidizing 
agent D, the rate would be given by 


ky’ =exp (mF (E4,-—Ea+Eop)/RT)-v. 


If n=m, 
ki/ky'=exp [(Eon—Eop)nF/RT] 
or In (ki/k1’') =_Eon—Eop |-nF/RT. 


(12) 


(13a) 
(13b) 


By substituting numerical values for the uni- 
versal constants, and letting »=2, as in some 
cases studied by Conant and Pratt, 


0.03 log (k:/k1’) =Eor—Eop. (14) 


This equation is identical with Eq. (11) of 
Conant and Pratt, which was found to give 
satisfactory agreement with the experimental 
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data. We see that this relationship is a direct 
consequence of the new theories of the absolute 
rate of chemical reactions, and that it may be 
derived without assuming the unimolecular 
mechanism that Conant e/ al. found it necessary 
to use. 

An essential distinction between the method 
we have used and that employed by Conant and 
his co-workers is involved in the nature of the 
equilibrium processes. Conant and Pratt,® for 
example, assume the following mechanism for 
the reduction of an azo dye, 


RN = NR’OH+ BH2@RNHNHR’OH+B 
O 
VA 
RNHNHR’/OH—RNH2+ R’ 
slow XQ 
NH 


The second of these reactions is the one that may 
be considered to go through an activated com- 
plex. But according to the method that we have 
used, it does not matter how many intermediate 
steps there are between this step and the initial 
reactants. As long as the activated complex is 
the same, the rate of the reaction will depend 
only on the potential of the reducing agent. 
Rice® has made clear the distinction between flux 
equilibrium and “‘back and forth’ equilibrium. 
In the theory of the activated complex the 
complex is in flux equilibrium with the normal 
reactants. This is why an apparently equilibrium 
method can give information about rates of 
reactions that are not in equilibrium. It is the free 
energy of the state of the system that is in flux 
equilibrium that determines the rate of the 
reaction. All the intermediate configurations 
that are in “back and forth’ equilibrium will 
drop out of the calculation. It is thus clear that 
the method of the activated complex is more 
complete and rigorous than that used by Conant, 
which is a special case, and that it will explain 
all the data that the latter will. This explicit 
definition of the nature of the activated complex 
also distinguishes this treatment from that based 
on the Bronsted-Bjerrum type of complex, which 
is not rigorously defined. 

Chow has found that it is possible to explain 


8 Conant and Pratt, J. Am. Chem. Soc. 48, 2464 (1926). 
90. K. Rice, J. Phys. Chem. 4, 53 (1936). 
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the experimental data for the kinetics of various 
oxidation-reduction reactions by means of a 
simple empirical relationship which he has pro- 
posed, namely 


rate « exp (nFE/3RT). (15) 


Using this equation, he has obtained good agree- 
ment with the data on the catalytic absorption 
of oxygen by linseed oil, and the oxidation of 
formate and oxalate ions by halogens in the dark. 
While the general form of this relation is that 
which would be expected from Eq. (8), it seems 
rather difficult to construct a mechanism for the 
formation of the activated complex which would 
give the exact dependence on the potential that 
is given by Chow. It is possible to construct 
several mechanisms that give the relation, rate 
«exp (nFE/2RT), but while these give good 
constants within a run, they do not give good 
results from run to run. From the data it would 
seem that a dependence on (nFE/3RT) is 
actually demanded. Such a relation can be 
obtained from a mechanism which demands the 
formation of a complex from three equivalents 
of the reagent. It would seem that more data are 


needed to help decide what is the nature of 
the activated complex. A step in this direction 
is afforded by the recent work of Abel and 
Schmid,’ who have obtained indirect spec- 
troscopic evidence of the presence of the ion 
[CO.CO,.-] in the oxidation of oxalate ion by 
iodine. 


III. Discussion 


Bancroft and Magoffin"! have recently dis- 
cussed the role played by excited energy levels in 
electrochemical processes. The method which 
they have used may be correlated with that of 
Conant by means of the approach used above. 
They state that the concept of an intermediate 
energy level which they use is quite different 
from the assumptions of a reversible electrolytic 
process and a practically irreversible nonelectro- 
lytic one made by Conant and his collaborators. 
From the treatment we have given, however, it 
may be seen that both methods are essentially 
the same, being qualitative descriptions, in 


10 Abel and Schmid, Naturwiss. 23, 501 (1935). 
1 Bancroft and Magoffin, J. Am. Chem. Soc. 57, 2561 
(1935). 
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different words, of the fact that the rate of re- 
action is determined by the concentration of an 
intermediate configuration of the system, which 
is in equilibrium with the reactants. 

The development of rigorous thermodynamics 
and the kinetic theory explanation of the rates 
of chemical reactions served to discredit the idea, 
prevalent since the time of van’t Hoff, that the 
rate of a reaction should be determined by the 
free energy change. The theories of reactions 
recently developed, which take into account all 
the degrees of freedom of the reacting molecules, 
especially in the form developed by Eyring, 
which makes explicit use of an activated com- 
plex, show that the rate of a reaction is deter- 
mined by a free energy change, which is usually 
not the same as the free energy change for the 
total reaction. The theory of Bronsted and 
Bjerrum, which introduced an_ intermediate 
complex, explained some of the relationships 
between rates and quantities usually associated 
with equilibrium measurements. La Mer has 
shown the relation between the free energy 
change and the concept of the critical complex.’ 
In these treatments, however, there still re- 
mained arbitrary constants. Wynne-Jones and 
Eyring and Evans and Polanyit have shown 
quantitatively the relation between the theory of 
Bronsted and Bjerrum and the general theory 
of the activated complex which applies to gas 
phase reactions, as well as to reactions in 
solution. The regularities observed for oxidation- 
reduction reactions and for acid and _ basic 
catalysis should no longer be considered as 
isolated phenomena, but as being manifestations 
of the existence of some critical configuration of 
the reactants. This relation between the two 
types of reaction helps to explain why Abel and 
Hilferding" were able to describe formally the re- 
action between oxalate ion and iodide ion as 
either an oxidation-reduction reaction or a re- 
action catalyzed by H*. Also, as Bancroft and 
Magoffin have remarked, there is no need to 
make any fundamental distinction between 
hydrogenation and reduction. It should now be 
possible to treat all reactions, homogeneous Or 
heterogeneous, by the same theoretical methods. 


12 Abel and Hilferding, Zeits. f. physik. Chemie A172, 
353 (1935). 
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Liquid nitric oxide consists largely of double molecules, as evidenced by the abnormally high 
entropy of vaporization and the high and variable specific heat. By making use of these data it is 
possible to calculate the extent to which single molecules exist in the liquid, the heat of dissocia- 
tion of the double molecules, and the value the specific heat would have were the association 
in the liquid complete. A statistical calculation is made of the entropy of vaporization, and 
the result compared with the experimental value. An estimate is made of the heat of dis- 


sociation in the gas phase. 





§1. 


HE existence of associated nitric oxide 

molecules has been inferred from the equa- 
tion of state! and from the abnormally high heat 
and entropy of vaporization,” and Johnston and 
Weimer have remarked that the abnormally high 
and rapidly varying specific heat of the liquid 
indicates the beginning of dissociation of the 
double molecules near the boiling point.* But 
though Eucken and d’Or' did attempt to make 
some estimate of the degree of association in the 
gas and the dissociation energy of the double 
molecules, the full possibilities of the available 
data have not, so far as I am aware, been 
realized. 

Let us first consider the entropy of vaporiza- 
tion,‘ A,.S. As we have stated, this is considerably 
higher than it would be were nitric oxide a 
normal unassociated liquid. The entropy of 
vaporization of this hypothetical normal liquid, 
which may be estimated by means of Hildebrand’s 
rule or by a comparison with normal liquids 
which boil around the same temperature, we 
shall call AyS. This is the entropy of vaporiza- 
tion of the ideal normal liquid to give vapor at 
the vapor pressure of the actual liquid, e.g., at 
one atmosphere if the temperature is the ob- 
served boiling point. Let N be Avogadro's 
number, NV, the number of (NO). molecules, and 


' Eucken and d’Or, Nachrichten Gesell. Wiss. Géttingen, 
Math.-phys. Klasse (1932), p. 107. See also Kassel, 
Kinetics of Homogeneous Gas Reactions, p. 170. 

* See e.g., Johnston and Giauque, J. Am. Chem. Soc. 51, 
3194 (1929), especially pp. 3208, 3213. 

* Johnston and Weimer, J. Am. Chem. Soc. 56, 625 
(1934). See footnote 19. 

_‘ This quantity is referred to unit amount (one mole) of 
nitric oxide (NO) whether single or combined with another 
NO to form a double molecule. 


Nz the number of NO molecules in one mole of 
nitric oxide. We shall assume that NO and 
(NO): form a perfect solution, so that the vapor 
pressure Py of the unassociated liquid bears 
to the actual vapor pressure, P,, the ratio* 
(Ni+ N2)/Ne. This characterization of the perfect 
solution includes the assumption that the vapor 
is a perfect gas. It may appear that this is not 
always justified, for at the boiling point of nitric 
oxide Py may be expected to be several atmos- 
pheres. All difficulties are avoided, however, if by 
Py and P, we understand the hypothetical vapor 
pressures which would obtain if the vapor were a 
perfect gas, and if all heats of evaporation are 
taken to be the heats for the formation of such a 
perfect vapor. The only error which will then 
occur is that arising from the assumption, which 
we shall make when we use the observed heat of 
evaporation, that the real vapor is a perfect gas 
at one atmosphere and that it is completely 
dissociated. This error is negligible. 

Now the heat of evaporation of the actual 
liquid,‘ A,.77, is related to the heat of evaporation,‘ 
AnH, of the normal liquid, by the equation 


A, =AnH+(Ni/N)AH, (1) 


where A/J is the change of heat content on 
dissociation per mole of (NO). Hence 


A,S=AnH/T+(Ni/N)SH/T. (2) 


But Ay///T is the entropy of vaporization of the 
normal liquid to vapor at its equilibrium pres- 
sure, Py; to get AyS we must add the change of 
entropy when the pressure is reduced to P,. Thus, 
since Py/P,=(NitN2)/No2, 


* See Hildebrand, Solubility, 2d Ed., p. 19. (We assume 
the vapor pressure of the associated liquid to be zero.) 
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AnS=AyH/T+RI1n ((Ni+Ne)/Ne]. (3) 
From (2) and (3) 


R In [ No/(Ni+ Ne) | 
=A,S—AnS—(Ni/N)AH/T. (4) 


Now the change of P, with the temperature will 
be given approximately by the Clausius- 
Clapeyron equation Rd In P,/dT=A,11/T? and 
similarly Rd In Py/dT=Ay/1/T?, so that 


Rd \n (P,/Px)/dT 
= Rd In [N2/(Ni+N2)]/dT 
=A,H/T?—AyH/T?=(Ni/N)AU/T?. (5) 


Eq. (4) is an integrated form of Eq. (5). 

Let us now see what information may be 
gleaned from the specific heat data. If we assume 
that the excess specific heat,‘ C,’=C,—Cp, 0, 
where C, is the observed specific heat‘ and C,, 9 
is the specific heat of the hypothetical completely 
associated liquid,‘ we have 


C,’ =(dN2/dT)(AH/2N). (6) 
Now from (5) we have 


d[Ne2/(Ni+Ne2) ]/dT 
=[N2/(Nit Ne) ](Ni/N)AH/RT?, 


which yields, remembering 2N,;+N2.=N=con- 


stant, ; 
dN» 2(Ni+ Ne) Ni 5 AH 


aT N? RT? 





(7) 





Now since N,/N and (N,+Ne2)/N are both 
almost exactly equal tof 3 and since one of them 
is less than $ by the same amount that the other 
exceeds it we can safely setd N2/dT = N.AH/2RT?, 
which yields, assuming AH constant, 


N2=b exp (—AH/2RT), (8) 


where b is a constant. Since C,’ is proportional to 
dN2/dT and hence to N2/T? we see from (8) that 


C,’=aT~ exp (—AH/2RT), (9) 


where a is a constant. We may write (9) in the 
form 


C,’=a' exp [—(AH—4RT;)/2RT], (10) 


T, being the value of 7 at the boiling point. 
a’=aT~ exp (—27,/T), which may for purposes 
¢ Since the liquid is almost completely associated, ac- 


cording to assumption, which turns out to be consistent 
with the results subsequently obtained. 
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of approximation be regarded as constant. If we 
plot log C,’ against 1/7 we should thus get a 
straight line, from the slope of which A/7 may 
be calculated (still assuming it is independent of 
temperature). Experimentally, however, we have 
C,, not C,’. In order to obtain C,’, therefore, we 
need to know the value of C,, o. We have tried 
calculating C,’ by setting C,, » equal to arbitrary 
values, 10, 11, and 12 calories per degree per 
mole of nitric oxide. This neglects the fact that 
Cy», 9 may not be exactly constant, and also the 
fact that the specific heat of a mixture of double 
and single molecules (as exists after some of the 
latter have been formed) is slightly different from 
that of the pure associated liquid, even forgetting 
the heat absorption due to further dissociation. 
In all cases a fairly good straight line is obtained 
by plotting log C,’ against 1/7. In Table I we 
have listed the values of A/7 corresponding to the 
various assigned values of C,, o. 

It is now possible finally to determine AJ/ by 
integrating Eq. (6). We get 

Tb 


@/aH) f C,'dT=(N2/N)s, (11) 
0 


(The subscript 5 indicates that a quantity is to 
be evaluated at the boiling point.) From Eggs. 
(9) and (11) we have 


Tb 
(N2/N)o= (20/aH) f T-? exp (—AH/2RT)dT 
0 


=4Ra(AH)~ exp (—AH/2RT>) 
=4RT,°C’ », »(AH)-*. (12) 
We can thus determine the value of (N2/N), that 
corresponds to any given value of AH. This is 
also given in Table I. We can now substitute 
(N2/N), back into Eq. (4) and determine the 
value of AH given by that equation (after 
evaluating A,S—AyS as discussed in the next 
paragraph). The values of AH thus determined 
are likewise given in Table I. Now in order for 
the calculation to be self-consistent the two 
values of AH must coincide. This occurs if 
Cp, o= 11.1, (N2/N),=0.062, and AH = 3910. 


TABLE I. (Energy units, calories per mole.) 











Cp, 0 4H, Ea. (10) (N2/N)p AH, Ea. (4 
10 3390 0.092 3840 
11 3840 0.064 3900 
12 4440 





0.042 4010 
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The entropy values in Eq. (4) are obtained as 
follows. A,.S is known from the work of Johnston 
and Giauque? to be 27.13 entropy units per mole 
of nitric oxide. AyS is estimated on the assump- 
tion that the normal unassociated nitric oxide 
would be like oxygen, which boils at 90.1°K with 
an entropy of vaporization of 18.1 entropy 
units.° Assuming that the difference between the 
specific heats of the unassociated nitric oxide gas 
and liquid is the same as for oxygen, about 6 
calories per mole per degree, we calculate that at 
the boiling point of nitric oxide AyS=16.3 
entropy units. The chief source of error in our 
results is presumably the uncertainty in A,S, 
which might well be incorrect by one entropy 
unit, or possibly even somewhat more. An error 
of one entropy unit would cause an error of 300 
calories in the value of A//. 


§2. 

It will now be of interest to attempt to 
compare the value assigned to A,S—AyS with 
the value which may be obtained theoretically 
with the aid of statistical mechanics and reason- 
able assumptions as to the properties of the 
associated molecules. Let us consider the entropy 
of a mixture of NO and (NO). molecules con- 
tained in a volume V. If we assume that there 
are no forces between the molecules the total 
entropy’ may be written in the form® 


S=N1Si+ N2S2+hk In (N!/Ni!Ne2!), (13) 


where S,; and S: are quantities which may be 
called the entropies per molecule of (NO): and 
NO, respectively, and k& is the Boltzmann 
constant, R/N. The last term in (13) is the 
entropy contribution due to the possibility of 
having different individual NO’s combined in the 
(NO)» molecules. 


The entropies S; (i=1 or 2) may be written in the 
form 


Si=k In (VB,C,D;E;) = S;°+k In V, (14) 


*Giauque and Johnston, J. Am. Chem. Soc. 51, 2300 
(1929), especially p. 2314. 

‘This definition of entropy (differing by the thermo- 
dynamic entropy by a constant amount, which does not 
affect entropy differences) and the notation used in this 
part of the paper follow Rice and Gershinowitz, J. Chem. 
Phys. 2, 853 (1934). The general form of entropy expression 
(Eqs. (14), etc.) is that which was used by Rice and 
Gershinowitz and independently suggested and discussed 
by Halford, J. Chem. Phys. 2, 694 (1934) and further 
(ane by Ejidinoff and Aston, J. Chem. Phys. 3, 379 
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where VB;, C;, and /; are the contributions due, respec- 
tively, to the translation, the rotation, and the vibration 
of the molecule, and D; is that due to electronic de- 
generacy, etc. We have 


B,=(2rem,kT /h?)*”, (15) 


where e is the base of natural logarithms, m; is the mass of 
the molecule, and h# Planck’s constant ; 


Cp= Cp Cress, 
Cy =r (Brel RT /h?o 7)", 


Cor = (Biel RT [Ho ni, 


(16a) 
with 
(16b) 


refer to different types of rotation, 
involving, respectively, Ji, li, - degrees of freedom, 
and having moments of inertia J;, Ji, «++ (Tie >Ijn >+ ++) 
and symmetry numbers o;, oi", ; and £; is a compli- 
cated but well-known function of the quantities hv;;/kT, 
where the »;; are the frequencies of the molecule. 


where Cy, Cin, ++° 


Eq. (13) applies to a dilute gas; it may also be 
applied to a liquid, provided the forces are of 
such a nature that they can be taken into account 
by assigning to V a certain effective value, not 
equal to the actual volume of the liquid, but in 
general much smaller, since a molecule cannot 
occupy the space already occupied by another 
molecule. This assumption, together with the 
assumption that V is proportional to Ni+ Ns, is 
equivalent to the assumption that NO and (NO). 
form a perfect solution, as we shall see. 

We consider the entropy change in the process 
whereby the liquid changes from a completely 
unassociated one with N,=0, No=N, and V= Vy 
to the actual liquid with Ni=N,;, Ne= Ne, and 
V=V. By Eg. (13) this entropy change will be 
(using Stirling’s approximation for factorials) 
Ay)S= N,S;°+ NoSe° — NS. — Nk In Vu 


+(Ni+Ne)k In V+ NR In N—Nik In My 
— Nok In No+(NitN2—N)k, 


where the S;° are defined in Eq. (14). Remem- 
bering now the relation between V and Vy (the 
effective volume being proportional to Ni+ Ne) 
and that 2N,+N2=WN this reduces to 


ApS = N,AS°—(Ni/N)R In V 
+R In (Ni+Ne2)—(Ni/N)RIin MN, 
—(N2/N)R In Ne—(Ni/N)R, 


where AS°=S,°—2S,°. 
constant is given by 


In K=In (Ni V/N2?) =AS°/R+AE/RT, 


where AE is the energy of dissociation at constant 


(17) 


Now the equilibrium 


(18) 

























































volume per mole of (NO)2. If we assume that, 
with change of V, the liquid acts like a perfect 
gas, then AT=AE+PAV=AE+RT, the extra 
energy being taken up by the liquid forces (P is 
internal pressure). Since it is obvious that 


AoS=AnS—A,S (19) 


it may be verified that Eqs. (17) and (18) are 
together equivalent to Eq. (4). This verifies the 
fact that the assumptions made about the liquid 
are equivalent to assuming that NO and (NO). 
form a perfect solution. 

We may now use Eq. (17) together with Eqs. 
(14), (15) and (16) to obtain a_ theoretical 
estimate of A S to compare with the ‘“‘experi- 
mental” value given by Eq. (19). In order to use 
Eq. (17), however, we need to know the value 
of V. An estimate of this quantity. (which may 
without appreciable error be set equal to Vy/2) 
may be obtained on the basis of the assumption 
that the quantity AyS is simply the change of 
entropy due to the change of volume when 
evaporation occurs from the normal unasso- 
ciated liquid, i.e., 


AnS=R In (V,./Vu)=RIn (V,/2V), (20) 


where V, is the molar volume of the vapor 
phase at atmospheric pressure. We may assume 
that wherever they occur in logarithms, V,=N; 
+N2=N/2. We then get from Eq. (17) 


ApS = NiAS®—(Ni/N)R In (V,/N) 
—(N2/N)R In 2+(Ni/N)(AnS—R) 
—(Ne2/N)R In (N2/N). (21) 


In order to evaluate AS® we have to consider 
the quantities involved in Eq. (14). The quan- 
tities pertaining to NO molecules are easily 
obtained. B. can be found immediately; to get 
C, we take’ J,= 1.64 X 10-**; D. = 2, on account of 
A-type doubling, since NO is in a II state (it is 
actually slightly greater than this, because the 
lowest state of NO is a doublet, but the upper 
doublet level is not very important at the boiling 
point of nitric oxide) ;§ Ez we assume to be equal 
to 1, which is a reasonably good approximation 
at this temperature. 


7Sponer, Molekiilspektren, 1 (Springer, Berlin, 1935), 
p. 36. 
8 oy Johnston and Chapman, J. Am. Chem. Soc. 55, 153 
1933). 
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In order to evaluate the quantities for (NO), 
we need to make some approximations about its 
structure. The unpaired electron of the NO 
molecule is a promoted 7 electron.’ The wave 
function presumably has two nodal planes, one 
containing the line of centers of the atoms, and 
one perpendicular to that line. The energy of 
dissociation of (NO)2 is probably sufficiently 
small so that the wave functions are not greatly 
distorted when the double molecule is formed. 
The formation of the double molecule presumably 
takes place in such a way that there is maximum 
overlapping of the wave functions, giving either a 
straight line structure O—N—O-—N or a square 
O-—N O-—N 

| |,or | |. (The selection of O—N—O-—N 
N-O O-N 

rather than N—O—O—N, or O-N-—N-(Q, is 
indicated by the fact that the thermodynamic 
entropy of solid nitric oxide does not become 
zero at 0°K, but approaches a value of & In 2 per 
molecule of (NO)2.; from this it seems likely that 
(NO)s has two distinct possible orientations of 
practically equal energy in the solid state.) We 
shall try both assumptions. 

Let us first assume the square molecule. The 
translational factor, B, is again easily obtained. 
To get C, we divide it into two parts, C;’ and 

i’, the first involving the rotation in the plane 
of the square, and the second involving the 
rotations about axes in that plane. We first note 

O-—N 
that if thestructureis | | o;’=2, but since the 

N-O 
square is probably not actually a true square, 
but a rectangle, we take o,’’=1. (If the other 
square structure is correct o;/=1 and o,"=2, 
which gives the same final result.) However, it 
should be sufficiently good to assume it to be a 
true square with the normal NO distance in 
order to make an estimate of the moments of 
inertia. We thus get J;’=4J_ and J,’’=2J/s. We 
assume D,=1. E, involves six translational 
degrees of freedom. Two of these will be like 
those of NO, and hence will introduce factors of 
1 in #;. The other four frequencies are un- 
doubtedly low, and we shall take E,=2!, which 
. » Mulliken, Rev. Mod. Phys. 4, 1 (1932). The rest of this 


paragraph is based on suggestions of Professor Mulliken, 
to whom I am grateful for discussing this problem with me. 
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corresponds to frequencies of approximately 
60 cm~. 

However, a correction to the above procedure 
is to be noted. The calculation of C; involves 
implicitly the assumption that the molecule 
(NO): can rotate freely in the liquid, which is 
probably not true. In order to correct for this we 
shall assume that the three degrees of freedom of 
rotation are restricted in the same way and to the 
same extent as the three translational degrees of 
freedom. In other words, in order to calculate C; 
we shall first multiply by the correction factor 
V/Vo, where Vo is the actual physical volume of 
the liquid, to be obtained from the observed 
density of the liquid.’ This gives V/Vo= Vu /2Vo 
=().058. Making this correction, then, we obtain 


AyS= —8.4. (22a) 


If the straight line model for the double 
molecule is correct there is one less rotational 
degree of freedom and one more vibration. We 
therefore take E,;=2°, and we shall calculate C, 
on the basis of the rough estimate 7;= 15/2. It is 
perhaps a little doubtful whether with only two 
degrees of freedom of rotation we should still 
correct C; by the same factor V/V )=0.058. 
However, we shall do this, thus obtaining 


AoS= —8.7 (22b) 


not very different from the value for the square 
model. 

Considering the uncertainties necessarily in- 
volved, these values seem to be in reasonably 
good agreement with the ‘“‘experimental”’ value, 
10.8. The discrepancy indicates that the motion 
of the double molecules is relatively more 
restricted than we have assumed. 


'° lt may possibly be that we are exaggerating this effect 
for the double molecule, square model, as compared with 
the single molecule, which is assumed to rotate freely in the 
liquid. For it is true that the extension in space of the 
square double molecule is not much greater in any direction 
than that of the single molecule, and so it may be that the 
double molecule should be assumed to rotate more freely 
than allowed by our assumptions. If this is the case the 
discrepancy between the theoretical and experimental 
values of ApS is increased. This may, perhaps, be an 
argument in favor of the straight line model. While, as 
noted, we may be guilty of exaggerating somewhat the 
effect of the neighbors on the rotation in that case also, 
the exaggeration cannot be very great, as considerable 
inhibition of rotation seems certain to occur. It should, 
perhaps, be remarked that this correction is completely 
consistent with our previous assumptions about the liquid 
only if it is independent of the amount of association in 
the liquid. 
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From the measurements of Johnston and 
Weimer’ or Eucken and d’Or,' it seems reasonable 
to infer that at the boiling point and atmospheric 
pressure the gas is about 1 percent associated. 
Using this estimate, we can apply Eq. (18) to 
calculate AE or AH for the gas phase. We shall 
insert into Eq. (18) the value of AS° from Eq. 
(19),* corrected, however, for the fact that in the 
gas phase the rotation of the double molecule is 
not restricted by the proximity of other mole- 
cules. This gives AE = 2630 calories per mole of 
(NO)s or AH = 2870. This differs from the value 
obtained for the liquid by about 1000 calories. 
It is perhaps, of some interest to note that the 
values for liquid and gas agree very well if we do 
not correct A.S° to allow for the free rotation in 
the gas. This suggests that possibly there is free 
rotation in the liquid, and that the restriction of 
the motion of the molecules is due to their 
structure and applies also to the gas phase. The 
theoretical calculations speak against this, how- 
ever, and it seems probable that the difference 
between the liquid and gas phases is real, though 
there is still the possibility that it is brought in 
by some of the assumptions and approximations 
we have made. These all seem quite reasonable, 
however, and it seems to me rather unlikely that 
all of the difference can be ascribed to this 
cause. It is, perhaps, not particularly astonishing 
in any event that the liquid forces should make a 
difference of 1000 calories in the heat of 
dissociation. 

lt is a pleasure to acknowledge the benefit of 
discussions with Professor J. H. Hildebrand. 

It may be well to add a few words regarding the con- 
nection of this article with another recently published." 
It may now be considered fairly definitely settled that the 
energy of dissociation of (NO): is sufficiently great so that 
it may well act as an intermediate in the trimolecular 
oxidation of nitric oxide. The opportunity is also taken 
to mention some earlier papers, which were overlooked in 
my work on the oxidation of nitric oxide. A number of 
years ago Baker” and Hasche"™ presented evidence that this 
reaction is more complicated than usually supposed. If this 
is true the discussion of the points brought out in my 


paper is somewhat less urgent, though the principles laid 
down are naturally unaffected. Bodenstein,!* however, 


*I.e., the value consistent with the ‘‘experimental”’ 
value, 10.8, of AoS. 

1 Rice, J. Phys. Chem. 4, 53 (1936). 

12 Baker, J. Chem. Soc. (London) 65, 611 
Hasche, J. Am. Chem. Soc. 48, 2253 (1926). 

13 Bodenstein, Zeits. f. physik. Chemie B20, 458 (1933). 


(1894). 
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believes, on the basis of his own experiments and those of 
his co-workers, that the reaction is simple—at least, he 
seems to feel that complications of the type envisaged by 
Baker and Hasche do not occur. More recently Boden- 
stein! has stated his belief that (NO). is an intermediate 
and has cited the work of Eucken and d'Or! and Johnston 
and Weimer’ as indicating its existence. He has also men- 
tioned experiments which indicate that the change of 


14 Bodenstein, Helv. Chim. Acta 18, 745 (1935). 
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temperature coefficient at low temperatures found by 
Briner, Pfeiffer and Malet!® and discussed in my article 
does not actually exist. 

It should, perhaps, be remarked that this correction is 
completely consistent with our previous assumptions about 
the liquid only if it is independent of the amount of associ- 
ation in the liquid. 


15 Briner, Pfeiffer and Malet, J. Chim. Phys. 21, 25 
(1924). 
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HE dielectric constants of highly polar sub- 
stances may not be determined in solution 
because of the polarity of the solvent necessary, 
nor in the vapor state at high temperatures 
because of the conductivity of the vapors, 
especially when in contact with metal surfaces. 
The only method available is the molecular ray 
method, and this method has been applied by 
Stern! and his collaborators to a number of 
substances under conditions which were not 
favorable to precise determinations of the mo- 
ment. Because of the complicated situation 
produced by the superposition of the transla- 
tional and rotational distributions on the space 
quantization of rotating dipoles in an electro- 
static field, a very exact quantitative measure- 
ment of beam intensity is required for a calcula- 
tion of the moment. Some years ago, one of the 
authors? found that the Taylor* hot filament 
detector worked as well for alkali halides as for 
alkali metals, and it was at once apparent that 
it would be possible to determine the moments of 
the alkali halides by the beam method. This re- 
search has been continued in various forms of 
apparatus for several years. The moments ob- 
tained were always of the same order of magni- 
tude, but there were so many sources of uncer- 
tainty that publication has been delayed until 
these were eliminated so far as possible. A 
preliminary report of this work was given at the 

1], Estermann and R. G. J. Fraser, J. Chem. Phys. 1, 
391 (1933). 

2W. H. Rodebush and W. F. Henry, Phys. Rev. 39, 


386 (1932). 
3 J. B. Taylor, Zeits. f. physik. Chemie 57, 242 (1929). 
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Oxford meeting of the Faraday Society in 1934. 
Since then, Schefferst has published results for 
three salts obtained by using a method similar 
to ours. Our results are not in agreement with 
those of Scheffers. Possible reasons for the dis- 
crepancy will be discussed later in this article. 


THE POLAR CHARACTER OF SALT MOLECULES 


It may bé¢ predicted from general considera- 
tions of a chemical character and perhaps more 
definitely from wave mechanics considerations 
that the alkali halides should have a higher 
dipole moment than is ordinarily observed in 
diatomic molecules. If one assumes, as_ has 
usually been done, that this moment is of the 
order of magnitude of 10 Debye units, then a 
very simple consideration shows that a con- 
siderable tendency to association should exist. 
The situation in the saturated salt vapor at 
temperatures below the melting point of the salt 
is analogous to the situation in dilute solutions 
in a solvent of low dielectric constant, where 
Kraus’ has found evidence of association. From 
familiar analogies one would expect this associa- 
tion to be limited to double molecules which 
would have zero moment, and one is thus led to 
the paradox that a high dipole moment should 
not be observed at all. This conclusion is evi- 
dently contrary to fact, and it is not possible to 
make any predictions as to association or how 
it should vary with the temperature of the 

4H. Scheffers, Physik. Zeits. 35, 425 (1934). 


®R. M. Fuoss and C. A. Kraus, J. Am. Chem. Soc. 57. 
1 (1935). 
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saturated vapor. It is possible that the associated 
molecules may have a greater moment than the 
single molecules. The moments reported here are 
the apparent moments as measured, and the 
question of association will be discussed else- 


where. 
METHOD AND APPARATUS 


The experiment consists in passing a beam of 
salt vapor through a Rabi type of electrostatic 
field and mapping the beam, with the field on 
and off by means of the positive ionization 
detector. An elaborate description of the appa- 
ratus is not necessary, since the various con- 
stituent parts of the apparatus have been de- 
scribed elsewhere. A schematic drawing of the 
furnace slit system field plates and detector is 
shown in Fig. 1. In order to simplify the problem 
of alignment, the furnace, slit system, and field 
plates were mounted upon a metal base which 
was removable as a unit from the vacuum 
chamber. The vacuum chamber consisted of a 
brass tube of 42 in. internal diameter and 30 in. 
length. The ends of this tube were closed by 
glass plates which were ground to give a vacuum- 
tight joint when lubricated with a suitable stop- 
cock grease. The vacuum was produced by glass 
pumps, using dibutyl phthalate as a vapor. 
These pumps were provided with the high speed 
turbine type of nozzles designed by Drs. Copley 
and Phipps’? in this laboratory. With these 
pumps it was possible to maintain at all times a 
pressure of about 10-® mm of mercury. 

The furnace of tungsten steel was similar to 
the one designed by Lewis.* Tungsten coils were 
suspended in holes drilled in the walls of the 
furnace, the heat from the coils being trans- 
mitted by radiation. Some difficulty was ex- 
perienced in early work with electrolysis of the 
salt. In order to avoid this, a special removable 





_—_— 


*R. G. J. Fraser, Molecular Rays (Cambridge Uni- 
versity Press, 1931). 
1938) J. Copley and T. E. Phipps, Rev. Sci. Inst. 6, 265 
(1935). 

*L. C. Lewis, Zeits. f. physik. Chemie 69, 786 (1931). 
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crucible was placed inside the furnace. This 
crucible was in contact with the furnace only at 
its base. The temperature of the furnace was 
controlled by a thermocouple, photoelectric cell, 
thyratron tube, set-up which proved to be highly 
satisfactory. In measuring intensities with the 
field, off and on, it is of the utmost importance 
to avoid temperature fluctuations. The tempera- 
ture variations were ordinarily less than 0.2 
degree. The accuracy with which the absolute 
temperature is known is of less importance for 
the measurement of the moments. 

The plates for the Rabi type field were 
mounted 1.65 mm apart on quartz rods of 
sufficient rigidity to maintain the separation 
against the mechanical forces due to the electro- 
static field. The electrostatic field was produced 
by the customary arrangement of transformer, 
kenotron and condenser. Elaborate tests were 
made to demonstrate the constancy of the 
field. 

The detector consisted of a tungsten filament, 
cylindrical plate, and high sensitivity galva- 
nometer. The detector filament was mounted 
eccentrically on a ground-glass joint, and the 
position of the filament was indicated by a 
protractor arm of 83 in. length. The filament was 
0.05 mm in diameter. So far as the conditions of 
this experiment are concerned, two slits are 
sufficient to define the beam, but it was found 
simpler to use three slits in actual practice. In 
the simple design of furnace used, the outside of 
the slit was somewhat cooler than the inside, so 
that difficulty was experienced in keeping the slit 
free from deposited salt when the slit was small 
enough to define the beam. Accordingly, an 
intermediate slit, designated as a ‘“‘heated slit,”’ 
was introduced. This slit was cut in a nickel 
cylinder on the inside of which a tungsten fila- 
ment was suspended so that the slit could be 
heated at frequent intervals to drive off salt 
deposits. It proved very satisfactory in operation. 

The distance between the second and third 
slit was 5.7 cm. The distance from the third 
slit to the point of incidence of the field was 
1.5 cm and the distance from the point of 
incidence to the detector filament was 7.6 cm. 
A magnetically controlled shutter was provided, 
so that the beam could be turned off and on 
at will. 
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The detection of beam intensity by the positive 
ion current is enormously complicated by the 
presence in the apparatus of a somewhat un- 
steady high potential on the field plates, which 
inevitably must affect the sensitive galvanometer 
required for detection. This difficulty was 
avoided by using the method introduced by 
Rabi® in which the salt molecules are deposited 
on a cold filament for a measured period of time. 
At the end of this period, the shutter is closed, 
the field on the plate switched off, and the fila- 
ment is flashed, the resulting positive ion current 
being then measured by the ballistic throw of 
the galvanometer. Experiment shows that, when 
more than about five percent of the filament 
surface is covered, not all the alkali atoms come 
off as ions, so that the deflection is no longer 
proportional to the time of exposure. The limiting 
time of exposure, for which the deflections were 
proportional to the time, varied, as might be 
expected, with the size of the salt molecule for a 
given beam intensity. Most of the readings were 
made with an exposure period of three minutes. 
The galvanometer had a ballistic sensitivity of 
4x10-". 


THEORETICAL CONSIDERATIONS 


The theory of the deflection of diatomic 
molecules in an inhomogeneous electrostatic field 
has been discussed by Fraser® and by Fraser and 
Estermann.” One is dealing here with a super- 
position of three distributions of the molecules in 
the field. These distributions are with respect to 
linear velocity, rotation and space quantization 
in the field. The moment that is effective in the 
deflection is the moment averaged over a com- 
plete rotation. Feierabend has made the very 
laborious calculation of the form of the dis- 
tribution curve for a beam of zero width. In the 
actual experiment, however, we find that there 
are two further complications; the finite slit 
width and the finite width of the detector 
filament. 

The latter quantity has an important bearing 
upon the design of the slit system, and will be 
considered first. When the beam is defined by 
two slits identical in size and shape, the intensity 

9]. 1]. Rabi and V. W. Cohen, Phys. Rev. 46, 707 (1934). 


10], Estermann and R. G. J. Fraser, J. Chem. Phys. 1, 
390 (1933). 
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in the beam plotted against a space coordinate 
taken perpendicularly to both beam and slits 
forms a trapezoid (Fig. 2). The flat part of the 
trapezoid cérresponds to the umbra, and the 
sides to the penumbra. This is for the undeflected 
beam. With the field on, the intensity curve is 
altered in a way that may be approximately 
described as a trapezoid of lower altitude and 
broader base. It turns out that the most im- 
portant measurements are the altitudes of the 
trapezoids in the two cases. In order that the alti- 
tude be measured with certainty, the flat part 
of the beam should be at least as wide as the de- 
tector filament. In order to obtain the maximum 
sensitivity, which means the maximum decrease 
in the altitude of the trapezoid when the field 
is on, the beam should not be wider than this. 
The best conditions, therefore, appear to be 
obtained when slits and detector filament are 


of the same width. In our experiments the slits 


were approximately 0.05 X0.5—0.8 mm. 

As the detector filament is moved through the 
beam, the “‘corners” of the trapezoid are rounded 
off and ‘‘tails’’ are produced by the overlapping, 
so that the actual curve of galvanometer de- 
flections obtained is that corresponding to the 
circles (Fig. 2). If the geometry of the slit system 
is known with sufficient exactness, the trape- 
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zoidal curve may be inferred from the actual 
curve obtained, but the only points on the two 
curves that correspond are the maximum and 
the mid-portions of the sides. If a sufficient 
number of readings are taken, a satisfactory 
correspondence may be established between the 
curve of galvanometer deflections and the true 
trapezoidal intensity curve, but experience has 
shown that the true intensity curve may be 
calculated from the geometry of the slits, pro- 
vided they are carefully constructed, with an 
accuracy quite equal to that obtained by careful 
mapping of the actual beam by galvanometer 
deflections. The measurement on the unde- 
flected beam is, therefore, in practice limited to 
the measurement of the maximum intensity. 

Similarly, for the deflected beam, it turns out 
that the only measurement of significance is the 
maximum intensity. Since, in deflection experi- 
ments of this character, the two maxima ap- 
parently coincide in space, when once the 
position of the maximum has been located a 
series of readings may be taken with the field 
alternately on and off without moving the 
filament. 

The intensity curve of the deflected beam may 
be considered as a distorted trapezoid lowered 
and broadened with sides which are no longer 
straight lines. 

The upper corners of the trapezoid will also 
be rounded off so that the top of the trapezoid 
is no longer flat, and the maximum reading 
obtained with the detector filament is not the 
actual maximum intensity of the beam, but 
something less. It is not possible to reconstruct 
the intensity curve, therefore, from galvanometer 
deflections with any certainty, except in the 
edges of the beam where the intensity is too low 
to permit accurate determinations. It is possible, 
however, to assume an approximate value of the 
moment and calculate the exact shape of the 
curve in the neighborhood of maximum intensity 
from theoretical considerations. The exact shape 
of the curve being known, it is, therefore, 
possible to calculate the exact maximum in- 
tensity from the galvanometer reading, and 
hence the true value of the moment as a second 
approximation. In our experiments, the actual 
galvanometer readings of maximum intensity 
were shown to be about one percent too low. 


THE ALKALI HALIDES 

The general shape of the intensity curve with 
the field on is given by the dotted line in Fig. 2. 
This shows the general characteristics of similar 
curves obtained in beam experiments in that the 
upper part of the curve does not appear changed 
in shape, but only decreased in height. 

The dipole moment is, therefore, calculated 
from the percentage decrease in maximum in- 
tensity of the beam when the field is applied 
and the accuracy depends upon the accuracy 
with which this percentage decrease is measured. 
The percentage decrease may be made greater in 
magnitude by decreasing slit width and in- 
creasing the distance between the field and de- 
tector filament, but this change is accompanied 
by a decrease in the magnitude of the deflection 
so that there is no advantage in changing 
dimensions of the apparatus beyond a certain 
point. Similarly, increasing the sensitivity of the 
galvanometer does not increase the reliability of 
the measurements. A possible way of increasing 
the precision of the measurements would be to 
increase the strength of the field on the Rabi 
plates, but there are practical limitations 
here too. 

We have finally to consider the effect of finite 
slit width on the calculation of the moment. 
The curve calculated by Feierabend, previously 
referred to, gives the distribution of the intensity 
in a beam from slits of infinitesimal width. We 
may represent the Feierabend curve as a function 
of the deflection in mm P,,). The intensity curve 
of the undeflected beam may be represented in 
a similar way as f(,), where s is now the displace- 
ment of the detector filament in mm. It has been 
shown by Estermann” that the intensity at the 
maximum of the deflected beam is given by 


the integral 
400 
r= fis) P sds. 


The integral must be obtained by quadrature, 
since the function P,,) is not known. The values 
for P;.), as given by Feierabend, are not given for 
actual deflection, but for an argument, 


o=s/So, 


where s is the actual deflection, and so is the 
deflection that a molecule would experience when 
in certain quantum states that are taken as an 
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DaTA FOR SAMPLE RuN KI 





.901°K 
83,700 volts/cm 


Temperature 
Field intensity 
Average maximum intensity 


undeflected beam 64.8 
Average maximum intensity 
deflected beam 55.3 
Distance field plates to detector 76 mm 
6 7° 48’ 
So 0.00263 mm 
mn 11.1X10718 Ls, 











TABLE OF RESULTS 

SALT TEMPERATURE uX 1018 RESULTS OF SCHEFFERS 
KCl 949°K 9.74 | 

“ 9.74 | 

” 9.25 } 6.3 (1023°K) 

. 9.40 | 

9.54 } 
KBr 916 10.8 

925 10.9 

KI 894 11.0 ) ° 
CsI 873 12.1 








average or standard state. s» is actually given by 
the equation 


So=p?E*l/8(kT)? tan 6, 


where kT is the Boltzmann factor, E the field 
intensity, / the distance from the point of in- 
cidence to the detector filament and @ is the angle 
between the incident beam and the field plates. 

uw is, of course, the dipole moment, and the 
calculation of the maximum intensity of the de- 
flected beam is made by assuming different 
values of uw until one is obtained which gives a 
value in agreement with the experimental 
results. 

By reference to Table I, it is seen that the 
values as calculated from our experiments are 
much greater than those obtained by Scheffers. 
The order of arrangement of the salts according 
to moment appears to be the same in the two 
experiments. A rough prediction of this order 
could be made from a consideration of the ionic 
radii. 

Scheffers’ results are at higher temperatures 
and pressures than ours. The discrepancy be- 
tween the two sets of experiments could be 
accounted for by association. Until the question 
of association is settled, one may not draw 
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conclusions as to the actual moment of the 
molecules. 

Scheffers used an inhomogeneous electrostatic 
field produced by a wire and concentric cylinder. 
This type field gives a greater deflection, but it 
is necessary to locate the beam with exactness 
with respect to the wire. The Rabi type of field 
does not introduce any serious uncertainty so 
long as the angle of incidence is not made too 
small, but it does not give as large a deflection. 
We have tested the Rabi field for ‘‘end effects” 
by changing the point of incidence of the beam. 
No change in the apparent moment was ob- 
served. We have not checked the calculation of 
Fraser that the effective square of the electric 
moment is 0.6 times the field between the plates. 
If this conclusion is in error, then our results 
would be in error by a corresponding amount. 

It is difficult to estimate errors in a measure- 
ment which involves so many variables in such 
a complicated way, but with a great many 
measurements made under different conditions 
we have never obtained more than a ten percent 
variation in the effective moment. 

In every case we tested the beams for ions by 
applying a ldw field to the Rabi plates, but only 
in the case of cesium fluoride did we observe an 
appreciable diminution of intensity. 

The purity of the salts used was tested. They 
were shown to be free from appreciable amounts 
of impurities. 

The fact the dipole moments obtained are of 
the order of magnitude expected and increase 
with the ionic radii may be considered as evi- 
dence favoring the conclusion that the salt 
vapors contain single molecules only. The density 
of salt vapors has been investigated at some 
length in this laboratory by Dr. V. Deitz, and the 
results of his investigations will be published 
elsewhere. 

The dipole moment obtained indicates that a 
mixture of the ionic and homopolar types of 
bonding prevails in salt molecules. If the bonding 
were purely ionic in character, a still larger 
dipole moment would be anticipated than has 
been observed. 

The authors wish to acknowledge their in- 
debtedness to Drs. M. J. Copley and T. E. 
Phipps of this laboratory for many helpful 
suggestions. 
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Determination of Molecular Structure by Electron Diffraction 


Puitie G. ACKERMANN* AND JOSEPH E. Mayer, Chemical Laboratory of The Johns Hopkins University 
(Received February 21, 1936) 


(1) An apparatus has been designed and built for the 
purpose of determining molecular structure by electron 
diffraction using 6400 volt electrons. (2) The apparatus 
was checked with carben tetrachloride whose structure is 
accurately known, and the interatomic distances in, and 
configurations of, a number of compounds determined. 
(3) Ethylene oxide was found to be a three membered 
ring with a carbon-carbon distance of 1.56+0.05A, a 
carbon-oxygen distance of 1.45+0.05A and a carbon- 
hydrogen distance of 1.05+0.07A. (4) Sulfur monochloride 
was found to have the two chlorine atoms attached to 


different sulfur atoms, the sulfur-sulfur distance was 
found to be 2.04+0.05A, the sulfur-chlorine distance of 
1.98+0.05A and an angle of 105+5° between bonds. 
(5) In acetaldehyde the carbon-oxygen distance was found 
to be 1.20+0.05A, the carbon-carbon distance 1.51+0.05A 
and the bond angle 122+5°. (6) Paraldehyde was found 
to have a puckered ring structure with a carbon-oxygen 
distance of 1.41+0.05A, a carbon-carbon distance of 
1.52+0.05A and the tetrahedral angle of 109° between 
bonds. 





INCE the original work of Mark and Wierl,! 

the results of many investigators? have 
shown that electron diffraction methods provide 
a very useful tool for determining molecular 
structure. The object of the experiments re- 
ported here was to extend this work and to 
determine if much lower voltages (6500 volts 
instead of 20 to 50 kilovolts) might not be used 
in electron diffraction apparatus for determining 
molecular structure. 

In the interpretation of the experimental 
results the complete formula for the coherent 
scattered intensity, given below, was not used. 
In this formula 





n sin Xi; 
I« y vivji— ’ (1) 
i, Xij 
where ¥i=(Z:—F;)/p? 
and xij= (4m sin $6/d)) 5; =ul;; 


Zand F being the atomic number and form 
factor for the ith atom, /;; the distance between 
the ith and jth atoms, the summation being 
taken over all the atoms of the molecule. As 
Pauling and Brockway’ have pointed out, when 





_* Adapted from a thesis submitted to The Johns Hopkins 
University in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy. 

1 Mark and Wierl, Physik. Zeits. 31, 366, 1028 (1930). 

* Especially Pauling and his co-workers: Pauling and 
Brockway, Proc. Nat. Acad. Sci. 19, 860 (1934); Brockway 
and Wall, J. Am. Chem. Soc. 56, 2373 (1934); Sutton and 
Brockway, J. Am. Chem. Soc. 57, 473 (1935); and also 
Dornte, J. Chem. Phys. 1, 566, 602 (1933) and de Laszlo, 
Proc. Roy. Soc. A146, 672 (1934). 

* Pauling and Brockway, J. Chem. Phys. 2, 867 (1934). 
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the maxima and minima of intensity are meas- 
ured by visual means, as was done here, the 
simpler formula, given below, usually suffices. 
n sin X;; 
/ % 
ae ) 22;-— (2) 


i, j=1 Xij 


where J’ =/-f(@) and f(@) is a monotonic function 
of the angle, 6. The difference between the two 
formulae is greatest at small angles, and pre- 
liminary calculations with a model of ethylene 
oxide similar to that finally used gave for the 
positions of the first maximum by the two 
formulae of 4.17 and 4.21X10~°, the difference 
being less than the experimental error. From the 
experimental curves given by Mott and Massey‘ 
we may conclude that the scattering by single 
atoms is still proportional to ((Z;— Fi)/u?)? for 
voltages as low as 2000 volts except for very 
heavy atoms or very large angles of scattering. 
Hence we may feel sure that our formula (1) 
and consequently (2) will be valid for the 6400 
volts used. 


DESCRIPTION OF APPARATUS AND METHOD 


The general design of the apparatus follows 
that of Hendricks,’ except that a filament was 
used as a source of electrons rather than a dis- 
charge tube. The filament was a flat spiral of 
platinum wire coated with barium oxide, and 


4Mott and Massey, The Theory of Atomic Collisions 
(Oxford University Press, 1933), p. 122. 

5 Hendricks, Maxwell, Moseley and Jefferson, J. Chem. 
Phys. 1, 549 (1933). 
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was so mounted that it could be moved from 
outside the apparatus to give the best electron 
beam intensity. Eastman process plates were 
used to record the scattered electrons, but 
because of the slight effect of the low voltage 
electrons on the plates it was found advisable to 
oil them. The distance from the jet through 
which the gas issues to the photographic plate 
was 7.45 cm. The whole apparatus was mounted 
so that the path of the electron beam made an 
angle of 45 degrees with the horizontal and was 
in the direction of the earth’s magnetic field so 
as to reduce the effect of this field on the beam as 
far as possible. 

The high potential was secured from a 110- 
7000 volt transformer using two General Electric 
PJ-24 tubes for full wave rectification. The 
pulsations were smoothed out by means of a 0.5 
microfarad condenser and a choke coil. The 
filament was heated by the current from another 
transformer and the filament temperature con- 
trolled by a variable resistance in the primary 
circuit of this transformer. The voltage applied 
across the apparatus was determined by meas- 
uring the current through a calibrated high 
resistance of eight megohms. 


EXPERIMENTAL RESULTS 


As a check on the apparatus and in order to 
become familiar with the interpretation of the 
actual photographs, it was decided to first study 
a compound whose structure was already known. 
Carbon tetrachloride was chosen for this because 
of the accurate data which are available for this 
compound. The C. P. material was distilled under 
a vacuum into the bulb from which it was to 
be vaporized and there several times alternately 
frozen and melted while under a vacuum to 
eliminate as much as possible of the dissolved air. 
During the run the material was kept at 0°C 
by means of an ice bath, and exposures of from 
4 to 32 seconds were taken. The results are tabu- 
lated in Table I below. The experimental figures 
represent the average of measurements on twelve 
plates, the deviation of any measurement from 
the average being not more than two percent. 
The accelerating voltage was 6400 volts which 
corresponds to a wave-length of 0.1531A. The 
calculated values are based on a chlorine-carbon 
distance of 1.77A which compares very favorably 
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TABLE I. Carbon tetrachloride. 








VALUES OF yp X 108 





Observed Calculated 
Maximum 2.41 2.69 
Minimum 3.76 3.74 
Maximum 4.80 4.79 
Minimum 5.97 5.96 
Maximum 7.16 7.16 








with the value of 1.76A found by Brockway and 
Wall? and by de Laszlo and Cosslet® and the 
value of 1.78A found by Braune and Knoke.’ 
Pauling and Brockway* have indicated that 
the interatomic distances measured by electron 
diffraction can be used to calculate a series of co- 
valent radii for the various atoms. Because of 
the fact that the structural formula for ethylene 
oxide is usually written as a three membered ring, 
it was thought that in this compound the carbon- 
oxygen distance might be slightly different from 
that in acyclic compounds. The ethylene oxide, 
taken from a commercial tank, was carefully 
distilled, the fraction taken for use having a 
boiling point range of less than one degree. 
The compound was kept at —40°C during the 
run, this temperature being found to give the 
optimum vapor pressure. Because of the small 
scattering power of ethylene oxide, exposures up 
to 300 seconds were necessary. The experimental 
values of uw are given in Table II and are the 


TABLE II. Ethylene oxide. 








VALUES OF yp X 108 
Observed Calculated 





Maximum 5.34 5.37 
Minimum 7.46 7.48 
Maximum 9.58 9.55 








average of measurements on nine plates. The cal- 
culated values are based on a carbon-oxygen 
distance of 1.45A and a carbon-carbon distance 
of 1.56. In compounds of low scattering power it 
may be necessary to take the scattering by the 
hydrogen atoms into account ; this was done here 
using a carbon-hydrogen distance of 1.05A and 
although the calculated values were not very 
sensitive to small changes in this distance, in- 
cluding these terms does have an appreciable 


6 de Laszlo and Cosslet, Nature 134, 63 (1934). 


7 Braune and Knoke, Zeits. f. physik. Chemie B21, 297 
(1933). 
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Fic. 1. Theoretical intensity curve for ethylene oxide. 


effect. The theoretical intensity curve is given in 
Fig. 1. The carbon-carbon distance found above 
is in good agreement with the distance of 1.54A 
found by Dornte.? The carbon-oxygen distance 
compares quite favorably with the distance of 
1.444 found by Sutton and Brockway?’ in 
dimethyl ether. These results, together with the 
fact that no alternative structure could be found 
which would give any agreement with the experi- 
mental points, indicates that ethylene oxide is a 
three-membered ring in which the interatomic 
distances are nearly the same as in acyclic 
compounds. Alternate structures in which the 
oxygen atom was at a greater distance from one 
carbon atom than from the other gave no agree- 
ment with experiment. 

The structure of the compound, sulfur mono- 
chloride, is not definitely established. The avail- 
able chemical evidence does not seem to be 
sufficient to distinguish between the two alter- 
native structural formulae given below, though 
it does seem to favor formula IT. 


cl 
Nga 
DSs=S I 


Cl 
Ns-S, Il. 
Cl Cl 


For this reason it was decided to study this 
compound by electron diffraction to see if more 
light could be thrown on the subject. Merck’s 
sulfur monochloride was carefully distilled in a 
vacuum from a receiver containing a quantity 
of free sulfur and activated bone charcoal to 
remove impurities. The material was kept at 
30°C during the runs and exposures of from 10 to 
200 seconds were taken. The results, given below, 
are the average of measurements on eleven 
plates. The calculated values were obtained by 
using formula II, a chlorine-sulfur distance of 
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1.98A, a sulfur-sulfur distance of 2.04A and an 
angle of 105° between bonds. In the structure IT 
for sulfur monochloride it is probable that there 
would be some free rotation about the sulfur- 
sulfur bond. Accordingly with the use of the 
above structure, distances, and angles, the theo- 
retical curve was calculated for three different 
positions of the two chlorine atoms. In Table ITI, 


TABLE III. Sulfur monochloride. 








VALUES OF yu X 108 





Observed Calculated 
A B C 
Maximum 4.12 4.20 4.13 4.17 
Minimum 5.35 5.44 5.50 5.44 
Maximum 6.68 6.64 6.70 6.63 
Maximum 7.83 8.10 7.86 7.60 
Minimum 8.90 9.04 8.93 9.07 
Maximum 10.21 10.27 10.16 10.32 








A gives the values of » when the shortest possible 
chlorine-chlorine distance consistent with the 
above structure was used, B the results when an 
average distance based on free rotation was 
used, and C the results when the longest possible 
distance was. used. It will be noted that B agrees 
best with the experimental results. When similar 
calculations are carried out using different values 
for the sulfur-sulfur and sulfur-chlorine distances, 
the results obtained do not agree at all with 
experiment for any relative positions of the 
atoms. The theoretical intensity curve for this 
structure is given in Fig. 2. 

While these results are not in such good agree- 
ment with the values of 2.05 and 2.10A, respec- 
tively for the two distances which one would 
obtain by the use of Pauling and Huggins® 
table of covalent radii, they are in favor of the 
second of the two formulae given for sulfur 
monochloride. 

Acetaldehyde and its polymer, paraldehyde, 
were also studied. The acetaldehyde was pre- 
pared by distilling it from C. P. paraldehyde to 
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Fic. 2. Theoretical intensity curve for sulfur monochloride. 


8 Pauling and Huggins, Zeits. f. Krist. 87, 205 (1934). 
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Fic. 3. Theoretical intensity curve for acetaldehyde. 


which a few drops of sulfuric acid had been 
added. The runs were made in the usual manner. 
The material was kept at — 35°C in the bulb and 
exposures up to 300 seconds were taken. The 
results, Table IV, are the averages of measure- 
ments on ten plates. The theoretical intensity 
curve is given in Fig. 3. The calculated values 
were obtained by using a carbon-carbon distance 
of 1.51A, and oxygen-carbon distance of 1.20A 
and an angle of 122° between the two bonds. 
Due to the fact that only two maximum could 
be measured, the structure could not be deter- 
mined uniquely, but that given does show good 
agreement with the experimental measurements 
and with other measurements of these distances. 

For paraldehyde the runs were made in the 
usual manner, the material being kept at 30°C 
during the runs and exposures up to 200 seconds 
were taken. The results, Table V, are the aver- 
ages of measurements on eleven plates. 

For the calculation of the values of u tabulated 
in Table V and the theoretical curve given in 
Fig. 4, a structural formula must be assumed for 


TABLE IV. Acetaldehyde. 








VALUES OF » X 108 





Observed Calculated 
Maximum 5.92 5.91 
Minimum 7.60 7.61 
Maximum 8.74 8.72 








TABLE V. Paraldehyde. 








VALUES OF yu X 103 





Observed P Calculated 
Maximum 3.40 3.43 3.44 
Minimum 4.41 4.45 4.20 
Maximum 5.69 5.67 5.79 
Minimum 7.47 7.46 7.47 
Maximum 9,24 9.22 9.16 
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Fic. 4. Theoretical intensity curve for paraldehyde. 


the molecule. The one used is the usual ring 
formula for paraldehyde with a puckered ring, 
i.e., the carbon and oxygen atoms are not all in 
the same plane. The carbon-oxygen distance was 
taken to be 1.41A, the carbon-carbon distance as 
1.52A with the angle between the bonds the 
tetrahedral angle of 109° 28’ 16”. The relative 
positions of all the carbon and oxygen atoms is 
given by the statement that if the plane formula 
. 
- 
0% So 
> | | 
| | 
CH HC, 
CH; SO” CH; 


be taken as the projection of the molecule on the 
plane of the paper, the projected ring would form 
a regular hexagon, and that if the ring carbon 
atoms were in the plane of the paper, the oxygen 
atoms would be a distance of 1.41 Xsin 19° 28’ 16” 
=0.47A above this plane and the three methy! 
carbon atoms a distance of 1.52 Xsin 19° 28’ 16” 
=(0.51A above this plane. The formula used to 
calculate the intensity then becomes 














sin 1.41 sin 2.3932y sin 2.3025u 
1.41 2.3932u 2.30254 
sin 3.6506u sin 2.7186 sin 4.0919p 
+432—__—_—_ 285 _ 
3.6506pu 2.7186u 4.0919u 
sin 1.524 sin 4.7847 
+21 + —— 
1.52u 4.7847 pu 


This structure gives the values of yw for the 
maxima and minima marked A in the table. 
If we use the similar structure in which the 
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methyl carbon atoms are directly below the other 
carbon atoms a distance of 1.52A, we get the 
values of » marked B in the table. The difference 
is not very marked but it is sufficient to make it 
probable that the methyl carbon atoms are, as 
we have sketched them, above the plane of the 
oxygen atoms. 


DISCUSSION OF RESULTS 


While an apparatus has been developed by 
which one may make measurements of inter- 
atomic distances by means of the diffraction of 
electrons of only 6400 volts energy, the accuracy 
of the results is not as great as that of other 
workers using higher voltages. In the experiments 
reported here, the accelerating voltage was 
known with an accuracy of one percent. The 
ammeter and resistance used were thus accurate, 
the choke coil and condenser were sufficient to 
smooth out the pulsations of the rectified 60- 
cycle alternating current to this degree of accu- 
racy, and the actual measurements showed that 
the variations in the input alternating current 
were no greater than this. 

The distance from the jet to the photographic 
plates could be measured quite accurately, the 
error here being certainly less than one percent. 
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The greater error is in the actual measurement of 
the positions of the maxima on the plates. This 
is especially true since the maxima are not sharp 
lines but are rather broad. The fact that the 
measurements from the individual plates do not 
differ from the average by more than 3 percent 
indicates that the maxima can be determined 
fairly accurately. The final values of yw are 
probably accurate to within two percent. The 
relation between the positions of the maxima 
and the values of the interatomic distances is, 
in general, quite complicated, so that for a given 
error in the former it is difficult to estimate the 
error in the latter. 

A real source of error in these experiments 
comes from the fact that only a few maxima and 
minima could be measured due to the limitations 
of the apparatus. When this is true, the fact that 
the observed and calculated values of uw agree 
very well is not a guarantee of like accuracy in 
the values of the interatomic distances. For 
small variations in the structure are best de- 
tected in the positions of the outer maxima, 
those beyond the third. It is difficult to say just 
how much effect this would have on the com- 
pounds studied here, but it is probable that the 
measurements recorded are accurate to +0.05A. 
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It is shown that the three-dimensional Patterson method 
can be so simplified by the use of the symmetry properties 
of the crystal under consideration that its use in deter- 
mining the positions of atoms in crystals is practicable. 
This method is then used to determine the positions of the 
heavy atoms in proustite, Ag;AsS3. The positions of the 
sulfur atoms are found by the use of the covalent atomic 
radii and assumptions concerning bond directions. The 
structure so derived is found to be compatible with 
the observed data. Proustite is found to contain pyramidal 
AsS;’" groups and silver atoms forming two bonds to 
sulfur in almost opposed directions. Pyrargyrite, AgsSbSs, 
is found to have almost the same structure as proustite. 
The analytical statement of the structures is as follows:— 
The space group, C3,°—R3c, is common to both proustite 
and pyrargyrite. The special positions of C3,° are (in 


hexagonal axes) 


(2a) 0,0, Z; 0,0, 3+2Z; : _ 
(6b) X, Y, Z; Y,X-Y, Z; Y-X, X,Z; 
Y, X,34+Z; X,X-Y,44+2Z; Y-X, Y, 3+2: 


and positions derived from these by the operations of the 
rhombohedral lattice. The parameter values are: 


Proustite Pyrargyrite 

ado = 10.74A, co=8.64A ao= 11.04A, co=8.71A 
2Asin 2a, Zs,=0.000 2Sb in 2a, Zgpy=0.000 
6Ag in 6b, Xag=0.246 6Ag in 6b, Xag=0.250 
Vag =0.298 Yag = 0.305 

Zag = 0.235 Zag = 0.210 

6S in 6b, Xs =0.220 6S in 6b, Xg =0.220 
Ys = 0.095 Ys =0.105 

Zs =0.385 Zs =0.355 
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INTRODUCTION 


HE determination of the inner structure of a 
crystal by means of x-rays has been, in 
most cases, a process of trial and error. Once the 
space group of the crystal is known, the direct 
use of the x-ray data is abandoned. A structure 
which satisfies the symmetry requirements of the 
space group is then assumed and the intensities 
of the diffracted x-ray beams are calculated on 
this basis. If these do not agree with the experi- 
mentally determined intensities, another struc- 
ture is assumed and the process is repeated. 
Systematic methods of carrying out this process 
have been devised which are remarkable for the 
rapidity and accuracy with which, in some cases, 
they lead to the elimination of every structure 
but the true one. In other cases, however, the 
systematic trial and error method breaks down 
and either fails completely or leads to the loss of 
much time. 

If a convenient direct method for calculating 
the structure of a crystal from the experimental 
data were available, these inefficiencies and 
failures would, of course, be eliminated. Such a 
method has not as yet appeared. However, a very 
close approximation to such a treatment is 
provided by the Patterson method when used as 
described in the next section. 

The determination of the structure of proustite 
by means of this method is described later in this 
paper and illustrates the advantages and limi- 
tations of the treatment. 


UsE OF THE THREE-DIMENSIONAL PATTERSON 
METHOD 


Patterson! has shown that the three-dimen- 
sional Fourier series 


P(x, y, 2)= x > | Foren |? 


=—o k=—o /=—o 


X cos 2r(hx+ky+/z) 


—where x, y and z are the coordinates of points 
in the crystal lattice measured in units of the 
primitive translations and Fy. is the amplitude 
of the diffracted beam ‘‘reflected’’ from the 
crystal plane with indices (hkl)—represents a 

* Contribution No. 535 from Gates Chemical Labora- 


tory, California Institute of Technology. 
1A. L. Patterson, Zeits. f. Krist. 90, 517 (1935). 
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function having maxima at vector distances from 
the origin equal to vector distances between pairs 
of maxima in the electron density in the crystal. 
Since there is a maximum in the electron density 
at the center of each atom, the coordinates of the 
maxima in P(x, y, 2) can be interpreted as the 
components of interatomic vectors. The evalua- 
tion of such a series at a sufficient number of 
points in the unit cell to show the true form of the 
function P(x, y, z) has never been carried out, the 
calculation involving a prohibitive amount of 
labor. Instead, the projection of P(x, y, z) on a 
plane perpendicular to one of the crystal axes has 
been used. If the projection is upon the plane 
perpendicular to the z axis, the projected function 
p(x, y) is calculated thus 


1 
p(x, y)= [ Pe. y, 2)dz 
“0 


@ @ @~ 


=D 2D 2X /Fury|? 


h=—a@ k=—ce |=—c 


1 
x f cos 2r(hx+ky+l/s)dz 
J0 


~ > > | Finxo)|? cos 2x(hx+ky). 


h=—o k=—oo 


p(x, y), which is represented by a two-dimen- 
sional Fourier series, has been evaluated com- 
pletely by several authors in the course of various 
crystal structure investigations. Its use has 
proved very fruitful, in some cases making 
possible the determination of structures which 
would have been impracticably difficult otherwise. 

The use of p(x, y) has, however, two disad- 
vantages. In the first place, there is a large 
probability that a given maximum is composed 
of the superposition of two or more maxima of 
P(x, y, z) which have projected onto nearly the 
same spot. Secondly, the resolution of maxima 
which are near together is not as good as can be 
obtained by the use of P(x, y, 2), since only 
a small fraction of the terms in P(x, y, 2) are used 
in p(x, y). The function P(x, y, z) would conse- 
quently be preferred if the labor involved in its 
calculation could be reduced. 

The use of the symmetry of the crystal under 
consideration often leads to the desired simplifi- 
cation. For example, let the crystal have a 
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twofold axis (which will be taken coincident with 
the 6 axis). Then, if there is an atom at the point 
(x, y, 2), there will be another atom, crystal- 
lographically equivalent to the first, at (2, y, Z). 
The vector between these atoms has the com- 
ponents (2x, 0, 2z) and consequently there will 
bea maximum in P(x, y, 2) at the point (2x, 0, 22). 
This maximum lies in the plane y=0. There will 
be a maximum in this plane for each (crystal- 
lographically) different kind of atom in the 
crystal. The x and y coordinates of all the atoms 
in the crystal can consequently be found by 
evaluating P(x, y, z) for y=0 only. In this case 
the expression for P(x, y, z) can be simplified as 
follows: 


@w ~o @w 


P(x,0.23= 4 Y ¥ | Farn!? 


h=—c k=—o I=—o0 


cos 2r(hx+/z) 


m@ 


=> > cos 2r(hx+/sz) 


h=—o l=—oo 
“l > | Four |?} | 
| k=—oo J 
Now let x | Fenny |? = Cnr 
k=—oo 


and we have 


@~ 


P(x,0,s)= > . Cy, cos 2r(hx+/z). 


h=—o I=—o 


P(x, 0,2) is a two-dimensional Fourier series of 
the same general type as p(x, y) and consequently 
its calculation is practicable. The same informa- 
tion about the x and z coordinates of the atoms is 
contained in P(x, 0, z) as in p(x, z), but in much 
clearer form; for (a) the maxima in P(x, 0, z) will 
be sharper, in general, since there are usually 
more terms in its Fourier series than in that of 
p(x, z) and (b) all the maxima corresponding to 
interatomic vectors not parallel to the plane 
y=0 appear in p(x,z), but are eliminated in 
P(x, 0, z). P(x, 0, z) can be used to determine the 
x and z coordinates of atoms in case the b axis is 
parallel to a two-, three-, four-, or sixfold 
symmetry axis of the crystal, for in all these 
Cases the vectors between equivalent atoms are 
parallel to the plane y=0. 


In case the } axis is a twofold screw axis the 
treatment of P(x, y, 2) is slightly changed. The 
coordinates of equivalent atoms are now (x, y, 2) 
and (Zz, y+3, Z) and the vectors between equiva- 
lent atoms have components (2x, —}3, 22). In this 
case P(x, y, z) must be evaluated for y= —} or, 
since P(x, y,z) has a center of symmetry, for 
y=43. The calculation simplifies as follows: 


@o @~ @ 


P(x, 2 2)= D YD dD | Forn!? 
h=—o k=—o |=—oo 
Xcos 2r(hx+k/2+/z) 
=>) 2 ZX (-1)*| Fuen|? 
h=—o k=—o /=—@ 


Xcos 2r(hx+/z) 


=> s Cy, cos 2r(hx+/z), 


h=—o l=—oo 


oo 


Cri= z. (—1)*| Foy |?. 


k=—oo 


where 


It is easily seen what changes in the calculation 
must be made for threefold, fourfold and sixfold 
screw axes. 

If the crystal contains a plane of symmetry, 
the form to which P(x, y,z) reduces is even 
simpler than in the case of an axis of symmetry. 
Suppose the plane is perpendicular to the b axis. 
Then there are equivalent atoms at (x, y, z) and 
(x, 97, z). The vector between these atoms has 
components (0, 2y,0) and the corresponding 
maximum in P(x, y, 2) lies on the b axis. To find 
this maximum it is only necessary to evaluate 
P(0, y, 0). The calculation is as follows: 


fee) 


P(0, y; 0) —_ i } > | Foy | 2 cos 2rky 


h=—w k=—o l=—c 


= > By cos 2rky 


k=—@ 


with B,.= > > | Foary | 2. 


h=—ao l=—ao 


This is a one-dimensional Fourier series and can 
be completely evaluated in a few hours. In the 
case of a glide plane the calculation is equally 
simple. In this case there are equivalent atoms at 
(x, y, 2) and (x, j, +4) (the glide is taken to be 
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along the c axis) and there will be a maximum in 
P(x, y,2) at (0,2y,4). The simplification of 
P(x, y, 2) takes place thus, 


P(O, y, 4)= > B,. cos 2rky 


with B,.= > D> (—1)!| Foren |*. 


h=—a l=—co 


In Table I is given a summary of the forms of 
P(x, y, 2) best suited to the various symmetry 
elements. It is to be noted that at least one of 
these forms is available for the study of any 
crystal that is not triclinic. 

TABLE I. Forms of P(x, y, 2) best suited to the determination 


of atomic coordinates in crystals having various 
symmetry elements. 








SYMMETRY ELEMENT ForM OF P(x, y, 2) 





(a) Axes parallel to b-axis 


(i) 2, 4, 42, 4, 6, 62* P(x, 0, 2) 
3, 3, 63 
(ii) 21, 41, 43, 61, 65 P(x, 4, 2) 
(iii) 31, 62, 64 P(x, 4, 2) 
(b) Planes perpendicular to b-axis 
(i) reflection planes P(O, y, 0) 
(ii) glide plane, glide =}a0 P(4, y, 0) 
(iii) ** sia ‘*  =$e0 P(O, y, }) 
(iv) =4(a0+co) P(4, y, 4) 
(v) =4}(a0+co) P(3, y, 4) 
(vi) =}(3a0+¢0) P(,, 4) 








In cases (a) x and z are determined for each atom. 

In cases (b) vy is determined for each atom. 

* The nomenclature of symmetry axes is taken from the ‘‘Inter- 
nationale Tabellen zur Bestimmung von Kristallstrukturen.”’ 


There are other cases in which the values of 
some of the atomic coordinates in a crystal can be 
determined by the evaluation of P(x, y, z) 
throughout a plane or along a line only. For 
instance, it might be known from symmetry 
considerations that the position (0, 0, 0) is occu- 
pied by an atom and that there are other atoms in 
the plane y=} with coordinates (say) (x1, 3, 21), 
(x2, 3, 22), etc. Then there will be maxima in 
P(x, y, 2) at just these points, for vectors from 
the origin to these atoms are also interatomic 
vectors. The evaluation of P(x,4,z) would 
determine the: x and z coordinates of all these 
atoms. Or, for another example, suppose again 
that (0, 0,0) is occupied by an atom and that 
other atoms must lie on the line x=}, z=}. The y 
coordinates of all these atoms could be found by 
calculating P(3, y, 3), for there will be a maxi- 
mum in this function at each value of y occupied 
by one of these atoms. Such examples could be 
multiplied almost indefinitely. 
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It is plain from the discussion above that the 
function P(x, y, 2) furnishes a powerful and 
simple method of attack on the structures of all 
crystals not belonging to the triclinic class. All 
that is required is a complete set of x-ray diffrac- 
tion data, for from these not only can the sym- 
metry of the crystal be determined, but also the 
values of | Fyx2 |? calculated. 

An illustration of the use of P(x, y,2z) is 
furnished by the determination of the structure 
of the mineral proustite, Ag3AsS3, which will be 
described below. 


DETERMINATION OF THE STRUCTURE OF 
PROUSTITE, Ags3AsSs3 


The proustite used in this determination came 
from Schneeberg in Saxony. It occurred in the 
form of tiny, brilliant red, hexagonal prisms 
measuring about 0.2 mm by 1.0 mm which were 
bounded by beautiful faces. 

Groth? describes proustite as ditrigonal pyram- 
idal (C3,) with the axial ratio c/a=0.8038. The 
density varies among different specimens from 
5.55 to 5.62. This assignment to a polar class is 
made on the basis of many observations of 
unsymmetrical face development and _ seems 
exceptionally well grounded. 

A Laue photograph taken perpendicular to a 
prism face showed a twofold axis of symmetry, 
another Laue photograph taken along a line 
bisecting the angle between two prism faces and 
perpendicular to the prism axis showed a plane of 
symmetry parallel to the prism axis, and a third 
Laue photograph taken along the prism axis 
showed a threefold axis and three planes of 
symmetry. The Laue symmetry is therefore Dx 
and the crystal must belong to one of the point 
groups C3,, D3, or Dzq on the basis of x-ray data 
alone. 

Measurements on oscillation photographs led 
to the dimensions a)=10.74A, co=8.64A for the 
hexagonal unit cell. The ratio co/a)=0.804 is in 
excellent agreement with the ratio c/a=0.8038 
obtained from crystallographic measurements. 

The two Laue photographs taken perpendicular 
to the prism axis were indexed by the aid of 
gnomonic projections and the wave-length of the 
radiation producing each reflection calculated. 


2 P. Groth, Chemische Kristallographie, Vol. II, p. 767. 
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It was then found that no planes reflected in the 
first order unless the indices satisfied the con- 
dition imposed by the rhombohedral lattice, 
I7H+K+L=0 mod 3, although many planes not 
satisfying this condition were in position to give 
first-order reflections. The lattice is consequently 
rhombohedral and its unit cell has one-third the 
volume of the hexagonal cell for which dimensions 
are recorded above. It is impossible to account 
for much of the observed data on the basis of a 
unit cell smaller than this and no data requiring a 
larger one has been found, although several 
hundred reflections on various photographs have 
been indexed. This rhombohedral cell is therefore 
the true unit cell of proustite. For convenience, 
however, the hexagonal cell will be used through- 
out this paper, much of the analytic geometry 
of the structure appearing simpler when referred 
to these axes. The limits of density given above 
require the number of formulas of Ag3AsS; in 
the hexagonal unit cell to lie between 5.89 and 
5.96. There are thus 6Ag;3AsS; in this cell. 


DETERMINATION OF THE SPACE GROUP 


No reflections of the type (OK-L), L odd, 
were observed in the first order on any photo- 
graph, although many planes of this type were in 
position to reflect on both the Laue and oscil- 
lation photographs. There are two space groups 
based on the rhombohedral lattice which require 
this absence: C3,°—R3c and Dsa—R3c. The 
hemihedral face development of proustite elimi- 
nates the latter. Consequently the space group of 
proustite is C3,°—R3c. 


DETERMINATION OF THE STRUCTURE 


The space group C;,°—R3c furnishes the 
following sets of equivalent positions : 


(2a) 0,0,Z; 0,0, 2Z+4; 
and (6b) X, Y,Z; ¥,X-Y,Z; -—X+4+Y,X,Z; 
X,X-Y,Z+43; Y,X,Z+3; -X+Y,Y,Z+4+}; 


and positiorts derived from these by the opera- 
tions of the rhombohedral lattice. 

There are 2 As, 6 Ag and 6S atoms to distribute 
among these positions. The arsenic atoms must 
be placed in (2a) and the parameter Za, may be 
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arbitrarily taken to be zero. The arsenic atoms 
are now fixed at (0,0,0) and (0, 0,4) in the 
hexagonal cell and the positions derived from 
these by the operations of the rhombohedral 
lattice. The silver and sulfur atoms must be in 
the general positions (6a) and there are conse- 
quently six parameters to be determined: 
aim Tie Las, and Xs, Ys, Zs. 

In order to determine these parameters, the 
Patterson method described in the previous 
section was applied. The crystal contains a glide 
plane and so it is necessary to evaluate the 
Patterson function, P(X, Y, Z), for points along 
a line only, in order to determine the distances 
of the atoms from the plane. By subtracting the 
coordinates of points produced from each other 
by the action of the glide plane perpendicular to 
the Y axis, it was found that there would be 
maxima in P(X, Y, Z) at (0, +(X-—2Y), 3); 
(0, +(—2X+Y), 3); (0, +(X+ Y), 4); for each 
kind of atom in the crystal. It was, therefore, 
only necessary to evaluate P(0, Y, 3) in order to 
determine X and Y for each atom. Now 


pm | Fux.1|? 


—o L=—o 


Ms 


PO, Y,3)= Lb 
H=-—0© 


=o K 


Xcos 2x(K VY+L/2) 


= - Cx cos 2rK Y 


K=-© 


with Cx= > a (—1)4| Fux.z|?. 
H 


=—o L=—o 


The coefficients Cx were calculated from the 
intensities of reflections on a series of oscillation 
photographs. A single small crystal was used 
throughout the series and all photographs were 
given the same exposure. The orientations of the 
crystal were so chosen that each possible reflec- 
tion should have appeared near the equator on at 
least one of the series. The intensities of reflec- 
tions were estimated by comparison with a series 
of spots of known exposure time on a strip of film, 
numbers proportional to the exposure times of 
the comparison spots being used to describe the 
intensities of the reflections. A complete search 
for reflections out to (sin @)/A=0.570 was made, 
sixty-eight forms being observed to reflect within 
this limit. The intensities were divided by the 
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Lorentz and polarization factors in order to 
obtain the values of | Fyx.,|?. No corrections for 
extinction or absorption were made, as there was 
no way of estimating their effects. 

The curve of P(0, Y, 4) plotted against Y 
appears in Fig. 1. The large maximum at Y=0 is 
due, of course, to the vector between the arsenic 
atoms at (0, 0,0) and (0, 0, 3). The maxima at 
&, 3 and 34 were considered to be due mainly to 
the vectors between silver atoms, for the maxima 
due to vectors between sulfur atoms would be 
expected to be only about one tenth the size of 
these. 

The axes can be so chosen that the inequalities 
Xagt Vag>2Vag—Xag>2Xag— Vag>0 are satis- 
fied. This choice and the positions of the maxima 
in P(0, Y, 3) lead to the equations 


Xagt Vag=3; 2Vag—Xag=3; 2X ag— Vag=i- 


These are solved by Xa,=4/18=0.22 and 
Yag=5/18 20.28. The silver atoms are now 
approximately located in the X, Y plane. The 
maxima in P(0, Y, 3) are fairly broad and conse- 
quently the values of Xa, and Ya, determined 
from these are not to be considered as accurate. 
However, a silver atom must lie near the line 
X =4/18, Y=5/18. 
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Fic. 1. P(O, Y, 3) for proustite plotted against Y. The 
large maximum at Y=0 is due to As-As interaction, the 
maxima at 0.17, 0.33, and 0.50 are due to Ag-Ag inter- 
actions. There is a center of symmetry in P(0, Y, 3) at 
Y =}, consequently only half the curve is plotted. 





HARKER 


There is an arsenic atom at (0, 0, 0) and there 
must therefore be a maximum in P(X, Y, Z) at 
(Xag, Yag, Zag). This maximum must lie near 
enough to the line X¥ =4/18, Y=5/18 so that 
P(4/18, 5/18, Z) should show a maximum at 
Z=Za,. P(4/18, 5/18, Z) was accordingly evalu- 
ated and the values plotted in Fig. 2 were 
obtained. The two sharp maxima at Z =0.227 and 
Z=0.743 can be due only to the vectors from the 
arsenic atoms at (0, 0,0) and (0,0, —3) to the 
silver atom. The values of Za, calculated from 
the positions of these two maxima are Za,=0.227 
and Z,,=0.243. The mean of these is Z4,=0.235 
with a probable error of about 0.008. The relation 
between the arsenic and silver atoms in the 
proustite structure is now approximately de- 
termined. 

Hofmann’ has stated that the AsS;’” and 
SbS;’” groups are pyramidal with bond angles of 
about 95° and bond lengths about equal to the 
covalent radius sums calculated from the table of 
Pauling and Huggins.‘ Corey and Wyckoff’ in 
the study of [Ag(NHs)2 ]2SO, and West® in the 
study of AgCN have shown that Ag’ forms two 
covalent bongs in opposed directions with a 
covalent radius of 1.36A. The assumption was 
accordingly made that the AsS;’” groups in 
proustite were of approximately the form de- 
scribed by Hofmann and that silver atoms formed 
two covalent bonds to sulfur in opposed direc- 
tions with a silver-sulfur distance of about 2.40A, 
which is the sum of the covalent radii. It was also 
assumed that sulfur would form bonds in 
approximately tetrahedral directions. It was 
found that these conditions could all be satisfied 
approximately by placing the sulfur atoms in 
proustite in the general positions (6a) with 
Xs20.22, Ys0.11, Zs=0.40. The complete 
structure of proustite was then roughly known. 

The parameters Xa, and Ya, were next fixed 
accurately by means of the intensities of the 
reflections (HH-0), H=1, 2, 3, 4, 5, 6; assuming 
the sulfur positions as above. These intensities 
are as follows: 

H ;e8 das é 
Observed Intensity of (HH-0) 0.50 0.40 0.45 0.10 0.15 0.20 


3 W. Hofmann, Zeits. f. Krist. 92, 161 (1935). : 
' 4 - Pauling and M. L. Huggins, Zeits. f. Krist. 87, 205 
1934). 
5 R. B. Corey and R. W. C. Wyckoff, Zeits. f. Krist. 87, 
264 (1934). 
6 C. D. West, Zeits. f. Krist. 90, 555 (1935). 
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Fic. 2. P(4/18, 5/18, Z) for proustite plotted against Z. The two large maxima are due to 
As-Ag interactions. 


Fig. 3 shows the values of Xa, and Ya, ruled out 
by the various comparisons. The center of the 
allowed region is at X,,=0.246, Ys, =0.298, and 
these values give calculated intensities in good 
agreement with the quantitative data with the 
sulfur parameters anywhere near Xs=2/9, 
Ys=1/9. A change in Xa, or Vay of 0.005 in 
either direction places the point Xa,, Ya, at, or 
beyond, the edge of the allowed region as shown 
in the figure. The limits of error in the determi- 
nation of Xa, and Ya, are consequently +0.005. 
If the Patterson method determination of Za, is 
accepted, the silver parameters are now com- 
pletely determined. The silver positions were 
accordingly assumed to be accurately known and 
the sulfur atoms placed so as to be 2.25A from 
arsenic atoms and 2.40A from silver atoms, these 
being the distances calculated from the covalent 
radii of these elements. The structure of proustite 


is now completely determined. The parameters 
are as follows: 


Xap =0.246+0.005, 
Ya" =0.298 +0.005, 
Zag =0.235 +0.008, 
Zas =0.000. 


Xg =0.220+0.10, 
Ys =0.095 +0.10, 
Zs =0.385 +0.10, 


Intensities of all reflections for which (sin @)/d 
=0.375 were calculated and compared with the 
observed data. The results of this calculation 
appear in Table II. The agreement between the 
observed and calculated values proves the 
structure of proustite described above to be 
correct. 

The structure of proustite being known, it is 
possible to discuss the positions at which maxima 
should fall in P(O, Y,3) and P(4/18, 5/18, Z). 
Suifur to sulfur distances should give maxima in 
P(0O, Y, 3) at Y=0.030, 0.315 and 0.345 corre- 
sponding to Xs—2Ys, Xst+ Vs and 2Xs— Ys, 


TABLE II. Observed and calculated intensities for proustite. In each square the top row gives the indices, and in the second 
row the first figure is the observed, the second figure is the calculated intensity. 








(12 -5) (31-5) 


0.9—0.87 0.45 —1.02 
(10 -4) (02 -4) (21-4) (13 +4) (40 -4) 
0.7-1.93  1.3—1.00 1.3—0.71 1.2—1.55 
(11-3) (22 +3) 
4.0 —4.39 0.25 —0.17 
P (01 -2) (20 -2) (12 -2) (31-2) (04 -2) 
40-0.16 4.5—-2.86  2.5—1.97 1.0—0.60 0.30-—0.20 1.3—1.40 


(21-1) (13-1) (32 -1) 
0.8 —0.77 


1.5—1.51 


(11-0) (03 -0) (22 -0) 
0.50-0.35 1.8—1.71 0.40—0.10 


2.0 —1.79 





(41 -0) 
0.40 —0.22 


(23 -5) 
0.15 —0.45 


(05 -4) (24 -4) 


(32 -4) 
0.9-1.48 0.10-—0.22 0.10—0.42 0.15-—-0.49 


(33 +3) (25 -3) (52 -3) 


(14 -3) (41 -3) 
0.30-—0.13 0.25-0.15 0.25-—0.05 0.15—0.37 0.15—0.29 
(23 +2) (50 -2) (42 +2) (15 +2) (34 +2) (61-2) 


1.2 —1.22 0.45-—0.62 0.40—0.33 0.10—0.14 0.15—0.14 


(24-1) (51-1) (43 +1) (16-1) 
0.40-—0.53 0.00-—0.07 0.30-—0.37  0.15—0.12 


(33 -0) (06 -0) (52 -0) (44 -0) 
0.45—0.24 0.60—1.20 0.40—0.54 0.10—0.04 
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Fic. 3. (left) Determination of Xagand Yag for proustite 
using comparisons of intensities of reflections (HH-0). 


Region A is ruled out by (55-0) calc. > J(66-0) 
ae B “ce a “ee “a (22-0) “é >I(33-0) 
“a C 4c ae oe ae (44-0) “ >I(22-0) 
“é D “a “a “a ae I(44-0) ae >J(55-0) 


X marks the point giving the best quantitative agreement. 


Fic. 4. (right) Determination of Xag and Yag for 
pyrargyrite by comparisons of intensities of reflections 


(HH-0) 


Region A is ruled out by (66-0) calc. < }7(55-0) 
oe B “ee “a aa oc (22-0) “oe > (33-0) 
4 Cc “oe ae “ce “a (44-0) “oe > 17(22-0) 
“oc D “ “oe “ce “ (66-0) ““ > (55-0) 


X marks the point giving the best quantitative agreement. 


respectively, silver to silver distances should give 
maxima at Y=0.194, 0.350, 0.456 and 0.544 
corresponding to 2Xag— Vag, 2Vag—Xag, 1—-Xag 
— Va, and Xag+ Va, respectively, and arsenic- 
arsenic distances account fora maximum at Y=0. 
No other interatomic vectors should give maxima 
in P(0O, Y, $). Comparing these figures with the 
curve for P(0, Y, 4) in Fig. 1, it is seen that the 
maximum at Y=0 must be due toan As—As anda 
S-S distance, that at Y=0.17 must be due toa 
Ag—Ag distance, that at Y=0.33 must be due to 
two S-S and one Ag—Ag distance, and that at 
Y=0.50 to two Ag—Ag distances. The small max- 
ima at Y=0.25 and Y=0.41 are not accounted 
for by any interatomic vectors and must conse- 
quently be due to the background fluctuation 
which occurs in any Fourier series treatment 
using only a finite number of terms. The small 
maxima in P(4/18, 5/18, Z) (see Fig. 2) must be 
due to this same background effect as only the 
Ag-—As vectors should produce maxima in this 
function. 

THE CRYSTAL STRUCTURE OF PYRARGYRITE 
Ag;3SbS3 


The pyrargyrite used in this study was from 
Freiberg in Saxony and occurred in the form of 
metallic appearing, hexagonal prisms measuring 
about 0.5 mm by 2.0 mm. 


Groth? assigned pyrargyrite to the ditrigonal 
pyramidal class (C3,) with the axial ratio 
c/a=0.7892. The density ranges from 5.75 to 
5.85. In this case, as in the case of proustite, the 
assignment to the polar class C3, is well grounded 
on observations of unsymmetrical face de- 
velopment. 

Measurements of reflections on oscillation 
photographs led to the axial lengths a9)=11.04A, 
co=8.71A for the hexagonal unit cell. The axial 
ratio is Co/ay9=0.789, in good agreement with the 
ratio c/a=0.7892 calculated from crystallo- 
graphic measurements. The number of formulas 
of Ag;SbS; in this unit cell was calculated from 
the density range given above to lie between 5.92 
and 6.03. There are thus 6 formulas of Ag:SbS; 
in the unit cell. 

Both the Laue and oscillation photographs of 
pyrargyrite show great similarity to the corre- 
sponding photographs of proustite. In particular, 
the conditions required by the rhombohedral 
lattice and the space group C3,°— R3c were never 
violated. These observations, together with the 





1.0 


hace 














0.25 0.20 
Zas 


Fic. 5. Determination of Zag for pyrargyrite by com- 
parisons of intensities. Jyx-.1, plotted on an arbitrary 
scale against Zag. The dotted line shows the chosen value. 
The observed intensities satisfy the inequalities: /(51-4) 
>1(35-4) >1(27-4) > (26-2) > 1(43-4). 
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CRYSTAL STRUCTURES 
similarity between the chemical compositions, 
crystallographic data, and unit cell dimensions of 
pyrargyrite and proustite, led to the assumption 
that these two substances have nearly the same 
structure, the chief difference being the substi- 
tution of the arsenic atoms in proustite by 
antimony atoms in pyrargyrite. That this as- 
sumption was justified is amply proved by the 
agreement between the calculated and observed 
intensities of reflection which resulted from its 
adoption. 

Gossner and Mussgnug’ made layer line 
measurements on oscillation photographs of 
pyrargyrite from which they determined the 
lengths of the hexagonal axes to be a9=11.06A 
and co= 8.84A in good agreement with the results 
obtained in this research. They also proved the 
lattice to be rhombohedral and determined the 
space group to be C3;,°—R3c. They did not 
determine the structure, however. 

The assumption that proustite and pyrargyrite 
have almost the same structure leads at once to 
the approximate structure for pyrargyrite: 


C3,.°; 2Sb in 2a, Zs, =0.000; 6 S in 66 Xs £0.22, 


6 Ag in 6b, Xa 20.25, Ys =0.10, 
Ya, 20.30, Zs =0.40. 
Dis =0.23 P 


The silver parameters Xa, and Ya, were 
determined, as in the case of proustite, by use of 
the intensities of the reflections (7H-0), H=1 to 
6. These intensities were as follows: 


H 123 4 5 6 
Observed Intensity of (HH-0) 0.0 5.0 6.0 0.40 1.5 0.80 
The regions of variation of Xa, and Ya, ruled 


7B. Gossner and F. Mussgnug, Centralbl. f. Min. (A), 
65 (1928). 
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out by the various comparison are shown in 
Fig. 4. Values of the parameters differing from 
Xag=0.250, Yag=0.305 by as much as 0.005 
gave poor agreement with the observed data, and 
these values were therefore accepted as being 
correct to +0.005. The parameter Za, was 
determined by comparisons, some of which are 
shown in Fig. 5, to be within 0.005 of Z,,=0.210. 
The sulfur atoms were then placed so as to make 
the SbS;’” groups pyramidal with the Sb-S 
distance 2.45A and the Ag-S distance 2.40A, in 
agreement with the covalent radii of these 
elements. The determination of the structure of 
pyrargyrite was then complete. The parameters 
are 


Xag=0.250+0.005, 
Yao =0.305 +0.005, 
Zag =0.210+0.005, 
Zsp = 0.000. 


Xs =0.220+0.010, 
Ys =0.105 +0.010, 
Zs =0.355 +0.010, 


Intensities were calculated on the basis of this 
structure for all the reflections for which calcu- 
lations had been made for proustite. 

The results of these calculations appear in 
Table III. The agreement between the observed 
and calculated intensities of reflection proves the 
above structure to be correct. 
DISCUSSION OF THE STRUCTURES OF PROUSTITE 
AND PYRARGYRITE 


The structure of proustite is shown in Fig. 6, 
and the structure of pyrargyrite is so similar that 
this figure may-be used for the discussion of both. 
Both structures contain (AgS)., groups of almost 
the same shape. In these, the sulfur atoms form 
the corners of a triangular spiral, while the silver 
atoms lie almost in the centers of the legs of the 
spiral. The sulfur to silver distance is 2.40A, the 


TABLE III. Observed and calculated intensities for pyrargyrite. In each square the top row gives the indices, and in the second 
row the first figure is the observed, the second figure is the calculated intensity. 





= 





(12 -5) (31-5) (23 -5) 
2.0 —2.94 1.2 —3.59 0.45 — 1.00 
(10 +4) (02 -4) (21-4) (13 -4) (40 ~y (32 -4) (05 -4) (24-4) 
10.0-12.50 3.0—1.09 6.0 —5.56 7.0 —6.04 4.5 —5.6 0.00—-0.40 0.00-—0.42 2.0—2.82 


(11-3) 
12.0 — 15.37 


(22 +3) (14 -3) (41 +3) (33 +3) (25 +3) (52 +3) 


0.9 — 1.08 0.45—0.40 0.45—0.43 0.50—0.44 0.60—0.68 0.45 —0.69 
(01 -2) (20 -2) (12 +2) (31-2) (04 Fy a (23 +2) (50-2) (42 -2) (15 +2) (34 +2) (61-2) 
6.0-3.60 12.0—12.28 10.0—10.50 4.0—3.19 1.5 —1.09 5.0 —4.73 5.5 —5.23 3.0 —2.49 2.5—2.00 0.50—0.27 1.0—0.91 
(21-1) (13 1) (32 -1) A at 1) (51-1) (43-1) (16 -1) 
5.0 —5.79 5.5 —6.66 2.0—1.61 1.5—1.82 0.00—0.19 1.3 —1.29 0.7 —0.80 
(11-0) (03 -0) (22 -0) (41-0) (33 -0) (06 -0) (52 Te (44-0) 
0.00-—0.13 8.0—10.0 3.0 — 1.39 1.5 —0.92 3.5 —1.75 3.0 —3.62 2.5 —3.1 0.40 —0.16 
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Fic. 6. The structure of proustite. Pyramids represent AsS;’”’ groups, balls represent silver atoms, and 
lines represent covalent bonds. The structure of pyrargyrite is almost the same. 


S—Ag-—S bond angle is about 165°, and the 
Ag—S—Ag bond angle is about 833°. One 
complete turn of the spiral is equivalent to a 
translation of co in the crystal. These spirals are 
fastened together by the As or Sb atoms to make 
up the structural frameworks of proustite or 
pyrargyrite, respectively. The sulfur to arsenic 
bond distance in proustite is 2.25A and the 
S—As-—S bond angle is 103°. The sulfur atom in 
proustite forms three bonds, two to silver and 
one to arsenic with these angles: As—S—Ag 
=101° and 1053°, Ag—S—Ag =833°. In pyrar- 
gyrite the antimony atoms form three bonds to 
sulfur 2.45A long at mutual angles of 953°. The 
sulfur atoms in pyrargyrite form three bonds 
with the angles: Sbh—S—Ag=97° and 111°, 
Ag—S—Ag=83}°. 

An interesting feature of these structures is the 
fact that they are composed of two _ inter- 
penetrating frameworks of the type described in 
the last paragraph which are apparently not 


connected by any bonds at all. The relation of 
the two frameworks can best be understood by a 
study of Fig. 6. The As or Sb atoms at (0, 0, 0) 
belong to one framework and the crystallographi- 
cally equivalent As or Sb atoms at (0, 0, 3) belong 
to the other. The two frameworks are in the 
relation to each other required by the opera- 
tion of the glide planes. One framework contains 
only right-handed (AgS),, spirals, the other only 
left-handed ones. Another example of two 
interpenetrating frameworks apparently not 
interconnected by chemical bonds is furnished by 
the structure of cuprite, Cu,O. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


Water Vapor Discharge and Hydrogen Peroxide Formation 


The discharge (with or without electrodes) through 
water vapor may produce a mixture of the following mole- 
cules and atoms: 


H, OH, Hz, O, O2, H2O (H2O2, HO», O2). 


Some years ago Geib and Harteck! showed that between 
30°K to 150°K and at pressures of 0.1-0.5 mm Hg hydrogen 
atoms and molecular oxygen react on the walls without 
activation energy with 100 percent efficiency to give pure 
hydrogen peroxide. If however the temperature is raised 
to 190°K, rather poor yields of H2O:2 are obtained. Thus 
we were astonished, when Rodebush and Wahl suggested 
that the formation of hydrogen peroxide in the electrode- 
less water vapor discharge was due to a bimolecular associ- 
ation of OH radicals. Recently I carried out experiments 
on this and using an electrode discharge I got the following 
results: 

Large quantities of HO» are collected in a trap cooled 
with liquid air both in the water vapor discharge and in 
the discharge through H», O2 mixture. On the other hand, 
if the products are frozen out at —80°C both in the dis- 
charge tube itself and 10 cm from it, not a singie drop of 
N/10 KMn0O, solution is needed to oxidize the hydrogen 
peroxide. 

These results now are an excellent agreement with those 
of Rodebush and Campbell,? who used an electrodeless 
discharge. 

As the result of all experiments it must be stated: (1) It 
is quite unimportant whether or not electrodes are used; 
(2) Hydrogen peroxide is formed .by the products of the 
water vapor discharge only at temperatures below — 100°C; 
(3) Hydrogen peroxide is formed by a mixture of H atoms 
and molecular oxygen and a mixture of H atoms and O 
atoms at temperatures below —100°C; (4) H atoms are 
destroyed by adsorption and recombination at a glass 
surface cooled with liquid air, but not with solid carbon 
dioxide.’ 

This suggests the following mechanism* of hydrogen 
peroxide-formation 


W+H2WH, We=vwall (glass surface, HO or H.0:); 


WH +0:—WHO;; 


WHO, +H—~W-+H.0.. 
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not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


It seems rather difficult to show, that beside this process 
the recombination of OH radicals to HO, plays an 
important role (perhaps also as wall reaction) in the 
formation of HO», as Rodebush and Wahl‘ suggest. Ex- 
periment only shows that no such formation can be 
detected at temperatures above — 100°C at low pressures. 

In this connection one should discuss again both the 
reaction 


H+H,0.—H,0+0H? 
and the Bonhoeffer and Pearson mechanjsm® 
20H—H,0+0, 


assigning to the latter some activation energy because of 
the work of Oldenberg.® 

Rodebush and Wahl don’t find O atoms’ in their experi- 
ments. Possibly they react very fast, perhaps at every 
10th collision, with OH to O.+H. 

K. H. Gers 
Physikalisch-Chemisches Institut der Universitat Leipzig 
April 22, 1936. 


1 Geib and Harteck, Ber. d. D. Chem. Ges. 65, 1551 (1932); Zeits. f. 
physik. Chemie A170, 1 (1934). 

2 Campbell and Rodebush, J. Chem. Phys. 4, 293 (1936). 

3 Geib, Zeits. f. physik. Chemie A169, 161 (1934). 

4 Rodebush and Wahl, J. Chem. Phys. 1, 696 (1933). 

5 Bonhoeffer and Pearson, Zeits. f. physik. Chemie B14, 5 (1931). 

6 Oldenberg, J. Chem. Phys. 3, 266 (1935). 

7 See Oldenberg, reference 6, p. 273. 


Entropy of Organic Compounds from Calorimetric Data. 
Lack of Equilibrium in Crystalline Tetramethylmethane 


Recently Kassel! has calculated the entropy of n-butane 
from spectroscopic data and found a value 3.76 e.u. higher 
than that of Huffman, Parks and Barmore? which he 
suggests may be due to lack of equilibrium in the solid at 
low temperatures or to an error in the value extrapolated 
below 90°K. As the total extrapolated value is only 11.7 
e.u., it is hardly possible that the error could be due to this 
cause unless a transition of rather large entropy change 
were involved. 

We have recently measured the entropy difference 
between ordinary crystalline tetramethylmethane at the 
absolute zero and the gas at 25°C from calorimetric meas- 
urements carried down to 14°K, the extrapolated entropy 
being only 0.66 e.u. While the calculations are not entirely 
complete a comparison of a preliminary value of this 
calorimetric entropy difference with the entropy of the 
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gas at 25°C shows the calorimetric value to be about 8 e.u. 
lower than the spectropscopic value. This indicates about 
8 e.u. zero-point entropy in the ordinary crystal of tetra- 
methylmethane at the absolute zero, due to randomness 
caused by lack of equilibrium. There is an abnormal rise 
in the heat capacity of the solid prior to a transition at 
140°K. In this region only is there evidence of hysteresis. 

Giauque and associates have found a similar state of 
affairs for carbon monoxide** and nitrous oxide*® and 
predicted the occurrence of the phenomenon in the case of 
organic compounds. Pauling has discussed the situation 
for some simple molecules.‘ 

In view of the large discrepancy in the case of tetra- 
methylmethane one must view entropies of organic com- 
pounds calculated from calorimetric data assuming a 
perfect crystal (i.e., zero entropy at the absolute zero) 
with considerable doubt even though they be based on 
heat capacities taken down to liquid hydrogen tempera- 
tures or lower. 

At present the only source of a reliable value or the 
entropy of neopentane is a calculation based on spectro- 
scopic data and a knowledge of the moments of inertia. 
Such a calculation yields a value of 79.85 e.u. at 25° and 
1 atmos. neglecting nuclear spin. This value is the trans- 
lational and rotational entropy 73.98 e.u.! plus the entropy 
of vibration (5.87 e.u.) calculated from the Raman fre- 
quencies determined by Rank and Bordner,' using the fol- 
lowing multiplicities: 335, (2); 414, (3); 733, (1); 925, (3). 

The complete calorimetric data will be submitted for 
publication in the near future. 

Joun G. Aston 
GeEoRGE H. MESSERLY 


School of Chemistry and Physics, 
The Pennsylvania State College, 
State College, Pennsylvania, 
April 23, 1936. 
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3» Blue and Giauque, ibid. 57, 991 (1935). 

4 Pauling, J. Am. Chem. Soc. 57, 2680 (1935). 

5 Rank and Bordner, J. Chem. Phys. 3, 248 (1935). 


The Raman Effect and the Dissociation of KHSO, 


It has been found by Rao! and others that a study of 
the Raman effect offers a means of determining the degree 
of dissociation of concentrated solutions, where other 
methods are not applicable, and thus supplies valuable 
information with respect to the physical and chemical 
properties of such solutions. Woodward,? Nisi,* and Bell 
and Jeppesen‘ have observed the Raman effect in KHSO, 
at one concentration but made no observations on the 
change in the degree of dissociation with the concentration. 


They found that HSO, ions occur in a concentrated 
aqueous solution of KHSO, and that these ions partially 


+ ao 
dissociate into H and SQ,. This is precisely what might 
be expected from chemical considerations which indicate 
that KHSO, would dissociate in two stages: 


+ ~ - + 
KHSO,—K+HSO,; HSO,—H+S0Q,. 
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Since other studies of the Raman effect in concentrated 
solutions indicate that normal salts are completely dis- 
sociated at all concentrations, and that in strong acids 
the degree of dissociation increases with dilution, it is of 
interest to consider what happens in the Raman spectrum 
of an acid salt like KHSO, when the concentration of the 
solution is changed. 

By means of a high speed spectrograph of low dispersion 
the Raman effect was photographed in four aqueous 
solutions of KHSO, of the following concentrations: 13, 
21, 28, 34 percent. The solutions were enclosed in the 
usual type of Wood tube and illuminated by mercury are 
lamps. The Raman spectra were excited by the 4358A and 
4047A lines. The frequencies of the observed Raman lines 
expressed in wave numbers (cm~') for concentrations of 
13 percent and 34 percent have been recorded in Table I. 


TABLE I. 
CONCENTRATION HSO, HSO, SOs ao —-- 
13% 1058 917 978 611 438 
34% 1040 906 975 599 436 


The ions to which these observed Raman lines have been 
attributed are indicated by the headings of the columns in 
the table. There is a change in frequency of the Raman 
lines with the dilution of the solution. A similar effect has 
been observed by Woodward and Horner® and also by 
Bell and Jeppesen‘ in the case of aqueous solutions of 
H.SQ,. The reality of such a change in frequency due to a 
change in concentration has been questioned by Mitra‘ 
but his observatjons on aqueous solutions of magnesium 
sulphate and lithium sulphate were concerned only with 


the Raman line 981 cm™ which has been attributed to SQ,. 
For this line the change in frequency due to the change in 
the concentration of the solution is small compared to the 
changes observed for the other lines. 

From a comparison of the intensities of the Raman lines 
as indicated by the microphotometer curves, the relative 
degrees of dissociation of the solutions with respect to the 
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degree of dissociation of 13 percent solution have been 
calculated and recorded in Table II. These results are also 


TABLE Ii. 


PERCENT CONCENTRATION 13 21 28 34 
RELATIVE DEGREE OF DISSOCIATION 1.00 0.84 0.68 0.82 


shown as a curve (Fig. 1) with the degree of dissociation 
plotted against the concentration. From this curve the 


progressive dissociation of the HSO, ion is clearly seen. 
No lines attributable to the undissociated KHSO, molecule 
were observed even in the most concentrated solution. 
The fact that no undissociated molecules are present in 
solution is characteristic of normal salts, while the fact 


that the HSQ, ion dissociates progressively is characteristic 
of acids. Thus the Raman method of investigation indi- 
cates that KHSO, has characteristics common to both acids 
and normal salts as would be inferred from chemical 
considerations. 

W. H. SHAFFER 

D. M. CAMERON 

Mendenhall Laboratory of Physics, 


Ohio State University, 
April 20, 1936. 
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Quantum Mechanics and the Third Law of 
Thermodynamics 


Rodebush! has recently discussed certain points con- 
cerning the third law of thermodynamics in terms of 
quantum mechanics. His recognition of the significance of 
quantum mechanics in considering the behavior of solu- 
tions and glasses at low temperatures we believe to be of 
interest, and we agree with the more important of his 
conclusions. However, there have arisen some questions of 
interpretation in which we differ from the statements of 
his paper. We have thought it of value to present here, as 
explicitly as possible and at the risk of some duplication of 
ideas, an outline of these points and a discussion of their 
bearing on the expected behavior of the systems under 
treatment. It is our opinion that the views presented below 
have been in the main accepted by investigators interested 
in the theory of the third law; however, so far as we are 
aware no equivalent discussion has been published. 

Without entering into the question of the relation 
between statistical mechanics and thermodynamics, we 
accept as giving a value for the entropy of the system in 
which we are interested the expression & In N, in which N 
is the number of independent wave functions for the system 
compatible with our knowledge as to its condition. Fol- 
lowing the plan adopted by Rodebush, we first discuss 
the case of an optically active substance.? An isolated 
molecule of such a substance in its lowest rotational and 
vibrational state and known to be in the dextro configura- 
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tion’ is to be represented by a single wave function corre- 
sponding to this configuration, and will have the quantum 
weight 1, and not 2. On the other hand, if the molecule is 
in one of the two lowest energy states corresponding to the 
accurate solution of the wave equation (as could be deter- 
mined by an energy measurement extending over a period 
of time long compared to the time required for spontaneous 
dextro-levo conversion), it is to be represented by a single 
wave function which is either symmetric or antisymmetric 
relative to inversion about a center, and which is formed 
by linear combination of the dextro and levo wave func- 
tions. It will then not be known whether the molecule is a 
dextro or a levo molecule. Only in the case that neither 
its optical activity nor its (exact) energy is known can we 
assign to it the quantum weight 2. 

By analogy with the simple case of the optical isomers 
Rodebush generalizes the many-minima problem exempli- 
fied by glasses and crystalline solutions with the following 
statement. ‘‘We have a number of configurations of the 
same energy separated by energy crests. By symmetric 
and antisymmetric linear combinations of wave functions,‘ 
we obtain a total number of wave functions equal to the 
degeneracy; that is to say, the total number of configura- 
tions having the same energy. Thus, instead of each con- 
figuration having a distinct energy, each configuration has 
now a large number of closely spaced energy levels.”’ 

Here we prefer to state the situation in a somewhat 
different way, as follows. The system in its low temperature 
condition and known to have a particular configuration® 
(one of N possible configurations) is to be represented by 
one of an appropriate set of N wave functions, correspond- 
ing to the N configurations. If not the configuration, but 
the energy of the system is to be investigated with great 
accuracy, by means of observations extending over a very 
long period of time, the set of wave functions of interest is 
that corresponding to accurate solution of the wave equa- 
tion (these wave functions being linear combinations of 
those corresponding to the N configurations). The energy 
value found by experiment would be one of the multiplet 


‘of N levels corresponding to the low temperature condition. 


We cannot, however, make use of the two sets of wave 
functions simultaneously, but must use either the con- 
figuration wave functions or the multiplet wave functions, 
or some other equivalent set. In other words, both the con- 
figuration and the exact energy of the system cannot be 
simultaneously known. If either is known a single wave 
function (of the set of NV) is sufficient to describe the system 
and the system will then have the quantum weight 1, and 
not N. However, if neither the configuration nor the exact 
energy of the system is known, we may use either set of V 
wave functions (or any other set formed from them by 
linear combinations) to describe it, and its quantum weight 
will then be N. 

The behavior of systems on cooling to a very low tem- 
perature (kT<hyv, where » is the frequency of the lowest 
oscillational state) may be outlined in the case of an ideal 
solution somewhat as follows.® If »* is the frequency of 
transition among the various configurations (as given by 
the corresponding resonance integral divided by h), the 
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energy spread of the multiplet levels will be of the order 
of hv*. If the cooling be done so slowly that a time very long 
compared to 1/»* is allowed at each temperature, the loss 
of heat will be such as to correspond to a loss of entropy of 
kin N, since the multiplet levels under such conditions 
may be regarded as discrete. On rapid cooling to the same 
low temperature, the energy of the system will be deter- 
minable only to within about h/2a7, where 7 is the time 
of observation at the temperature in question. If this un- 
certainty in energy is large compared to hy* there will still 
remain N states (N independent wave functions) of the 
system falling within the observable energy limits, and the 
entropy will be k In N. Only after a long time of observa- 
tion at this temperature can these limits be narrowed suf- 
ficiently to permit the statement that a single state is 
represented and the entropy is zero. 

A real solution or glass differs from an ideal solution or 
glass, discussed above, in that one of its configurations, 
corresponding to a pure crystalline substance or a mixture 
of pure crystalline substances, possesses an energy lower 
than any other. If such a system be subjected to very slow 
cooling, or be held for a long time at a low temperature as 
discussed above, it will ultimately attain this configuration 
(mainly by transition through excited vibrational states 
rather than by penetration of potential barriers), and will 
therefore have zero entropy. On cooling sufficiently rapidly 
to exclude effectively transitions between configurations, 
or on holding at the low temperature for a time shorter than 
required for these transitions to occur, the multiplicity N 
will be retained by the system, its entropy then remaining 
greater than zero. 

It may be pointed out that the discussion of the entropy 
of a real solution or glass on the basis of quantum mechan- 
ics, aS just given, is identical with the discussion of the 
problem in the pre-quantum-mechanical era. On the other 
hand, the quantum mechanics has introduced a significant 
change in the discussion of the entropy of an ideal solution 
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or glass. Before the development of quantum mechanics 
there was no conceivable mechanism leading to the loss 
of the final & In N of entropy at very low temperatures for 
an ideal solution or glass, defined as having N configura- 
tions with exactly the same energy. The quantum-mechan- 
ical mechanism leading to the loss of this entropy is the 
resonance of the system among the N configurations. 

Rodebush has also implied that the accuracy with which 
very low temperatures can be measured is restricted by the 
uncertainty principle and by the nature of the substance 
under investigation. However, the accuracy of a tempera- 
ture measurement is not limited in a serious way by the 
uncertainty principle for energy, inasmuch as the relation 
between the uncertainty in temperature and the length of 
time involved in the measurement depends on the size of 
the thermometer, and the uncertainty in temperature can 
be made arbitrarily small by sufficiently increasing the size 
of the thermometer; we assume as the temperature of the 
substance the temperature of the surrounding thermostat 
with which it is in either stable or metastable equilibrium, 
provided that thermal equilibrium effective for the time of 
the investigation is reached. 

Linus PAULING 
Gates Chemical Laboratory, 
California Institute of Technology, 


E. D. EASTMAN 


Department of Chemistry, 
University of California, 
March 30, 1936. 


1W. H. Rodebush, J. Chem. Phys. 2, 669 (1934). 

2 For a thorough discussion of this see F. Hund, Zeits. f. Physik 43, 
805 (1927). 

3 We shall use the word configuration in this paper to represent the 
approximate relative spatial arrangement of the nuclei of the atoms in 
the system under consideration (neglecting the effect of the vibrational 
metion on the nuclear positions). 

4 At the request of Professor Rodebush, the wording of this sentence 
has been changed slightly from that in his paper. 

5 We assume each configuration to be in its lowest vibrational state. 

6 We believe that the ideas discussed in this connection parallel those 
¢ — but have thought it worth while to give them in greater 
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